Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 858,425 G→T 8.1% L213M (CTG→ATG) ‡ ybiU ← DUF1479 domain‑containing protein YbiU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913858,4250GT 8.1% 362.0 / 14.4 110L213M (CTG→ATG) ‡ybiUDUF1479 domain‑containing protein YbiU
Reads supporting (aligned to +/- strand):  ref base G (43/57);  new base T (5/4);  total (49/61)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 123 total aligned reads displayed.

GCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAT  >  NC_000913/858276‑858574
                                                                                                                                                     |                                                                                                                                                     
gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                       <  2:1532916/150‑1 (MQ=255)
gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                       >  1:543535/1‑150 (MQ=255)
gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                       <  1:1326593/150‑1 (MQ=255)
 cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCTTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCGATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAtgcagc                                                                                                                                                      <  2:1486949/150‑1 (MQ=255)
 cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                      >  2:1244252/1‑150 (MQ=255)
  ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGcagcc                                                                                                                                                     <  1:505101/150‑1 (MQ=255)
    tGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGcagccag                                                                                                                                                   <  1:226171/150‑1 (MQ=255)
    tGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGcagccag                                                                                                                                                   <  2:292396/150‑1 (MQ=255)
      gAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCg                                                                                                                                                 <  2:680408/150‑1 (MQ=255)
        aTGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGta                                                                                                                                               <  1:876774/150‑2 (MQ=255)
        aTGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGta                                                                                                                                               <  1:1275527/150‑2 (MQ=255)
         tGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTc                                                                                                                                              >  2:1245223/1‑150 (MQ=255)
          gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                             >  2:1162700/1‑150 (MQ=255)
          gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                             <  2:526068/150‑1 (MQ=255)
           tCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCa                                                                                                                                            >  1:30268/1‑150 (MQ=255)
            ccGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                           <  2:1315688/150‑1 (MQ=255)
             cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatc                                                                                                                                             <  1:1368960/147‑4 (MQ=255)
             cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGt                                                                                                                                          >  1:417436/1‑150 (MQ=255)
             cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGACATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatc                                                                                                                                             >  2:1368960/1‑144 (MQ=255)
              gAAACACGGAACATTTGGTGTTGTTGTCCACCGTGTACTCTTCAACTTCGGTACTTTGTGCCGCATGCCAGGGATCATATTGCTCCAGTTTGCGATTAAAGACGTTGGCGAAATCGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                         <  2:71588/150‑8 (MQ=255)
                  cacGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcccc                                                                                                                                     <  1:1187944/150‑1 (MQ=255)
                   acGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCg                                                                                                                                    <  2:656276/150‑1 (MQ=255)
                   acGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCg                                                                                                                                    >  1:1563729/1‑150 (MQ=255)
                     ggAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGa                                                                                                                                  >  1:390486/1‑150 (MQ=255)
                     ggAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGa                                                                                                                                  >  1:119523/1‑150 (MQ=255)
                          aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                             >  2:968493/1‑150 (MQ=255)
                          aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                             <  1:1026204/150‑1 (MQ=255)
                           tttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGt                                                                                                                            <  1:1435597/150‑1 (MQ=255)
                            ttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTa                                                                                                                           <  2:1512178/150‑1 (MQ=255)
                              ggtgggGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATg                                                                                                                         <  1:734903/150‑1 (MQ=255)
                                tggtgTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATgcg                                                                                                                       <  2:562090/150‑1 (MQ=255)
                                    gttgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTcc                                                                                                                   <  1:273116/150‑1 (MQ=255)
                                    gttgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTcc                                                                                                                   >  2:1572291/1‑150 (MQ=255)
                                     ttgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCa                                                                                                                  <  2:1563729/150‑1 (MQ=255)
                                       gtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAg                                                                                                                <  1:1255759/150‑1 (MQ=255)
                                       gtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAg                                                                                                                >  1:533734/1‑150 (MQ=255)
                                       ctCCACCCTGTAAACTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTCCCATTAAAGACGTTGGCGAAAACGCCCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGACGGAATGCGCTCCaac                                                                                                                >  2:1578167/2‑149 (MQ=255)
                                         ccACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCTCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTc                                                                                                                    <  2:36889/144‑1 (MQ=255)
                                         ccACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCACCGGAGTCGGTATGCGCTc                                                                                                                    >  1:36889/1‑144 (MQ=255)
                                          cACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                             <  1:1570405/150‑1 (MQ=255)
                                           aCCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCt                                                                                                            <  2:931095/150‑1 (MQ=255)
                                              gtgtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTg                                                                                                         <  2:1199518/150‑1 (MQ=255)
                                                   ctctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGt                                                                                                    <  1:1245223/150‑1 (MQ=255)
                                                    tctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTc                                                                                                   >  1:1651713/1‑150 (MQ=255)
                                                     ctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCg                                                                                                  <  1:1004807/150‑1 (MQ=255)
                                                      ttCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGt                                                                                                 >  2:1361302/1‑150 (MQ=255)
                                                        cAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGACAACGCGCTGATACGCTGGCAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGGATGCGCGCCACGACCTTTGGAggtggtgc                                                                                               >  2:874354/1‑148 (MQ=255)
                                                        cAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTc                                                                                               >  2:1122810/1‑150 (MQ=255)
                                                           cTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCgg                                                                                            >  2:213819/1‑150 (MQ=255)
                                                             tCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGc                                                                                          <  1:1543070/150‑1 (MQ=255)
                                                                        tgCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcggcggcg                                                                               <  1:379267/150‑1 (MQ=255)
                                                                         gCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGAAGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcggcggcgg                                                                              >  2:98330/1‑150 (MQ=255)
                                                                          ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                             <  2:407281/150‑1 (MQ=255)
                                                                           cGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGAt                                                                            >  2:1352264/1‑150 (MQ=255)
                                                                           cGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGAt                                                                            <  2:417436/150‑1 (MQ=255)
                                                                             cATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATAc                                                                          >  1:1138569/1‑150 (MQ=255)
                                                                             cATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATAc                                                                          >  1:1447576/1‑150 (MQ=255)
                                                                                gCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGTAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGTTCGGACGGCGGCGGATACGGt                                                                       <  1:214780/150‑1 (MQ=255)
                                                                                 ccAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTc                                                                      >  2:892740/1‑150 (MQ=255)
                                                                                       aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                                >  2:1228752/1‑150 (MQ=255)
                                                                                          atatTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGAt                                                             <  1:1543669/150‑1 (MQ=255)
                                                                                           tatTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATg                                                            <  2:553406/150‑1 (MQ=255)
                                                                                             ttGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGAc                                                          <  1:702277/150‑1 (MQ=255)
                                                                                              tGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACg                                                         <  1:448496/150‑1 (MQ=255)
                                                                                              tGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACg                                                         >  2:1288304/1‑150 (MQ=255)
                                                                                               gcgcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGc                                                        >  2:189634/1‑150 (MQ=255)
                                                                                                      atagCCAATAAAAACCTAGTCGAAAACCCGAAGATACGCTCGAAGCACAAAACGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATc                                                 >  2:1083863/4‑150 (MQ=255)
                                                                                                       ttGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATcc                                                >  1:726935/1‑150 (MQ=255)
                                                                                                          ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                             <  2:201939/150‑1 (MQ=255)
                                                                                                          ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                             >  1:631862/1‑150 (MQ=255)
                                                                                                            aTTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGTGTAGATGACGCTCACATCCGGGtt                                           <  2:533734/150‑1 (MQ=255)
                                                                                                            aTTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGtt                                           <  2:484969/150‑1 (MQ=255)
                                                                                                             ttAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTa                                          <  1:1228752/150‑1 (MQ=255)
                                                                                                              tAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGa                                                              >  1:550939/1‑129 (MQ=255)
                                                                                                              tAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGa                                                              <  2:550939/129‑1 (MQ=255)
                                                                                                               aaaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTaaa                                        <  1:1000936/150‑1 (MQ=255)
                                                                                                                      ttGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGc                                 >  2:1281658/1‑150 (MQ=255)
                                                                                                                        ggCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTGGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAACCCATTGCtt                               <  2:931222/150‑1 (MQ=255)
                                                                                                                          cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                             >  2:1117269/1‑150 (MQ=255)
                                                                                                                           gAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTcc                            <  1:896529/150‑1 (MQ=255)
                                                                                                                            aaaaCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCa                           >  2:3670/1‑150 (MQ=255)
                                                                                                                              aaCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATc                         <  1:713769/150‑1 (MQ=255)
                                                                                                                               aCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATg                                                            <  2:1074705/114‑1 (MQ=255)
                                                                                                                                cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGCGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                       >  2:145602/1‑150 (MQ=255)
                                                                                                                                    cTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGGCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTc                   <  2:321982/150‑1 (MQ=255)
                                                                                                                                     tggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCa                  >  1:1639345/5‑150 (MQ=255)
                                                                                                                                      ggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCaa                 >  2:391041/4‑150 (MQ=255)
                                                                                                                                      gATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCaa                 >  1:941799/1‑150 (MQ=255)
                                                                                                                                        tACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAAt               >  2:83589/1‑150 (MQ=255)
                                                                                                                                         aCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATATGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg              <  1:98330/150‑1 (MQ=255)
                                                                                                                                          cGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTATATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg              <  1:1451715/149‑1 (MQ=255)
                                                                                                                                          cGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTATATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg              >  2:1451715/1‑149 (MQ=255)
                                                                                                                                          cGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGt             >  2:1547672/1‑150 (MQ=255)
                                                                                                                                               gAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCcaca        <  1:213819/150‑1 (MQ=255)
                                                                                                                                               gAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCcaca        <  1:1117269/150‑1 (MQ=255)
                                                                                                                                                aaGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAg       >  2:1655760/1‑150 (MQ=255)
                                                                                                                                                  gcagcCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGAc     <  1:1667141/150‑1 (MQ=255)
                                                                                                                                                    agccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGa   >  2:1624278/1‑150 (MQ=255)
                                                                                                                                                    agccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGa   >  1:1026078/1‑150 (MQ=255)
                                                                                                                                                     gccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAt  >  1:1060033/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGAT  >  NC_000913/858276‑858574

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