Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,399,905 A→T 5.1% I25N (ATT→AAT)  mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,9050AT 5.1% 287.5 / 4.8 79I25N (ATT→AAT) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base T (2/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

TTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTT  >  NC_000913/3399757‑3400051
                                                                                                                                                    |                                                                                                                                                  
ttCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                                    >  1:1279961/1‑149 (MQ=255)
ttCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                                    <  2:1279961/149‑1 (MQ=255)
 tCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                                   <  4:930913/149‑1 (MQ=255)
   tAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCATTATTATTGc                                                                                                                                                 >  2:571719/1‑149 (MQ=255)
      ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCATTATTATTGCCGa                                                                                                                                              <  1:571719/149‑1 (MQ=255)
      ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGa                                                                                                                                              <  3:77108/149‑1 (MQ=255)
      ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGa                                                                                                                                              >  4:77108/1‑149 (MQ=255)
              gcgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                                    <  4:417984/135‑1 (MQ=255)
              gcgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                                    >  3:417984/1‑135 (MQ=255)
                                 tACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGc                                                                                                                   >  1:931614/1‑149 (MQ=255)
                                    gCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGaa                                                                                                                >  1:919807/1‑149 (MQ=255)
                                     ccATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         <  1:1215482/139‑1 (MQ=255)
                                     ccATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         >  2:1215482/1‑139 (MQ=255)
                                       aTCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCt                                                                                                             <  3:1435564/149‑1 (MQ=255)
                                          cAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGc                                                                                                                                                 <  2:1256929/110‑1 (MQ=255)
                                          cAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGc                                                                                                                                                 >  1:1256929/1‑110 (MQ=255)
                                           aGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAg                                                                                                         >  2:848792/1‑149 (MQ=255)
                                                 tCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACgg                                                                                                   <  1:848792/149‑1 (MQ=255)
                                                 tCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACgg                                                                                                   <  2:931614/149‑1 (MQ=255)
                                                        gagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATaa                                                                                                                      <  4:192867/123‑1 (MQ=255)
                                                        gagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATaa                                                                                                                      >  3:192867/1‑123 (MQ=255)
                                                          gCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                                   <  2:458567/92‑1 (MQ=255)
                                                          gCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                                   >  1:458567/1‑92 (MQ=255)
                                                           cATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         <  2:596949/117‑1 (MQ=255)
                                                           cATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         >  1:596949/1‑117 (MQ=255)
                                                             ttGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGc                                                                                                                                                 >  4:657596/1‑91 (MQ=255)
                                                             ttGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGc                                                                                                                                                 <  3:657596/91‑1 (MQ=255)
                                                             ttGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAgcgc                                                                                                                                          >  2:1335037/1‑98 (MQ=255)
                                                             ttGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAgcgc                                                                                                                                          <  1:1335037/98‑1 (MQ=255)
                                                             ttGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTATCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGAGACGGGCCACGGCTaaa                                                                                       >  2:236831/1‑149 (MQ=255)
                                                              tGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTaaaa                                                                                      >  1:907841/1‑149 (MQ=255)
                                                                gAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         >  4:1344264/1‑112 (MQ=255)
                                                                gAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                         <  3:1344264/112‑1 (MQ=255)
                                                                   aCAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                                 <  2:907841/149‑1 (MQ=255)
                                                                         tAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACgg                                                                                            >  3:178019/1‑132 (MQ=255)
                                                                         tAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACgg                                                                                            <  4:178019/132‑1 (MQ=255)
                                                                          aGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                          >  3:671856/1‑149 (MQ=255)
                                                                          aGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                          >  4:852877/1‑149 (MQ=255)
                                                                          aGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                          >  2:1348013/1‑149 (MQ=255)
                                                                           gAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATaa                                                                         >  4:1138637/1‑149 (MQ=255)
                                                                                gACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                    <  3:852877/149‑1 (MQ=255)
                                                                                gACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                    <  1:1348013/149‑1 (MQ=255)
                                                                                  cTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                  <  4:671856/149‑1 (MQ=255)
                                                                                  cTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                  <  3:1138637/149‑1 (MQ=255)
                                                                                      cggcggTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                                 >  3:332917/1‑130 (MQ=255)
                                                                                      cggcggTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                                 <  4:332917/130‑1 (MQ=255)
                                                                                              tCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCAc                                                                                              >  2:1229912/1‑109 (MQ=255)
                                                                                              tCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCAc                                                                                              <  1:1229912/109‑1 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATaa                                                                         >  4:1109006/1‑127 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATaa                                                                         <  3:1109006/127‑1 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                  >  3:1125863/1‑134 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                  <  4:1125863/134‑1 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                        >  1:903150/1‑144 (MQ=255)
                                                                                                 atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                        <  2:903150/144‑1 (MQ=255)
                                                                                                   ataAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGAc                                                 >  4:1082322/1‑149 (MQ=255)
                                                                                                    taAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            >  1:1066462/1‑137 (MQ=255)
                                                                                                    taAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            <  2:1066462/137‑1 (MQ=255)
                                                                                                        tACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                    <  4:144431/125‑1 (MQ=255)
                                                                                                        tACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                    >  3:144431/1‑125 (MQ=255)
                                                                                                          cGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggcag                                          <  2:919807/149‑1 (MQ=255)
                                                                                                            gATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACg                                                                                             >  3:1065573/1‑96 (MQ=255)
                                                                                                            gATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACg                                                                                             <  4:1065573/96‑1 (MQ=255)
                                                                                                                  aCTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                        <  1:1338533/127‑1 (MQ=255)
                                                                                                                  aCTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                        >  2:1338533/1‑127 (MQ=255)
                                                                                                                    tGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggc                                            >  2:1225830/1‑137 (MQ=255)
                                                                                                                    tGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggc                                            <  1:1225830/137‑1 (MQ=255)
                                                                                                                      aaCGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGaca                              <  1:236831/149‑1 (MQ=255)
                                                                                                                       aCGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAgg                                                                                                                    >  3:1376384/1‑62 (MQ=255)
                                                                                                                       aCGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAgg                                                                                                                    <  4:1376384/62‑1 (MQ=255)
                                                                                                                            cccAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTccc                        >  1:519166/1‑149 (MQ=255)
                                                                                                                              cAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcct                      >  3:910844/1‑149 (MQ=255)
                                                                                                                                     cTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCa               <  3:1082322/149‑1 (MQ=255)
                                                                                                                                      tGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCAt              >  1:1379149/1‑149 (MQ=255)
                                                                                                                                      gctCGGCATTATTATTGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                          >  4:375528/3‑89 (MQ=255)
                                                                                                                                      gctCGGCATTATTATTGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                          <  3:375528/87‑1 (MQ=255)
                                                                                                                                       gTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATc             <  4:910844/149‑1 (MQ=255)
                                                                                                                                          ggCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            >  4:1382586/1‑99 (MQ=255)
                                                                                                                                          ggCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            <  3:1382586/99‑1 (MQ=255)
                                                                                                                                           gCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            <  4:651679/98‑1 (MQ=255)
                                                                                                                                           gCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtc                                                            >  3:651679/1‑98 (MQ=255)
                                                                                                                                           gCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGaca                              <  1:469548/128‑1 (MQ=255)
                                                                                                                                           gCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGaca                              >  2:469548/1‑128 (MQ=255)
                                                                                                                                            cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcc                   <  1:715137/138‑1 (MQ=255)
                                                                                                                                            cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcc                   >  2:715137/1‑138 (MQ=255)
                                                                                                                                            cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGa          >  3:475064/1‑147 (MQ=255)
                                                                                                                                            cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGa          <  4:475064/147‑1 (MQ=255)
                                                                                                                                             aataataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtctc                                                          >  1:552275/1‑98 (MQ=255)
                                                                                                                                             aataataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtctc                                                          <  2:552275/98‑1 (MQ=255)
                                                                                                                                               taataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTc                                                                             >  4:977156/1‑77 (MQ=255)
                                                                                                                                               taataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTc                                                                             <  3:977156/77‑1 (MQ=255)
                                                                                                                                                  taatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCtt  >  3:1082623/1‑149 (MQ=255)
                                                                                                                                                   aatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTa                                                                                                          >  4:1371705/1‑44 (MQ=255)
                                                                                                                                                   aatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTa                                                                                                          <  3:1371705/44‑1 (MQ=255)
                                                                                                                                                    atGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGa                                <  4:212154/117‑1 (MQ=255)
                                                                                                                                                    atGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGa                                >  3:212154/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTT  >  NC_000913/3399757‑3400051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑