Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 968,409 T→A 6.1% L15Q (CTA→CAA)  lpxK → lipid A 4'kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,4090TA 6.1% 517.6 / 17.5 148L15Q (CTA→CAA) lpxKlipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base T (62/77);  new base A (4/5);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 150 total aligned reads displayed.

ATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTC  >  NC_000913/968262‑968557
                                                                                                                                                   |                                                                                                                                                    
aTCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTa                                                                                                                                                     >  1:1271355/1‑149 (MQ=255)
aTCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTa                                                                                                                                                     >  1:628310/1‑149 (MQ=255)
    tGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCGTTGTGGCGGCTATtgc                                                                                                                                                 >  1:708044/1‑149 (MQ=255)
      gtcgtcGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctg                                                                                                                                               <  1:1207430/149‑1 (MQ=255)
      gtcgtcGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctg                                                                                                                                               <  3:1517779/149‑1 (MQ=255)
       tcgtcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctgc                                                                                                                                              <  2:1585920/149‑1 (MQ=255)
       tcgtcgAGGATGGGGTCATTGTGGACCGCGGTACGCATAACGTTTTGCTTGAGCACCGCGGCGTTTACGCGAAATTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctgc                                                                                                                                              <  4:935356/149‑1 (MQ=255)
          tcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCAc                                                                                                                                           <  1:917863/149‑1 (MQ=255)
             aGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaa                                                                                                                                        <  4:210656/149‑4 (MQ=255)
             aGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaa                                                                                                                                        <  3:1138407/149‑4 (MQ=255)
              ggATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGATTCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                       <  3:501723/149‑5 (MQ=255)
               gATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTcc                                                                                                                                      <  2:929811/149‑1 (MQ=255)
                aTGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCt                                                                                                                                     <  3:1089462/149‑1 (MQ=255)
                   gtgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGt                                                                                                                                  >  3:746830/1‑149 (MQ=255)
                    tgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGtt                                                                                                                                 <  3:455563/149‑1 (MQ=255)
                    tgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                       >  3:1449220/1‑139 (MQ=255)
                       cATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                              <  4:1432291/149‑1 (MQ=255)
                       cATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                              <  3:875286/149‑1 (MQ=255)
                        aTTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                             >  1:1435264/1‑149 (MQ=255)
                         ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                            >  1:436595/1‑149 (MQ=255)
                         ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                            >  3:660448/1‑149 (MQ=255)
                         ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                            >  3:584169/1‑149 (MQ=255)
                            tgGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGcc                                                                                                                         >  3:347078/1‑149 (MQ=255)
                              gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                       >  2:440393/1‑149 (MQ=255)
                                aCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGt                                                                                                                     <  3:538396/149‑1 (MQ=255)
                                 cgcgGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTg                                                                                                                    >  2:932886/1‑149 (MQ=255)
                                  gcgGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGa                                                                                                                   <  4:479741/149‑1 (MQ=255)
                                    ggTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGt                                                                                                                 <  3:1265087/149‑1 (MQ=255)
                                     gTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTg                                                                                                                >  2:1065859/1‑149 (MQ=255)
                                      tACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTgg                                                                                                               <  2:453045/149‑1 (MQ=255)
                                         gCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGgcgc                                                                                                            >  3:959087/1‑149 (MQ=255)
                                            tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                                                                                         >  4:1264173/1‑149 (MQ=255)
                                             aaCGATTTGCTTGAGCACCGCGGCGTTTACGTGCAACTTCACAAAATGCAGTTTGACCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTTTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                        <  1:457259/149‑1 (MQ=255)
                                             aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                             <  4:400730/128‑1 (MQ=255)
                                             aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                        >  1:1061472/1‑149 (MQ=255)
                                             aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                        >  4:86705/1‑149 (MQ=255)
                                              aCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATcc                                                                                                       <  3:726659/149‑1 (MQ=255)
                                                gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                     >  1:301762/1‑149 (MQ=255)
                                                gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                     >  1:697437/1‑149 (MQ=255)
                                                gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                     >  4:837286/1‑149 (MQ=255)
                                                  tttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                >  4:1556724/1‑120 (MQ=255)
                                                  tttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                <  3:1556724/120‑1 (MQ=255)
                                                     gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                <  4:347078/149‑1 (MQ=255)
                                                     gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                <  2:621578/149‑1 (MQ=255)
                                                     gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                <  2:1100036/149‑1 (MQ=255)
                                                       ttGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                              >  2:659279/1‑149 (MQ=255)
                                                       ttGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                              >  1:895320/1‑149 (MQ=255)
                                                         gAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                            <  2:81125/149‑1 (MQ=255)
                                                           gCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACtctc                                                                                                                                       <  3:38874/104‑1 (MQ=255)
                                                             aCCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAAc                                                                                        >  4:1228619/1‑149 (MQ=255)
                                                             aCCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAAc                                                                                        >  1:1483683/1‑149 (MQ=255)
                                                                   gcgTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTaa                                                                                     >  1:568831/1‑146 (MQ=255)
                                                                      tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                             <  2:1141997/103‑1 (MQ=255)
                                                                      tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                             >  1:1141997/1‑103 (MQ=255)
                                                                      tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGa                                                                               >  1:67570/1‑149 (MQ=255)
                                                                           gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                             <  3:227084/98‑1 (MQ=255)
                                                                           gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcgc                                                                          >  1:536161/1‑149 (MQ=255)
                                                                             gcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcgcgc                                                                        >  3:11585/1‑149 (MQ=255)
                                                                               aaCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGc                                                                                               >  3:1161464/1‑124 (MQ=255)
                                                                                  ttCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGc                                                                   <  3:1264173/149‑1 (MQ=255)
                                                                                   tCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTaaaa                                                                                   >  1:764340/1‑132 (MQ=255)
                                                                                   tCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTaaaa                                                                                   <  2:764340/132‑1 (MQ=255)
                                                                                    cacaAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGt                                                                 <  2:866304/149‑1 (MQ=255)
                                                                                      caAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGc                                                               <  3:837286/149‑1 (MQ=255)
                                                                                          aTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAATCTAAAACTGAAGCGCGCCTGGCGTGccccc                                                           <  4:746830/149‑1 (MQ=255)
                                                                                                 ttGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCgg                                                    <  2:1203027/149‑1 (MQ=255)
                                                                                                 ttGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCgg                                                    <  4:1083138/149‑1 (MQ=255)
                                                                                                   ggCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGtt                                                  <  4:1134627/149‑1 (MQ=255)
                                                                                                   ggCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGtt                                                  <  1:1264280/149‑1 (MQ=255)
                                                                                                      cggacaTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTc                                               <  1:259121/144‑1 (MQ=255)
                                                                                                      cAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTc                                               >  4:1301186/1‑149 (MQ=255)
                                                                                                         tGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtg                                            <  3:45240/149‑1 (MQ=255)
                                                                                                           aTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggt                                          >  1:1449239/1‑149 (MQ=255)
                                                                                                            tCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                         <  2:1449239/149‑1 (MQ=255)
                                                                                                             cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                        <  4:584169/149‑1 (MQ=255)
                                                                                                             cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                        <  2:1439987/149‑1 (MQ=255)
                                                                                                              gAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGTCCTCGTACCGGTTGTCGtggttgg                                       >  4:1480073/1‑149 (MQ=255)
                                                                                                              gAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttgg                                       <  4:111885/149‑1 (MQ=255)
                                                                                                               aaaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                      >  4:372796/1‑149 (MQ=255)
                                                                                                                   aaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAAtc                                  <  1:932886/149‑1 (MQ=255)
                                                                                                                              gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                                                                                 <  3:1061565/75‑1 (MQ=255)
                                                                                                                              gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                        >  1:614797/1‑132 (MQ=255)
                                                                                                                              gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                        <  2:614797/132‑1 (MQ=255)
                                                                                                                              gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggc                       >  2:940540/1‑149 (MQ=255)
                                                                                                                                 aaTCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                        <  4:705676/129‑1 (MQ=255)
                                                                                                                                 aaTCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggcggc                    <  2:301762/149‑1 (MQ=255)
                                                                                                                                   tCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCaa                  <  4:1157460/149‑1 (MQ=255)
                                                                                                                                   tCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCaa                  >  1:1006914/1‑149 (MQ=255)
                                                                                                                                     aggagagtgGCAGCAATAGCCGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACg                >  4:515120/7‑149 (MQ=255)
                                                                                                                                      ccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                      <  3:914222/126‑1 (MQ=255)
                                                                                                                                      ccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                      >  4:914222/1‑126 (MQ=255)
                                                                                                                                      ccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                         >  4:190902/1‑123 (MQ=255)
                                                                                                                                      ccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                         <  3:190902/123‑1 (MQ=255)
                                                                                                                                       gagagtgGCAGCAATAGCCGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGa              >  2:540577/5‑149 (MQ=255)
                                                                                                                                        tttGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                              <  1:938039/68‑1 (MQ=255)
                                                                                                                                           gtgGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGaaaaa          <  1:1568242/149‑1 (MQ=255)
                                                                                                                                           gtgGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGaaaaa          <  1:1442477/149‑1 (MQ=255)
                                                                                                                                                ggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCg     <  2:364740/149‑1 (MQ=255)
                                                                                                                                                   tATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTc  <  1:1318196/149‑1 (MQ=255)
                                                                                                                                                   tATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTc  <  3:1466480/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTC  >  NC_000913/968262‑968557

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