Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,303,658 T→A 5.5% Y254N (TAT→AAT)  oppB → oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,303,6580TA 5.5% 247.7 / 4.7 92Y254N (TAT→AAT) oppBoligopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (45/42);  new base A (2/3);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

TTATTCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGACCC  >  NC_000913/1303512‑1303803
                                                                                                                                                  |                                                                                                                                                 
ttATTCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATg                                                                                                                                                >  1:775786/1‑150 (MQ=255)
ttATTCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGGCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATg                                                                                                                                                <  2:917530/150‑1 (MQ=255)
tcaTTCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATg                                                                                                                                                >  1:703421/3‑150 (MQ=255)
   ttCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGtt                                                                                                                                             <  1:263002/150‑1 (MQ=255)
   ttCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGtt                                                                                                                                             <  1:99279/150‑1 (MQ=255)
     cGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTTTGGTTTg                                                                                                                                           >  1:1257529/1‑150 (MQ=255)
     cGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTg                                                                                                                                           >  1:1257524/1‑150 (MQ=255)
      gTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATAAATGGtttcg                                                                                                                                          <  1:650890/150‑3 (MQ=255)
       tACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTTCTCTCTTATATGGGCCCCGCTTTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATAAATGGtttcga                                                                                                                                         <  1:1392047/150‑4 (MQ=255)
        attGCCCGGGCGATAGGGTTACCTATGATGCGTATCATTTTTCGTCACGCATTAAAACCTGCTCTTTTTCCCGTGCGCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATAAATGGtttcgat                                                                                                                                        <  1:25348/148‑5 (MQ=255)
         cTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCgg                                                                                                                                       <  1:568175/150‑1 (MQ=255)
         cTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTGGCggg                                                                                                                                       >  2:1530497/1‑150 (MQ=255)
           gCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCCCGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCgggg                                                                                                                                     >  1:784197/1‑150 (MQ=255)
           gCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTTTCGAAACCATTTATGGTTTGCCgggg                                                                                                                                     >  1:300200/1‑150 (MQ=255)
           gCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCgggg                                                                                                                                     >  1:634725/1‑150 (MQ=255)
           gCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCgggg                                                                                                                                     >  1:536726/1‑150 (MQ=255)
             ccGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGAt                                                                                                                                   >  1:1645440/1‑150 (MQ=255)
             ccGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTAGGGTTTGCCGGGGAt                                                                                                                                   >  2:130392/1‑150 (MQ=255)
               gggCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTg                                                                                                                                 <  1:1448748/150‑1 (MQ=255)
                 gCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGa                                                                                                                               >  2:1430238/1‑150 (MQ=255)
                  cGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGAc                                                                                                                              <  1:292081/150‑1 (MQ=255)
                   gAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACa                                                                                                                             >  1:1355946/1‑150 (MQ=255)
                     aaGGGTTACCTATGCGTCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAAc                                                                                                                           <  2:1463604/150‑1 (MQ=255)
                       gggTTACCTATTCGGCGGATCATTTTACGTCACGCTTTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTg                                                                                                                         <  1:1050294/150‑1 (MQ=255)
                       gggTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTg                                                                                                                         >  1:604519/1‑150 (MQ=255)
                                 tGCGGCGGCTCATTTGACGTCACGCATTAACTCCGGCTCTGGTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATg                                                                                                               <  2:1205863/150‑1 (MQ=255)
                                    ggctgATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGgcg                                                                                                            <  2:1279613/150‑1 (MQ=255)
                                       ggATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCAt                                                                                                         >  1:429083/1‑150 (MQ=255)
                                        gATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAATCCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTTATGGCGCAtt                                                                                                        >  1:1025543/1‑150 (MQ=255)
                                         aTCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTg                                                                                                       >  2:402006/1‑150 (MQ=255)
                                         aTCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTg                                                                                                       <  1:538711/150‑1 (MQ=255)
                                            aTTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCAATGAAc                                                                                                    >  2:392306/1‑150 (MQ=255)
                                            aTTTTACGTCAAGCATTAAAACGTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGAATTTGTGTGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTTATGGCGCATTGAAc                                                                                                    <  2:719413/150‑1 (MQ=255)
                                                tACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGGTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTg                                                                                                <  2:698894/150‑1 (MQ=255)
                                                  cGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGAc                                                                                              <  1:1350168/150‑1 (MQ=255)
                                                  cGTCAAGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGAc                                                                                              <  1:1356923/150‑1 (MQ=255)
                                                   gTCACGCATTAAAACCTGCTCTGTTACCCTTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACt                                                                                             <  2:1404419/150‑1 (MQ=255)
                                                   gTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACAGTTCGTTAATGGCGCATTGAACCGTGACt                                                                                             >  1:522837/1‑150 (MQ=255)
                                                   gTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGCGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACt                                                                                             >  2:1426531/1‑150 (MQ=255)
                                                     cACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTAt                                                                                           >  2:165126/1‑150 (MQ=255)
                                                     cACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAAACGTGACTAt                                                                                           >  1:773419/1‑150 (MQ=255)
                                                        gCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTcc                                                                                        <  1:1430238/150‑1 (MQ=255)
                                                                  cTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCTCATTGAACCGTGACTATTCCTTAGTGTTaa                                                                              >  1:1070248/1‑150 (MQ=255)
                                                                       ctGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTg                                                                         >  2:749191/1‑150 (MQ=255)
                                                                         gtaaCCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATATATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGAc                                                                       <  1:1265346/147‑1 (MQ=255)
                                                                            aCCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGAGTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCAt                                                                    <  2:300200/150‑1 (MQ=255)
                                                                               cGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCtt                                                                                      <  1:1102146/129‑1 (MQ=255)
                                                                               cGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCtt                                                                                      >  2:1102146/1‑129 (MQ=255)
                                                                                gTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCtg                                                                <  2:237929/150‑1 (MQ=255)
                                                                                   ctctctTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCtggtt                                                             >  2:1207470/1‑150 (MQ=255)
                                                                                    tctctTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCtggttg                                                            <  2:1645440/150‑1 (MQ=255)
                                                                                       ctTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTg                                                         <  2:478496/150‑1 (MQ=255)
                                                                                        ttATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGc                                                        >  1:776586/1‑150 (MQ=255)
                                                                                         tataTGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCt                                                       <  2:372553/150‑1 (MQ=255)
                                                                                         tataTGGGCCCCGCATTTGTCGGCATTATTACCGGTGCTATGGTTATAGAAACAATTTATGGTTTGCCGGCGATTGGACAAATCTTCGTTAATTGCGCATTGCACCGTGACTATTCATAAGAGTTAAGCCTGACCAACCTGGATGGTCCt                                                       >  1:853092/1‑150 (MQ=255)
                                                                                          atatGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCtt                                                      <  2:175805/150‑1 (MQ=255)
                                                                                           tatGTGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCttt                                                     <  2:339938/150‑1 (MQ=255)
                                                                                             tGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTaa                                                   <  2:634725/150‑1 (MQ=255)
                                                                                              gggCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAAc                                                  <  1:1571860/150‑1 (MQ=255)
                                                                                                gcaacGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCAGTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCAGTGAACCGTGTCTAGTCCTTAGTGTTAAGCCTGACAATCCTGGTTGGTGCTTTAACCa                                                >  2:456935/5‑150 (MQ=255)
                                                                                                 ccccGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCAt                                               <  1:706180/150‑1 (MQ=255)
                                                                                                   ccGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGTTGCTTTAACCAttt                                             >  1:878520/1‑150 (MQ=255)
                                                                                                     cactttGGCGGGATTATTACGCGGGCTCTGGTTATCGACACCATTTCTGGTTTGCCGGGGATTGGACAACTGTTCGTTAATTGCGCATTGACCCGGGACTATGACTTAGTGGTAAGCCTGCCCATCCTGGTTGGTGCTTTAACCATTTTg                                           <  2:301692/147‑1 (MQ=255)
                                                                                                      cttttGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGt                                          <  1:1198524/148‑1 (MQ=255)
                                                                                                      cATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGt                                          <  2:1355946/150‑1 (MQ=255)
                                                                                                      cATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGt                                          >  2:771738/1‑150 (MQ=255)
                                                                                                        tttGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGttt                                        <  2:938145/150‑1 (MQ=255)
                                                                                                        tttGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGttt                                        <  1:1428147/150‑1 (MQ=255)
                                                                                                         ttGTCGGCATTATTTCCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGAGTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTa                                       <  2:201022/150‑1 (MQ=255)
                                                                                                           gTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGACCCGTGACTATTCCTTAGTGTTAAGCCTGACCATCATGGTTGGTGCTTTAACCATTTTGTTTatt                                     >  1:1602263/1‑148 (MQ=255)
                                                                                                              ggCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGcc                                  >  2:691113/1‑150 (MQ=255)
                                                                                                               gCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCa                                 >  1:1680338/1‑150 (MQ=255)
                                                                                                                cATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCAt                                >  1:1476159/1‑150 (MQ=255)
                                                                                                                   tattaCCTGTTCTAGGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGACCCGTGACTATTCCTTAGTGTTAAGCCGGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGt                             <  2:1013305/150‑1 (MQ=255)
                                                                                                                    attaCCGGTTCTATGGTTATCGACACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAAGGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTc                            >  1:1530400/1‑150 (MQ=255)
                                                                                                                    attaCCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTc                            >  2:813175/1‑150 (MQ=255)
                                                                                                                    attaCCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTc                            >  2:519800/1‑150 (MQ=255)
                                                                                                                        ccGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGAtg                        <  2:1509119/150‑1 (MQ=255)
                                                                                                                        ccGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCAATGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGAtg                        <  1:392306/150‑1 (MQ=255)
                                                                                                                              cTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCtat                  <  1:794956/150‑1 (MQ=255)
                                                                                                                                  ggTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGc              >  1:313878/1‑150 (MQ=255)
                                                                                                                                   gTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCg             <  2:1484527/150‑1 (MQ=255)
                                                                                                                                   gTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACAATTTTGTTTAATGCCATTGTCGATGTGCTATATGCg             >  1:859746/1‑150 (MQ=255)
                                                                                                                                     tATCGAAACCATTTATGGTTTGCCGGGGTTTGGACAACTGTGCGTTAATGGCGCCTTGAACCGTGACTATTCCTTAGTGTAAAGCCTGACCATCCTGGGTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGt           <  2:1054849/150‑1 (MQ=255)
                                                                                                                                      cggcaaaCCATAAATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCTTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCAATTGTCGAAGTGCTATATGCTGtt          >  2:1130461/5‑150 (MQ=255)
                                                                                                                                      cggcaaaCCATAAATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGtt          >  1:1270785/5‑150 (MQ=255)
                                                                                                                                      aTCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGtt          >  1:688325/1‑150 (MQ=255)
                                                                                                                                         gAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATc       >  2:56977/1‑150 (MQ=255)
                                                                                                                                         gAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATc       >  2:984010/1‑150 (MQ=255)
                                                                                                                                           aacattttATGGTTTGCCGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGAc    <  1:1663874/145‑1 (MQ=255)
                                                                                                                                            aCCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGAc    <  2:1070248/150‑1 (MQ=255)
                                                                                                                                             ccATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGAcc   >  1:1455844/1‑150 (MQ=255)
                                                                                                                                              cATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGAccc  >  2:1161798/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTATTCGTACTGCCCGGGCGAAAGGGTTACCTATGCGGCGGATCATTTTACGTCACGCATTAAAACCTGCTCTGTTACCCGTGCTCTCTTATATGGGCCCCGCATTTGTCGGCATTATTACCGGTTCTATGGTTATCGAAACCATTTATGGTTTGCCGGGGATTGGACAACTGTTCGTTAATGGCGCATTGAACCGTGACTATTCCTTAGTGTTAAGCCTGACCATCCTGGTTGGTGCTTTAACCATTTTGTTTAATGCCATTGTCGATGTGCTATATGCGGTTATCGACCC  >  NC_000913/1303512‑1303803

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