Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,140 T→C L53L (TTG→CTG)  rnlB → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,1400TC98.8% 304.7 / ‑4.2 80L53L (TTG→CTG) rnlBCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (42/37);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CCGGAATTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCTTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAATCAT  >  NC_000913/2766997‑2767286
                                                                                                                                               |                                                                                                                                                  
ccGGAATTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGGGACCACTGGTGAAATCCTTTTTGACCTCCTGTgcg                                                                                                                                              <  2:715609/150‑1 (MQ=255)
 cGGAATTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCt                                                                                                                                            >  2:849940/1‑151 (MQ=255)
 cGGAATTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCt                                                                                                                                            >  2:314189/1‑151 (MQ=255)
    aaTTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAAc                                                                                                                                         <  1:1218877/151‑1 (MQ=255)
        aaCGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTc                                                                                                                                     >  2:409852/1‑151 (MQ=255)
        aaCGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTc                                                                                                                                     <  2:284591/151‑1 (MQ=255)
           gTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCag                                                                                                                                  <  2:646089/151‑1 (MQ=255)
             ctcAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCag                                                                                                                                  >  2:553762/1‑149 (MQ=255)
              tcAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGt                                                                                                                               >  1:694041/1‑151 (MQ=255)
               cAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGt                                                                                                                               <  2:251968/150‑1 (MQ=255)
               cAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTg                                                                                                                              >  1:810580/1‑151 (MQ=255)
                 ggggCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGa                                                                                                                            <  1:13902/151‑1 (MQ=255)
                    gCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAAcc                                                                                                                         <  1:981153/151‑1 (MQ=255)
                    gCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAAcc                                                                                                                         >  2:1168672/1‑151 (MQ=255)
                      tctAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGc                                                                                                                       >  1:625216/1‑151 (MQ=255)
                      tctAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGc                                                                                                                       >  1:1016570/1‑151 (MQ=255)
                       ctAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCt                                                                                                                      >  2:589489/1‑151 (MQ=255)
                         aaaaGCTTGCGTTATGGCGACCGGCTTCGACAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTGTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTc                                                                                                                    <  2:385894/151‑1 (MQ=255)
                         aaaaGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGc                                                                                                                       <  2:148948/148‑1 (MQ=255)
                         aaaaGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGc                                                                                                                       <  1:606156/148‑1 (MQ=255)
                         aaaaGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTc                                                                                                                    <  1:998932/151‑1 (MQ=255)
                         aaaaGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTc                                                                                                                    >  1:811664/1‑151 (MQ=255)
                            aGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGtt                                                                                                                 >  2:774450/1‑151 (MQ=255)
                             gCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGtt                                                                                                                 >  2:359708/1‑150 (MQ=255)
                              cTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTAc                                                                                                               >  2:919057/1‑151 (MQ=255)
                                gTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACg                                                                                                              <  1:1017908/150‑1 (MQ=255)
                                gTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACg                                                                                                              >  2:1045915/1‑150 (MQ=255)
                                gTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACg                                                                                                              <  2:82384/150‑1 (MQ=255)
                                     tATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCTTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGga                                                                                                         <  1:1211386/150‑1 (MQ=255)
                                      aTGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCACTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGgaa                                                                                                        >  1:931530/1‑150 (MQ=255)
                                      aTGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTCTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGgaa                                                                                                        >  1:1028879/1‑150 (MQ=255)
                                           gACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGgaaatgaa                                                                                                   >  1:652247/1‑150 (MQ=255)
                                            aCCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGgaaatgaaat                                                                                                 <  2:375204/151‑1 (MQ=255)
                                             ccGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGga                                                                                                         <  1:503809/142‑1 (MQ=255)
                                             ccGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGga                                                                                                         >  2:503751/1‑142 (MQ=255)
                                                  ttCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGcc                                                                                           >  2:1083927/1‑151 (MQ=255)
                                                     gagaACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGCAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCgg                                                                                         >  1:221494/1‑150 (MQ=255)
                                                         aCCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCAt                                                                                    <  2:1185484/151‑1 (MQ=255)
                                                              ttGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCt                                                                                 <  1:200196/149‑1 (MQ=255)
                                                               tGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTgg                                                                               >  2:865248/1‑150 (MQ=255)
                                                                gTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGa                                                                              <  2:4470/150‑1 (MQ=255)
                                                                gTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGAt                                                                             <  1:1384358/151‑1 (MQ=255)
                                                                gTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGAt                                                                             <  2:696013/151‑1 (MQ=255)
                                                                  ccTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATAc                                                                           <  1:937805/151‑1 (MQ=255)
                                                                   cTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACg                                                                          >  2:519854/1‑151 (MQ=255)
                                                                      aGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGcc                                                                       <  1:323517/151‑1 (MQ=255)
                                                                        ttAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCa                                                                      >  1:659405/1‑150 (MQ=255)
                                                                            tGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAatt                                                                 >  2:1147697/1‑151 (MQ=255)
                                                                             gAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAatt                                                                 >  2:356157/1‑150 (MQ=255)
                                                                                aTTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATtg                                                             >  2:255341/1‑151 (MQ=255)
                                                                                   gagaCGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATtgatg                                                          <  2:1329888/151‑1 (MQ=255)
                                                                                    agaCGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATtgatg                                                          <  1:1046011/150‑1 (MQ=255)
                                                                                    agaCGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATtgatg                                                          >  1:947061/1‑150 (MQ=255)
                                                                                     gaCGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAg                                                        <  2:124617/151‑1 (MQ=255)
                                                                                      aCGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGc                                                       >  1:422082/1‑151 (MQ=255)
                                                                                         aaaCTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAgg                                                    >  2:130678/1‑151 (MQ=255)
                                                                                         aaaCTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAgg                                                    <  1:845273/151‑1 (MQ=255)
                                                                                          aaCTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGa                                                   <  2:634907/151‑1 (MQ=255)
                                                                                           aCTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGAc                                                  >  1:1384592/1‑151 (MQ=255)
                                                                                            cTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACg                                                 >  1:908047/1‑151 (MQ=255)
                                                                                                cTGCCCTACTGGGTAGTGAGACCACTGGAGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGCTCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTc                                              >  1:1151781/1‑150 (MQ=255)
                                                                                                  gCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTcc                                             <  1:8833/149‑1 (MQ=255)
                                                                                                  gCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCg                                            <  1:1082873/150‑1 (MQ=255)
                                                                                                  gCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGa                                           <  1:1121339/151‑1 (MQ=255)
                                                                                                   cccTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGac                                          >  2:519885/1‑151 (MQ=255)
                                                                                                       aCTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGAcaca                                       <  1:163905/150‑1 (MQ=255)
                                                                                                         tGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACatat                                    >  2:205970/1‑151 (MQ=255)
                                                                                                              aGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGa                               <  2:48937/151‑1 (MQ=255)
                                                                                                              aGTGAGACCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGa                               <  1:850026/151‑1 (MQ=255)
                                                                                                                  agaCCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCa                            >  1:1003248/1‑150 (MQ=255)
                                                                                                                    aCCACTGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAAc                         >  1:1120380/1‑151 (MQ=255)
                                                                                                                         tGGTGAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACt                    >  2:1197333/1‑151 (MQ=255)
                                                                                                                             gAAATCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACtttt                 >  2:1270568/1‑150 (MQ=255)
                                                                                                                              aaaTCCTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTtttac               >  1:227374/1‑151 (MQ=255)
                                                                                                                               aaTCCTTTTTGAGCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGACGGATCATGCTGGAGACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTtttact              >  2:104231/1‑151 (MQ=255)
                                                                                                                                  ccTTTTTGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTtttacttt            <  2:1037218/150‑1 (MQ=255)
                                                                                                                                      tttGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGt        >  2:374379/1‑150 (MQ=255)
                                                                                                                                      tttGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGt        <  1:1221691/150‑1 (MQ=255)
                                                                                                                                       ttGACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTa       <  1:1168769/150‑1 (MQ=255)
                                                                                                                                         gACCTCCTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAat     <  2:1028781/150‑1 (MQ=255)
                                                                                                                                           cctcctGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAatcat  <  2:693963/151‑1 (MQ=255)
                                                                                                                                           cctcctGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAatcat  >  1:838804/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCGGAATTAACGTCTCAGGGGCTCTAAAAGCTGTCGTTATGGCGACCGGCTTCGAGAACCCTTTGTCCTCAGTTAATGAAATTGAGACGAAACTTTCTGCCCTACTGGGTAGTGAGACCACTGGTGAAATCCTTTTTGACCTCTTGTGCGCTAACGGTCCAGAGTGGAACCGCTTCGTTACGCTGGAAATGAAATATGGCCGGATCATGCTGGATACGGCCAAAATTATTGATGAGCAGGACGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAATCAT  >  NC_000913/2766997‑2767286

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