Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,020 C→T 5.0% intergenic (+35/+16) csdE → / ← tcdA CsdA‑binding activator; Fe‑S protein/tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,945,0200CT 5.0% 324.4 / 3.5 100intergenic (+35/+16)csdE/tcdACsdA‑binding activator; Fe‑S protein/tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA
Reads supporting (aligned to +/- strand):  ref base C (52/42);  new base T (2/3);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 105 total aligned reads displayed.

CACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCG  >  NC_000913/2944873‑2945169
                                                                                                                                                   |                                                                                                                                                     
caccacTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCAt                                                                                                                                                     >  2:53537/1‑150 (MQ=255)
  ccacTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGTGTCGGATAAGACGTTCGCGTCTCATcc                                                                                                                                                   <  2:1235416/150‑1 (MQ=255)
  ccacTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATcc                                                                                                                                                   >  2:1387834/1‑150 (MQ=255)
   cacTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCg                                                                                                                                                  <  2:58808/150‑1 (MQ=255)
    acTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGa                                                                                                                                                 <  1:415091/150‑1 (MQ=255)
    acTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGa                                                                                                                                                 >  2:6848/1‑150 (MQ=255)
     cTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCGAc                                                                                                                                                <  1:810299/150‑1 (MQ=255)
      tGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGAcc                                                                                                                                               <  1:963849/150‑1 (MQ=255)
       ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAATCCGTTCGCGTCTCATACGAcca                                                                                                                                              >  1:1339751/1‑149 (MQ=255)
       ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                                                                                                                                                  <  1:438640/146‑1 (MQ=255)
       ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                                                                                                                                                  >  2:438640/1‑146 (MQ=255)
       ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt                                                                                                                                              <  2:1630678/150‑1 (MQ=255)
        gCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg                                                                                                                                             >  1:737792/1‑150 (MQ=255)
        gCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg                                                                                                                                             >  1:1123921/1‑150 (MQ=255)
        gCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGGCGGATGAGACGCGAACGTCTTATCCGAcccg                                                                                                                                             <  1:171106/150‑3 (MQ=255)
         cATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGa                                                                                                                                            >  1:1057428/1‑150 (MQ=255)
          aTTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGAt                                                                                                                                           >  2:1152571/1‑150 (MQ=255)
             gTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt                                                                                                                                        >  2:311957/1‑150 (MQ=255)
             gTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt                                                                                                                                        <  1:922715/150‑1 (MQ=255)
                    taGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAccc                                                                                                                                 <  2:1651093/149‑1 (MQ=255)
                    gaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAccc                                                                                                                                 >  2:905126/1‑150 (MQ=255)
                    gaGCTGGGAGTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAccc                                                                                                                                 <  1:1087450/150‑1 (MQ=255)
                          ggATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGAcgc                                                                                                                           >  1:1214961/1‑150 (MQ=255)
                              tACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGcgccg                                                                                                                       >  2:617423/1‑150 (MQ=255)
                               aCGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGcgccgc                                                                                                                      >  2:242178/1‑150 (MQ=255)
                                  tGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCttt                                                                                                                   <  1:311957/150‑1 (MQ=255)
                                   gcgcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCAATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTc                                                                                                                  >  2:1398305/1‑150 (MQ=255)
                                    cgcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCg                                                                                                                 <  2:1457055/150‑1 (MQ=255)
                                     gcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGc                                                                                                                <  1:856692/150‑1 (MQ=255)
                                          ttAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatca                                                                                                           <  1:1152571/150‑1 (MQ=255)
                                           tAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatcat                                                                                                          >  1:370367/1‑150 (MQ=255)
                                            aGCGCCTCCCGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGAGCCGCTTTCGCAatcatc                                                                                                         >  2:1006974/1‑150 (MQ=255)
                                            aGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCatcatc                                                                                                         <  1:739866/150‑1 (MQ=255)
                                             gcgcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGGCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtct                                                                                                        <  2:1133005/150‑1 (MQ=255)
                                               gcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                               <  2:520922/141‑1 (MQ=255)
                                               gcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                               >  1:520922/1‑141 (MQ=255)
                                                ccgcACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGGTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtcttct                                                                                                     <  2:175828/147‑1 (MQ=255)
                                                ccTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtcttct                                                                                                     <  2:1214961/150‑1 (MQ=255)
                                                 cTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATcttctt                                                                                                    >  2:1478058/1‑150 (MQ=255)
                                                   cACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCa                                                                                                  <  1:765936/150‑1 (MQ=255)
                                                   cACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCa                                                                                                  >  2:718945/1‑150 (MQ=255)
                                                         gCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCAt                                                                                            <  1:1554719/150‑1 (MQ=255)
                                                            aGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATgag                                                                                         >  2:930449/1‑150 (MQ=255)
                                                              gggTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGaa                                                                                       >  2:1396268/1‑150 (MQ=255)
                                                              gggTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCATTCGCCATCATCTTCTTCAGCGCATGAGaa                                                                                       >  1:565207/1‑150 (MQ=255)
                                                                 ttAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAAcc                                                                                    <  1:1660018/150‑1 (MQ=255)
                                                                 ttAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAAcc                                                                                    <  2:722788/150‑1 (MQ=255)
                                                                   aaaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGc                                                                                  <  1:728698/150‑1 (MQ=255)
                                                                   aaaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGc                                                                                  <  2:1491089/150‑1 (MQ=255)
                                                                   aaaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGc                                                                                  <  1:6848/150‑1 (MQ=255)
                                                                      tGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAGCCGcaac                                                                               <  1:455285/150‑1 (MQ=255)
                                                                       gCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGcaaca                                                                              <  1:905126/150‑1 (MQ=255)
                                                                         gTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACaaa                                                                            >  1:650218/1‑150 (MQ=255)
                                                                           tAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTCCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAAcc                                                                          <  1:537793/150‑1 (MQ=255)
                                                                             aGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCa                                                                        >  2:1690559/1‑150 (MQ=255)
                                                                               cGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCaaa                                                                      >  2:932191/1‑150 (MQ=255)
                                                                                gAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAg                                                                     >  1:1424931/1‑150 (MQ=255)
                                                                                 aGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAgg                                                                    <  1:415360/150‑1 (MQ=255)
                                                                                    cGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTgg                                                                 <  2:1057428/150‑1 (MQ=255)
                                                                                     gATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGc                                                                >  2:1713419/1‑150 (MQ=255)
                                                                                      aTTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCg                                                               <  2:737792/150‑1 (MQ=255)
                                                                                      aTTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCg                                                               <  2:135733/150‑1 (MQ=255)
                                                                                      aTTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCg                                                               <  2:1123921/150‑1 (MQ=255)
                                                                                         ttCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTc                                                            <  2:795437/149‑1 (MQ=255)
                                                                                          tCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCa                                                           <  1:1479215/150‑1 (MQ=255)
                                                                                           cGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCAc                                                          <  2:633639/150‑1 (MQ=255)
                                                                                            gCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCCTCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCAcc                                                         >  1:442925/1‑150 (MQ=255)
                                                                                            gCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCAcc                                                         <  1:443660/150‑1 (MQ=255)
                                                                                             cTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCa                                                        >  2:1330705/1‑150 (MQ=255)
                                                                                               gCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGAACTGATTGTTTAACCCTGACGCGCCGCATTCGCCATCATCTTATTCAGCGCAAGAGAAACCGCAACAAAAACAAATGTGGCGGTCACCATc                                                      >  1:602559/1‑150 (MQ=255)
                                                                                                  acgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGtt                                                   >  1:456945/1‑150 (MQ=255)
                                                                                                     aaGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGcc                                                >  2:1652852/1‑150 (MQ=255)
                                                                                                     aaGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGcc                                                <  1:25571/150‑1 (MQ=255)
                                                                                                      aGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGGCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCcg                                               >  1:1382621/1‑150 (MQ=255)
                                                                                                        cAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCcgcg                                             <  2:705259/150‑1 (MQ=255)
                                                                                                          ggTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGcc                                           <  2:270450/150‑1 (MQ=255)
                                                                                                           gTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACATGATTGTTTAACCCTGACGCGCCGCTTTCGCCAGCATCTTCTTCAGCGCATGAGAAACCGCGACAAGCCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCa                                          >  1:326576/1‑150 (MQ=255)
                                                                                                           gTTTAAGCACAGAACCGGGTCAGTTAATACGATGGCGTCTAAAACCACCTGATTTTTTAAACCAGACGCGCCGCTTTACCCCTAATCTTCTTAAGCCCATGAGAACACGCAACAAACCCAAAAGGGGCGGTGACCAAAGTTTGCGCGCCa                                          >  2:1526099/1‑150 (MQ=255)
                                                                                                              tAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCCGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCaaag                                       >  2:179956/1‑149 (MQ=255)
                                                                                                                 gCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCACAGGTGGCGGTCCCCATCGGTGCCGCGCCAAAGCCa                                    >  2:1036869/1‑150 (MQ=255)
                                                                                                                 gCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCa                                    <  2:274361/150‑1 (MQ=255)
                                                                                                                  cacaGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAg                                   >  1:653488/1‑150 (MQ=255)
                                                                                                                       aaCCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGca                              >  1:876638/1‑150 (MQ=255)
                                                                                                                        atcaggTCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGcac                             >  1:773346/5‑150 (MQ=39)
                                                                                                                        atcaggTCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGcac                             >  1:1275681/5‑150 (MQ=39)
                                                                                                                        aCCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAGCTTCTTCAGCGCAAGAGAAAAAGCAACAAACCCCAAGGTGGCGGTCACCGTCGGTGCCGCGCCAAAGCCAGATGcgc                             >  2:1597590/1‑148 (MQ=255)
                                                                                                                          cGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACaa                           >  1:1673820/1‑150 (MQ=255)
                                                                                                                           gggTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAat                          >  2:422696/1‑150 (MQ=255)
                                                                                                                             gTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCACTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatcc                        >  1:53539/1‑150 (MQ=255)
                                                                                                                             gTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCACTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatcc                        >  1:1483979/1‑150 (MQ=255)
                                                                                                                                ggATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGAGGCACAatccatc                     >  1:586992/1‑150 (MQ=255)
                                                                                                                                    aaGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAACCCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCt                 >  1:1707378/1‑150 (MQ=255)
                                                                                                                                     aGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCtt                >  1:176860/1‑150 (MQ=255)
                                                                                                                                       aCGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCg              <  2:370367/150‑1 (MQ=255)
                                                                                                                                        cGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCgg             >  1:1214297/1‑150 (MQ=255)
                                                                                                                                               gTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCt      >  1:345475/1‑150 (MQ=255)
                                                                                                                                                tctcATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTg     >  2:1385982/1‑150 (MQ=255)
                                                                                                                                                  tcATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGcc   >  1:1656312/1‑150 (MQ=255)
                                                                                                                                                   cATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCg  >  2:1053901/1‑150 (MQ=255)
                                                                                                                                                   cATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCGTCTGCCg  >  2:1151295/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCG  >  NC_000913/2944873‑2945169

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