Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,499,391 G→T intergenic (+190/‑57) cysG → / → yhfL siroheme synthase/DUF4223 domain‑containing lipoprotein YhfL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,499,3910GT97.9% 160.3 / ‑5.3 47intergenic (+190/‑57)cysG/yhfLsiroheme synthase/DUF4223 domain‑containing lipoprotein YhfL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (25/21);  minor base A (0/1);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGG  >  NC_000913/3499162‑3499633
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
tgGCCTGCGCGATACACTGAACTGGTTCTCCAACCATTAATTTACACCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCCGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTAACCGTTACAGCGCCTTACAAAATTAAGTCATTCGCTCACTGCCTACAATTAAATCTTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCatcgaccagat                                                                                                                                                                                                                                >  2:342029‑M1/1‑239 (MQ=255)
       cgcgATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgcc                                                                                                                                                                                                                         >  1:35303‑M1/1‑232 (MQ=255)
             aaaCTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcg                                                                                                                                                                                                                   <  2:452870‑M1/250‑25 (MQ=255)
             aaaCTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcg                                                                                                                                                                                                                   <  1:468556‑M1/250‑25 (MQ=255)
                    actgGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggc                                                                                                                                                                                                           <  1:134024‑M1/251‑33 (MQ=255)
                               aaCCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttaaatccgt                                                                                                                                                                                                >  2:192202‑M1/1‑208 (MQ=255)
                                                         aaTTCTGGCCTTCTTAATGACTCTTTTTATCTATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgtttt                                                                                                                                                                      >  2:233216‑M1/1‑182 (MQ=255)
                                                                 ccTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATCTTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaccgcgg                                                                                                                                                              >  2:331378‑M1/1‑174 (MQ=255)
                                                                           gACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtca                                                                                                                                                     <  1:551009‑M1/250‑87 (MQ=255)
                                                                                         atTCAGCAAATAAAACAAGGCCCCTTAATAAATATTATCTGCTAATTGGTTAGATTCTCCGTTCCAGCGACTTGAAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTAGCGTTCTTAACGGGCAGTCGTCGTATATTCctcgaccagataactgccattgcgcccacggcggtcaatccgtatcgctccgcgtgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaac                                                                                                                                      <  2:243967‑M1/251‑102 (MQ=255)
                                                                                                           aGCCCCTTAATAAATATTATCTGCTAATGTGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaccgcggcctgcgtcaggccaacggatgaacggcattcatttgccagct                                                                                                                    <  2:425259‑M1/251‑120 (MQ=255)
                                                                                                                    ataaataTTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatgg                                                                                                           >  2:169155‑M1/1‑123 (MQ=255)
                                                                                                                          attatCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggagttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttc                                                                                                      >  2:139658‑M1/1‑117 (MQ=255)
                                                                                                                                 gCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgatt                                                                                              >  1:477143‑M1/1‑110 (MQ=255)
                                                                                                                                  cTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtctcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaccggatgaacggcattcatttgccagctcatctatggttttcaatctattc                                                                                             >  2:552364‑M1/1‑109 (MQ=255)
                                                                                                                                    aaTGGTTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattac                                                                                            <  1:225013‑M1/250‑144 (MQ=255)
                                                                                                                                       gggTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgca                                                                                        <  1:449261‑M1/251‑148 (MQ=255)
                                                                                                                                        ggTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttaaatcgattactgca                                                                                        >  2:2293‑M1/1‑103 (MQ=255)
                                                                                                                                                  aCCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcaccc                                                                              >  2:199765‑M1/1‑93 (MQ=255)
                                                                                                                                                        aCAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgcctccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaccaaat                                                                       >  2:25652‑M1/1‑87 (MQ=255)
                                                                                                                                                            cGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgctcttttcacccaacaaataacg                                                                   <  1:102098‑M1/251‑169 (MQ=255)
                                                                                                                                                            cGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacg                                                                   >  2:238987‑M1/1‑83 (MQ=255)
                                                                                                                                                                          aaGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcccccaacaaataacgggcca                                                              >  2:371708‑M1/1‑69 (MQ=255)
                                                                                                                                                                          aaGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggcca                                                              <  1:371708‑M1/242‑174 (MQ=255)
                                                                                                                                                                            gTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataaagggccagcgtcgtgatg                                                   >  2:555058‑M1/1‑67 (MQ=255)
                                                                                                                                                                                  cGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacg                                             >  1:138110‑M1/1‑61 (MQ=255)
                                                                                                                                                                                      cACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgcct                                          <  2:138110‑M1/250‑194 (MQ=255)
                                                                                                                                                                                         tGACGCCAATTAAATATTACTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataaatgccattgcgcacacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgacagctcatctatggttttcaatcgattaccgcacttttcacacaacaaataacggtccagcgtcatgaggacgacgacttctt                                      >  2:492099‑M1/1‑54 (MQ=255)
                                                                                                                                                                                          gACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttt                                     <  1:435442‑M1/251‑199 (MQ=255)
                                                                                                                                                                                          gACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggacaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttt                                     >  1:346387‑M1/1‑53 (MQ=255)
                                                                                                                                                                                              ccAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgacaagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccaacgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttg                                 >  2:402310‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                    aaaTATTCCTGTATCGTTCTTAACGAGCAGTATATGTATATGCctctaccagataactgccattgcgtccacggcgttcaatccggatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggtgttcaagcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgtcttcttttttgataaac                           <  2:367163‑M1/251‑209 (MQ=255)
                                                                                                                                                                                                        atTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcgtttttgataaacgtt                        <  1:139658‑M1/250‑212 (MQ=255)
                                                                                                                                                                                                         ttCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctctaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgagtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattacctcacttttcacccaacaaataacggaccagcgtcgtgatgacgacgccttcttttttgataaacgttcg                      >  1:6518‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                          tCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcg                      <  2:304407‑M1/250‑214 (MQ=255)
                                                                                                                                                                                                           ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgac                    <  1:438171‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                           ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataacagacattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcatcgtgatgacgacgccttcttttttgataaacgttcgac                    >  1:576452‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                           ccTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCcccgaccagataactgccattgcgcccacggcggtcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgac                    <  1:102707‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataactggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcata               <  1:494835‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgac                                                 <  1:439503‑M1/217‑187 (MQ=255)
                                                                                                                                                                                                                aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgac                                                 >  2:439503‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                  cGTTCTTACCGAGCCGTCTTCGTATATTCctcgaccagataaccgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttctcccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataa              <  2:124706‑M1/250‑222 (MQ=255)
                                                                                                                                                                                                                  cGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccccggcgttcaatccgtagcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataat             <  2:419986‑M1/251‑223 (MQ=255)
                                                                                                                                                                                                                       ttAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttcataaacgttcgcagcataataaatt        >  2:540152‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                            gAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttccg   >  1:529341‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                             aGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttcagt  <  1:559017‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                              gCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcatactaaattttcagt  >  2:424195‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
TGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGG  >  NC_000913/3499162‑3499633

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: