Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,399,899 A→T 5.1% I27N (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,8990AT 5.1% 293.2 / 4.3 79I27N (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 = 9.99e-01

TTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGA  >  NC_000913/3399751‑3400043
                                                                                                                                                    |                                                                                                                                                
ttCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaa                                                                                                                                                  <  4:1330952/149‑1 (MQ=255)
     gTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaa                                                                                                                                                  >  4:354415/1‑144 (MQ=255)
     gTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaa                                                                                                                                                  <  3:354415/144‑1 (MQ=255)
      ttCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                            <  2:1279961/149‑1 (MQ=255)
      ttCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                            >  1:1279961/1‑149 (MQ=255)
       tCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                           <  4:930913/149‑1 (MQ=255)
         tAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCATTATTATTGc                                                                                                                                         >  2:571719/1‑149 (MQ=255)
            ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCATTATTATTGCCGa                                                                                                                                      <  1:571719/149‑1 (MQ=255)
            ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGa                                                                                                                                      >  4:77108/1‑149 (MQ=255)
            ttGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGa                                                                                                                                      <  3:77108/149‑1 (MQ=255)
                  aCGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataat                                                                                                                                              >  4:803925/1‑135 (MQ=255)
                  aCGCGAGGCCAGCGACTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaacaat                                                                                                                                              <  3:803925/135‑1 (MQ=255)
                    gcgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                            >  3:417984/1‑135 (MQ=255)
                    gcgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataa                                                                                                                                            <  4:417984/135‑1 (MQ=255)
                           cAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataa                                                                                                                                               <  2:557526/125‑1 (MQ=255)
                           cAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataa                                                                                                                                               >  1:557526/1‑125 (MQ=255)
                           cAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaa                                                                                                                                                  <  3:409325/122‑1 (MQ=255)
                           cAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaa                                                                                                                                                  >  4:409325/1‑122 (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)
                                                             ggagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataat                                                                                                                                              <  1:1327939/92‑1 (MQ=255)
                                                             ggagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataat                                                                                                                                              >  2:1327939/1‑92 (MQ=255)
                                                              gagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATaa                                                                                                              >  3:192867/1‑123 (MQ=255)
                                                              gagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATaa                                                                                                              <  4:192867/123‑1 (MQ=255)
                                                                gCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                           >  1:458567/1‑92 (MQ=255)
                                                                gCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCaataataat                                                                                                                                           <  2:458567/92‑1 (MQ=255)
                                                                 cATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                 >  1:596949/1‑117 (MQ=255)
                                                                 cATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGaa                                                                                                                 <  2:596949/117‑1 (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                                                          <  3:1138637/149‑1 (MQ=255)
                                                                                        cTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                          <  4:671856/149‑1 (MQ=255)
                                                                                            cggcggTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                         <  4:332917/130‑1 (MQ=255)
                                                                                            cggcggTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                         >  3:332917/1‑130 (MQ=255)
                                                                                                    tCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCAc                                                                                      <  1:1229912/109‑1 (MQ=255)
                                                                                                    tCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCAc                                                                                      >  2:1229912/1‑109 (MQ=255)
                                                                                                       atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATaa                                                                 <  3:1109006/127‑1 (MQ=255)
                                                                                                       atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATaa                                                                 >  4:1109006/1‑127 (MQ=255)
                                                                                                       atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                          <  4:1125863/134‑1 (MQ=255)
                                                                                                       atatAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                          >  3:1125863/1‑134 (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                                                                                     <  4:1065573/96‑1 (MQ=255)
                                                                                                                  gATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACg                                                                                     >  3:1065573/1‑96 (MQ=255)
                                                                                                                        aCTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                >  2:1338533/1‑127 (MQ=255)
                                                                                                                        aCTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGt                                                <  1:1338533/127‑1 (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                                                                  <  3:375528/87‑1 (MQ=255)
                                                                                                                                            gctCGGCATTATTATTGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATa                                                                  >  4:375528/3‑89 (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                      >  2:469548/1‑128 (MQ=255)
                                                                                                                                                 gCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGaca                      <  1:469548/128‑1 (MQ=255)
                                                                                                                                                  cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcc           >  2:715137/1‑138 (MQ=255)
                                                                                                                                                  cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcc           <  1:715137/138‑1 (MQ=255)
                                                                                                                                                  cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGa  <  4:475064/147‑1 (MQ=255)
                                                                                                                                                  cAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGa  >  3:475064/1‑147 (MQ=255)
                                                                                                                                                   aataataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtctc                                                  <  2:552275/98‑1 (MQ=255)
                                                                                                                                                   aataataatGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATtctc                                                  >  1:552275/1‑98 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGA  >  NC_000913/3399751‑3400043

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: ‑