Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 968,415 T→C 6.3% L17P (CTG→CCG)  lpxK → lipid A 4'kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,4150TC 6.3% 455.7 / 13.3 142L17P (CTG→CCG) lpxKlipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base T (59/74);  new base C (4/5);  total (63/79)
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.85e-01

Only 100 of 150 total aligned reads displayed.

GTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTC  >  NC_000913/968268‑968563
                                                                                                                                                   |                                                                                                                                                    
gtcgtcGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctg                                                                                                                                                     <  1:1207430/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                                                                                                                                              <  3:1138407/149‑4 (MQ=255)
       aGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaa                                                                                                                                              <  4:210656/149‑4 (MQ=255)
       aGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaa                                                                                                                                              <  4:1133773/149‑4 (MQ=255)
         gATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTcc                                                                                                                                            <  2:929811/149‑1 (MQ=255)
             gtgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGt                                                                                                                                        >  3:746830/1‑149 (MQ=255)
              tgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGtt                                                                                                                                       <  3:455563/149‑1 (MQ=255)
              tgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                             <  4:1449220/143‑5 (MQ=255)
               gtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTTCCACTCTCCTGGTTg                                                                                                                                      >  3:1083138/1‑149 (MQ=255)
                 cATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                                    <  3:875286/149‑1 (MQ=255)
                  aTTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                   >  2:768316/1‑149 (MQ=255)
                  aTTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                   >  1:1435264/1‑149 (MQ=255)
                   ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                                  >  3:660448/1‑149 (MQ=255)
                        gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                             >  2:440393/1‑149 (MQ=255)
                        gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                             <  1:859389/149‑1 (MQ=255)
                          aCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGt                                                                                                                           <  3:538396/149‑1 (MQ=255)
                            gcgGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGa                                                                                                                         <  4:479741/149‑1 (MQ=255)
                              ggTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGt                                                                                                                       <  3:1265087/149‑1 (MQ=255)
                               gTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTg                                                                                                                      <  3:576747/149‑1 (MQ=255)
                                tACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTgg                                                                                                                     <  2:453045/149‑1 (MQ=255)
                                   gCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGgcgc                                                                                                                  >  3:959087/1‑149 (MQ=255)
                                      tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                                                                                               <  4:178577/149‑1 (MQ=255)
                                      tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                                                                                               <  3:896841/149‑1 (MQ=255)
                                      tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                                                                                               >  4:1264173/1‑149 (MQ=255)
                                       aaCGATTTGCTTGAGCACCGCGGCGTTTACGTGCAACTTCACAAAATGCAGTTTGACCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTTTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                              <  1:457259/149‑1 (MQ=255)
                                       aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                   >  3:400730/1‑128 (MQ=255)
                                       aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                              >  4:86705/1‑149 (MQ=255)
                                       aaCGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                              >  1:1061472/1‑149 (MQ=255)
                                          gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                           >  1:185466/1‑149 (MQ=255)
                                          gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                           >  4:837286/1‑149 (MQ=255)
                                          gATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGt                                                                                                           >  1:301762/1‑149 (MQ=255)
                                            tttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                      <  3:1556724/120‑1 (MQ=255)
                                            tttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                      >  4:1556724/1‑120 (MQ=255)
                                               gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                      >  3:1154253/1‑149 (MQ=255)
                                               gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                      <  2:1100036/149‑1 (MQ=255)
                                                cTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGc                                                                                                     <  3:1500416/149‑1 (MQ=255)
                                                 ttGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                                    >  2:659279/1‑149 (MQ=255)
                                                  tGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCta                                                                                                   >  1:1303685/1‑149 (MQ=255)
                                                   gAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                                  >  3:499438/1‑149 (MQ=255)
                                                   gAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                                  <  2:81125/149‑1 (MQ=255)
                                                     gCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACtctc                                                                                                                                             <  3:38874/104‑1 (MQ=255)
                                                     gCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACtctc                                                                                                                                             >  4:38874/1‑104 (MQ=255)
                                                             gcgTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTaa                                                                                           <  2:568831/146‑1 (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)
                                                                    cgcgCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcg                                                                                 <  2:628310/149‑1 (MQ=255)
                                                                     gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                   >  4:227084/1‑98 (MQ=255)
                                                                     gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                   <  3:227084/98‑1 (MQ=255)
                                                                     gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcgc                                                                                >  1:578225/1‑149 (MQ=255)
                                                                       gcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcgcgc                                                                              >  3:311241/1‑149 (MQ=255)
                                                                         aaCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGTGCGATCCGTCTTTGc                                                                                                     <  4:1161464/124‑1 (MQ=255)
                                                                         aaCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGc                                                                                                     >  3:1161464/1‑124 (MQ=255)
                                                                          aCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTg                                                                           >  1:1439987/1‑149 (MQ=255)
                                                                            ttCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGc                                                                         <  2:1061472/149‑1 (MQ=255)
                                                                             tCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTaaaa                                                                                         >  1:764340/1‑132 (MQ=255)
                                                                               acaAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTg                                                                      >  1:364740/1‑149 (MQ=255)
                                                                                caAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGc                                                                     <  1:467146/149‑1 (MQ=255)
                                                                                    aTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAATCTAAAACTGAAGCGCGCCTGGCGTGccccc                                                                 <  4:746830/149‑1 (MQ=255)
                                                                                    aTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGccccc                                                                 <  3:755066/149‑1 (MQ=255)
                                                                                         gTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTAcc                                                            >  3:1327396/1‑149 (MQ=255)
                                                                                           ttGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCgg                                                          <  2:1203027/149‑1 (MQ=255)
                                                                                            tGGCAAATGATCGAAAACATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGTTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAGCTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGt                                                         <  1:1458608/149‑1 (MQ=255)
                                                                                             ggCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGtt                                                        <  1:1264280/149‑1 (MQ=255)
                                                                                             ggCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGtt                                                        <  2:886585/149‑1 (MQ=255)
                                                                                                cAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTc                                                     >  4:1301186/1‑149 (MQ=255)
                                                                                                     aTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggt                                                >  1:1449239/1‑149 (MQ=255)
                                                                                                      tCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                               <  2:1449239/149‑1 (MQ=255)
                                                                                                       cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                              <  3:184336/149‑1 (MQ=255)
                                                                                                       cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                              <  4:584169/149‑1 (MQ=255)
                                                                                                       cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                              <  2:1439987/149‑1 (MQ=255)
                                                                                                        gAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttgg                                             <  4:111885/149‑1 (MQ=255)
                                                                                                         aaaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                            >  4:372796/1‑149 (MQ=255)
                                                                                                          aaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTa                                           <  1:659279/149‑1 (MQ=255)
                                                                                                             aaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAAtc                                        <  3:1228619/149‑1 (MQ=255)
                                                                                                                        gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                                                                                       <  3:1061565/75‑1 (MQ=255)
                                                                                                                        gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                                                                                       >  4:1061565/1‑75 (MQ=255)
                                                                                                                        gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                              <  2:614797/132‑1 (MQ=255)
                                                                                                                        gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggc                             >  2:940540/1‑149 (MQ=255)
                                                                                                                           aaTCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                              >  3:705676/1‑129 (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)
                                                                                                                                ccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                            >  4:914222/1‑126 (MQ=255)
                                                                                                                                 cTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGa                    <  1:925126/149‑1 (MQ=255)
                                                                                                                                  tttGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                                    >  2:938039/1‑68 (MQ=255)
                                                                                                                                  tttGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                                    <  1:938039/68‑1 (MQ=255)
                                                                                                                                     gtgGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGaaaaa                <  1:1568242/149‑1 (MQ=255)
                                                                                                                                      tgGCAGCAATAGCCGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAc               >  4:1392535/1‑149 (MQ=255)
                                                                                                                                        gcggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAccc             <  4:499438/149‑1 (MQ=255)
                                                                                                                                          ggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCg           <  2:364740/149‑1 (MQ=255)
                                                                                                                                             tATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTc        <  3:1466480/149‑1 (MQ=255)
                                                                                                                                                tgctgcCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGtt     <  2:536161/149‑1 (MQ=255)
                                                                                                                                                tgctgcCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGATCGtt     <  4:61560/149‑1 (MQ=255)
                                                                                                                                                   tgcCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTTCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTc  >  2:1086120/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTC  >  NC_000913/968268‑968563

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: