Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,013 C→A 5.0% intergenic (+28/+23) 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,0130CA 5.0% 312.3 / 4.1 99intergenic (+28/+23)csdE/tcdACsdA‑binding activator; Fe‑S protein/tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA
Reads supporting (aligned to +/- strand):  ref base C (51/42);  new base A (2/3);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 103 total aligned reads displayed.

GGCACAGTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGG  >  NC_000913/2944865‑2945158
                                                                                                                                                    |                                                                                                                                                 
ggCACAGTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTcg                                                                                                                                                  >  2:455285/1‑150 (MQ=255)
ggCACAGTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTcg                                                                                                                                                  >  2:732182/1‑150 (MQ=255)
ggCACAGTCACCAATGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTcg                                                                                                                                                  >  2:443660/1‑150 (MQ=255)
  cacaGTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                                                                                                                                                <  1:732182/150‑1 (MQ=255)
        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                                                                                                                                      >  2:6848/1‑150 (MQ=255)
            acTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGa                                                                                                                                      <  1:415091/150‑1 (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:1123921/1‑150 (MQ=255)
                gCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg                                                                                                                                  >  1:737792/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                                                                                                                             <  1:922715/150‑1 (MQ=255)
                     gTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt                                                                                                                             >  2:25571/1‑150 (MQ=255)
                     gTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt                                                                                                                             >  2:311957/1‑150 (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                                                                                                    >  1:520922/1‑141 (MQ=255)
                                                       gcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                    <  2:520922/141‑1 (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)
                                                                        gTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAAc                                                                          >  1:721786/1‑150 (MQ=255)
                                                                        gTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGAAAAGACGATCGCGTCTCATCCGCCCTGATTGTTAAACCCTGCCGCGCCGCTTTCGCCAACATCTTCTTCAGCGCATGAGAAAc                                                                          >  1:307810/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:135733/150‑1 (MQ=255)
                                                                                              aTTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCg                                                    <  2:1123921/150‑1 (MQ=255)
                                                                                              aTTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCg                                                    <  2:737792/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                                     <  1:25571/150‑1 (MQ=255)
                                                                                                             aaGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGcc                                     >  2:1652852/1‑150 (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)
                                                                                                                      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:1275681/5‑150 (MQ=39)
                                                                                                                                atcaggTCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGcac                  >  1:773346/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)
                                                                                                                                            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)
                                                                                                                                                    |                                                                                                                                                 
GGCACAGTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGG  >  NC_000913/2944865‑2945158

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