Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 842,009 (GA)4→3 coding (47‑48/261 nt) mcbA ← DUF1471 domain‑containing protein McbA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913842,0030G.98.8% 362.5 / ‑2.7 82coding (54/261 nt)mcbADUF1471 domain‑containing protein McbA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (30/51);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913842,0040A.98.8% 360.6 / ‑2.6 82coding (53/261 nt)mcbADUF1471 domain‑containing protein McbA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (30/51);  minor base C (1/0);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGA  >  NC_000913/841865‑842136
                                                                                                                                          ||                                                                                                                                    
cgcgCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTga                                                                                                                                      >  1:1316542/1‑140 (MQ=255)
    ccTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagaga                                                                                                                                  >  2:359631/1‑140 (MQ=255)
      ttGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagatt                                                                                                                                >  1:552441/1‑138 (MQ=255)
      ttGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagatt                                                                                                                                <  1:1373259/140‑3 (MQ=255)
      ttGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagatt                                                                                                                                <  2:841130/140‑3 (MQ=255)
      ttGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagatt                                                                                                                                <  1:25303/140‑3 (MQ=255)
           cgcgcgCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagattccgct                                                                                                                           <  1:727325/140‑8 (MQ=255)
           cgcgcgCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTgagagattccgct                                                                                                                           <  1:319578/140‑8 (MQ=255)
               cgcTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAgg                                                                                                                     >  1:777555/1‑140 (MQ=255)
                   ttttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGGTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACgg                                                                                                                 >  2:1412370/1‑140 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAGCGGTACCGGCGGGCCGTAGCTGACCGCTGGCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  2:1348892/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  1:1225722/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  2:1448870/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  1:485942/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  2:467999/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  2:60394/140‑1 (MQ=255)
                     ttCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTg                                                                                                               <  2:1003399/140‑1 (MQ=255)
                       cTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGc                                                                                                             <  1:78717/140‑1 (MQ=255)
                       cTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGc                                                                                                             <  1:900007/140‑1 (MQ=255)
                       cTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCAAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGc                                                                                                             <  1:1496919/140‑1 (MQ=255)
                        tGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCa                                                                                                            <  1:888411/140‑1 (MQ=255)
                        tGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCa                                                                                                            <  1:72624/140‑1 (MQ=255)
                          ccAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAAt                                                                                                          <  1:1020228/140‑1 (MQ=255)
                          ccAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAAt                                                                                                          <  1:357890/140‑1 (MQ=255)
                          ccAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAAt                                                                                                          <  1:677556/140‑1 (MQ=255)
                            aGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCa                                                                                                        >  2:515782/1‑140 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTTCAGGCTACTCATCGGGTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGAGGCAAGCag                                                                                                       >  2:265323/1‑140 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCGTCGGTTGTGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  1:400312/140‑1 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  2:189705/140‑1 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  2:1429444/140‑1 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  1:673266/140‑1 (MQ=255)
                             gTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  1:873433/140‑1 (MQ=255)
                              tttATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCagt                                                                                                      >  1:1334907/1‑140 (MQ=255)
                               ttATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCagta                                                                                                     <  2:607091/140‑1 (MQ=255)
                                       cAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAgg                                                                                             <  2:867058/140‑1 (MQ=255)
                                       cAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAgg                                                                                             <  1:203368/140‑1 (MQ=255)
                                        aGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGc                                                                                            >  2:739470/1‑140 (MQ=255)
                                        aGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGGAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGGGGCAATCAGTAGTGTGAGGc                                                                                            >  1:1348537/1‑140 (MQ=255)
                                              cTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACtttt                                                                                      <  2:1045951/140‑1 (MQ=255)
                                               ttAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACttttt                                                                                     >  1:338387/1‑140 (MQ=255)
                                                 aGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCa                                                                                   >  1:449941/1‑140 (MQ=255)
                                                   gTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATa                                                                                 <  2:1316542/140‑1 (MQ=255)
                                                       gATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCAT‑‑GAGAGATTCCGCTCAGGACGGTGGCTATCAGTAGTGTGAGGCACTTTTTCATAATGa                                                                             >  1:430284/1‑140 (MQ=255)
                                                        aTGCGCCAGTTGCCGAAACGGTACCGGCGGGGCGTAGCTGACCGCTGTCCAGGTAACTCCTCGGTAGCGCGGCGTAAGCCGTGAGGGATTCCGCTCAGGGCGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGAc                                                                            >  2:1101801/1‑140 (MQ=255)
                                                        aTGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGAc                                                                            <  1:628145/140‑1 (MQ=255)
                                                         tGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACa                                                                           <  1:265323/140‑1 (MQ=255)
                                                         tGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACa                                                                           <  1:1397093/140‑1 (MQ=255)
                                                             ccAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGGGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATccc                                                                       >  2:327550/1‑139 (MQ=255)
                                                             ccAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGGGGCAATCAGTAGTGGGAGGCACTTTTTCATAATGACATCCt                                                                       >  2:600658/1‑140 (MQ=255)
                                                                gTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCAt                                                                    <  2:498485/140‑1 (MQ=255)
                                                                    ccGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAAt                                                                <  2:1334907/140‑1 (MQ=255)
                                                                        aaCGGTACCGGGGGGCCGTAGCTGACCGCTGTCCAGGTTTCTCATCGGTTGCGGGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTATTTCATAATGACATCCTCATAAATGCat                                                            <  2:1278920/140‑2 (MQ=255)
                                                                        aaCGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACaa                                                            <  2:74083/140‑1 (MQ=255)
                                                                          cGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACg                                                          >  1:1315755/1‑140 (MQ=255)
                                                                            gTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACggg                                                        <  2:1182700/140‑1 (MQ=255)
                                                                             tACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGt                                                       >  1:91627/1‑140 (MQ=255)
                                                                             tACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCACAAATACAACGGGt                                                       <  2:966716/140‑1 (MQ=255)
                                                                              aCCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTa                                                      >  1:1471972/1‑140 (MQ=255)
                                                                              aCCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAg                                                     >  1:315113/1‑140 (MQ=255)
                                                                               ccGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAg                                                     <  1:18239/140‑1 (MQ=255)
                                                                                cggcggGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGt                                                    <  2:777555/140‑1 (MQ=255)
                                                                                 ggcggGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGc                                                  <  2:315113/140‑1 (MQ=255)
                                                                                  gcggGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGc                                                  >  2:1395612/1‑140 (MQ=255)
                                                                                   cggGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGcc                                                 <  1:934583/140‑1 (MQ=255)
                                                                                     ggCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTg                                               >  2:1168115/1‑140 (MQ=255)
                                                                                      gCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGa                                              <  1:1456363/140‑1 (MQ=255)
                                                                                      gCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGa                                              <  1:307692/140‑1 (MQ=255)
                                                                                       ccGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGAt                                             <  1:1421197/140‑1 (MQ=255)
                                                                                       ccGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGAt                                             <  1:1332277/140‑1 (MQ=255)
                                                                                       ccGTAGCTGACCGCTGTCAAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGAt                                             <  1:916074/140‑1 (MQ=255)
                                                                                        cGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATa                                            <  2:355965/140‑1 (MQ=255)
                                                                                          tAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAg                                          <  2:1020926/140‑1 (MQ=255)
                                                                                             cTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCC‑‑GCGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAgtgt                                       >  1:1442731/1‑140 (MQ=255)
                                                                                               gACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAg                                     <  2:1171618/140‑1 (MQ=255)
                                                                                                aCCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAgg                                    >  2:1280641/1‑140 (MQ=255)
                                                                                                aCCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAgg                                    >  1:1474306/1‑140 (MQ=255)
                                                                                                 ccGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGt                                   <  2:1471972/140‑1 (MQ=255)
                                                                                                  cGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTc                                  >  1:930331/1‑140 (MQ=255)
                                                                                                         cAGGTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTgcg                           <  2:1193340/140‑1 (MQ=255)
                                                                                                           ggTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGcgcg                         <  1:714424/140‑1 (MQ=255)
                                                                                                            gTTACTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGa                        >  2:1092651/1‑140 (MQ=255)
                                                                                                               aCTCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGc                     >  2:138437/1‑140 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       >  2:272283/1‑56 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCag                                                                                                       <  1:272283/56‑1 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCaa                   >  1:815353/1‑140 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCaa                   >  1:1449466/1‑140 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCaa                   >  1:291740/1‑140 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCaa                   >  2:794281/1‑140 (MQ=255)
                                                                                                                 tCATCGGTTGCGCGGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCaa                   >  1:702689/1‑140 (MQ=255)
                                                                                                                           cgcgGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCaaa         >  2:1168509/1‑140 (MQ=255)
                                                                                                                             cgGCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCaaaaa       <  2:1431195/140‑1 (MQ=255)
                                                                                                                               gCATAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAAt     <  2:552441/140‑1 (MQ=255)
                                                                                                                                  tAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGa  >  2:173548/1‑140 (MQ=255)
                                                                                                                                  tAAGCCGT‑‑GAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGa  <  1:515782/140‑1 (MQ=255)
                                                                                                                                          ||                                                                                                                                    
CGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGA  >  NC_000913/841865‑842136

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