Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,128,199 G→A 5.5% A121T (GCT→ACT)  murJ → putative lipid II flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,128,1990GA 5.5% 233.8 / 2.4 73A121T (GCT→ACT) murJputative lipid II flippase
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base A (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

TAAAAGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAG  >  NC_000913/1128051‑1128347
                                                                                                                                                    |                                                                                                                                                    
taAAAGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAg                                                                                                                                                      >  5:217811/1‑149 (MQ=255)
taAAAGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAg                                                                                                                                                      >  7:140786/1‑149 (MQ=255)
 aaaaGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGc                                                                                                                                                     >  8:215768/1‑149 (MQ=255)
  aaaGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCt                                                                                                                                                    <  5:64189/149‑1 (MQ=255)
    aGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGa                                                                                                                                                  <  8:15583/149‑1 (MQ=255)
    aGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGa                                                                                                                                                  <  7:382901/149‑1 (MQ=255)
     gTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGa                                                                                                                                                  >  4:118546/1‑148 (MQ=255)
            gggTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTg                                                                                                                                          >  3:159700/1‑149 (MQ=255)
             ggTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGc                                                                                                                                         <  2:25368/149‑1 (MQ=255)
             gctGACGACGCCACGCTGGTCCTTCTAGCTAATCTTTCTGGCCGGCTGCCACTTGCGCTCGCGGTTGTGACGGGCGCAGGCATGCTCGCCGCACCGTGCGTGATCATGCTGACCTCGCCAGGCTTCGCTGACACAGCTGACAAATTTGc                                                                                                                                         <  2:134328/147‑1 (MQ=255)
                gAAGCCGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCt                                                                                                                                      >  2:334672/1‑149 (MQ=255)
                   gACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGAc                                                                                                                                   >  6:334626/1‑149 (MQ=255)
                         aCGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGAccagc                                                                                                                               >  3:84299/1‑147 (MQ=255)
                         aCGCGGGGCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGAccagcca                                                                                                                             <  4:325010/149‑1 (MQ=255)
                           gcgGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACcagccagc                                                                                                                           >  4:268335/1‑149 (MQ=255)
                                tCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGAGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGctact                                                                                                                       >  8:267110/1‑148 (MQ=255)
                                      tctcTTATGTTTCTGGCCTGCTGACACTTGCGCGGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGAtt                                                                                                                <  6:20055/149‑1 (MQ=255)
                                        tctTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGCTCATGGTGCCCGCGCCAGCCTTCGCTGACACAGCTGCCACATTTGCCCTGACCAGCCAGCTACTAAAGCTTAc                                                                                                              >  6:226627/1‑149 (MQ=255)
                                          ttATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGTCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAACTGACAAAt                                                                                                                                             <  1:228734/116‑1 (MQ=255)
                                          ttATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAACTGACAAAt                                                                                                                                             >  2:228734/1‑116 (MQ=255)
                                          ttATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAACTGACAAAt                                                                                                                                             >  2:222704/1‑116 (MQ=255)
                                          ttATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAACTGACAAAt                                                                                                                                             <  1:222704/116‑1 (MQ=255)
                                           tATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGATGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCtt                                                                                                           >  4:173470/1‑149 (MQ=255)
                                            aTGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCttt                                                                                                          >  1:374040/1‑149 (MQ=255)
                                            aTGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCt                                                                                                            <  1:334672/147‑1 (MQ=255)
                                             tGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTc                                                                                                         <  6:311181/149‑1 (MQ=255)
                                                 tCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCta                                                                                                     >  6:61465/1‑149 (MQ=255)
                                                  cTGGCCTGCTGACACTTGCGCTGTCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCta                                                                                                     <  3:115746/148‑1 (MQ=255)
                                                   tGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTAc                                                                                                              >  2:177267/1‑138 (MQ=255)
                                                   tGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCtata                                                                                                   >  3:47454/1‑149 (MQ=255)
                                                   tGGCCTGCTGACACTTGCGCTGGCGGTTGAGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTAc                                                                                                              <  1:177267/138‑1 (MQ=255)
                                                          ctgACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTg                                                                                              <  3:47062/147‑1 (MQ=255)
                                                          ctgACACTTGCGCTGGCGGTTGTGACGGGCGCTGGCATGCTCGCCGCACCGTGGGGGATCATGGTGACAGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTg                                                                                              <  4:237453/147‑1 (MQ=255)
                                                            gACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGa                                                                                          >  3:315496/1‑149 (MQ=255)
                                                             acacTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGAt                                                                                         >  4:101826/1‑149 (MQ=255)
                                                             acacTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGAGCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGAt                                                                                         <  2:374040/149‑1 (MQ=255)
                                                                   gcgcTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTccc                                                                                    <  6:377188/148‑1 (MQ=255)
                                                                   gcgcTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTcc                                                                                     <  4:47454/147‑1 (MQ=255)
                                                                    cgcTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCt                                                                                   <  7:63825/148‑1 (MQ=255)
                                                                     gcTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTgg                                                                                 >  5:376214/1‑149 (MQ=255)
                                                                       tGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCAAGGTGACCGAGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCACTATATCGTGCTGATCTCCCTGGCg                                                                               >  1:333830/1‑149 (MQ=255)
                                                                           ggTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGCCCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGc                                                                           >  3:210290/1‑149 (MQ=255)
                                                                             ttGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCtg                                                                         <  5:293923/149‑1 (MQ=255)
                                                                             ttGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCtg                                                                         <  5:32280/149‑1 (MQ=255)
                                                                               gtgACGGTCGCTGGCATGCTCGCCGCACAGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCtgg                                                                        <  3:309153/148‑1 (MQ=255)
                                                                                tgACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCtggtg                                                                      <  5:334626/149‑1 (MQ=255)
                                                                                tgACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGGGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCtggtg                                                                      <  7:39274/149‑1 (MQ=255)
                                                                                  aCGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGCTGACCGCGCCAGGCTACGCGGACACAGCGGACAAATTTGCCCTCACCAGCCCTCTACTAAAGATTACAAAACCCTATATCTTGCTGATCTCCCTGGCTGCGCtcctgg                                                                     >  4:35643/1‑143 (MQ=255)
                                                                                    ggTCGCTGGCATGCTCGCCACACCGTGGGTGATAAAGGTGAACGAGCCTGTCTTCGCTGAAGCAGCTGACAAATTTGCCCTGACCAGCCAGCTACAAAAGATTAACTTTACCTATATCATACTGATCTCCCAGGCTTCGCTGGTGGTAg                                                                  >  4:94172/1‑149 (MQ=255)
                                                                                       ggtcaccATGCTCGCCGCACCGTGGGGGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGa                                                               <  1:176781/143‑1 (MQ=255)
                                                                                        gCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACCAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGa                                                               <  4:84299/148‑1 (MQ=255)
                                                                                            gCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCt                                                           <  3:101826/148‑1 (MQ=255)
                                                                                            gCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCt                                                           <  4:315496/148‑1 (MQ=255)
                                                                                            gCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCt                                                           >  8:215818/1‑148 (MQ=255)
                                                                                                gCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTAACCCTATATCTTGCTGATCTCCCTCGCGTCGCTCGTGGGAGCGATTCTGAATa                                                      >  7:242738/1‑149 (MQ=255)
                                                                                                  tCGCCGCCCCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACg                                                    >  6:307244/1‑149 (MQ=255)
                                                                                                    gccgcACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGt                                                   <  8:140786/148‑1 (MQ=255)
                                                                                                    gccgcACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTg                                                  >  3:356739/1‑149 (MQ=255)
                                                                                                        cACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAAc                                              <  3:268335/149‑1 (MQ=255)
                                                                                                           cGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGc                                           <  7:215768/149‑1 (MQ=255)
                                                                                                             tGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCtt                                         <  4:331156/149‑1 (MQ=255)
                                                                                                                 tGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTtctc                                      >  4:177774/1‑148 (MQ=255)
                                                                                                                   aTCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGAt                                   >  3:30526/1‑149 (MQ=255)
                                                                                                                    tCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGAtt                                  >  3:109937/1‑149 (MQ=255)
                                                                                                                       tGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCg                               >  8:242856/1‑149 (MQ=255)
                                                                                                                            aCCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGt                           <  5:61465/148‑1 (MQ=255)
                                                                                                                             ccGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGttt                         <  4:210290/149‑1 (MQ=255)
                                                                                                                                cgcCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCt                      <  3:173470/149‑1 (MQ=255)
                                                                                                                                   cAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGAtt                                  <  5:360275/134‑1 (MQ=255)
                                                                                                                                   cAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGAtt                                  >  6:360275/1‑134 (MQ=255)
                                                                                                                                   cAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTcc                    >  6:261152/1‑148 (MQ=255)
                                                                                                                                          cGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTCTATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACt              <  5:226627/147‑1 (MQ=255)
                                                                                                                                           gCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGc            <  4:155138/148‑1 (MQ=255)
                                                                                                                                                 acaGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACAt     <  6:217811/149‑1 (MQ=255)
                                                                                                                                                  caGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATc    <  7:267110/149‑1 (MQ=255)
                                                                                                                                                    gCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAg  >  6:149431/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAAAGTAAGCAGGGTGAAGACGCCACGCGGGTCTTTGTCTCTTATGTTTCTGGCCTGCTGACACTTGCGCTGGCGGTTGTGACGGTCGCTGGCATGCTCGCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAG  >  NC_000913/1128051‑1128347

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