Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,000,308:1 +G 7.7% coding (749/765 nt) yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,000,3081.G 7.7% 377.7 / 28.1 91coding (749/765 nt)yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base . (42/42);  new base G (2/5);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTG  >  NC_000913/4000167‑4000455
                                                                                                                                              |                                                                                                                                                   
gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaag                                                                                                                                              <  2:478611/149‑2 (MQ=255)
gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaag                                                                                                                                              <  4:196928/149‑2 (MQ=255)
 cACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaagg                                                                                                                                             <  2:69647/149‑3 (MQ=255)
   ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                            <  4:564873/148‑4 (MQ=255)
   ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                            <  4:389684/148‑4 (MQ=255)
   ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                            >  3:564873/1‑145 (MQ=255)
   ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggt                                                                                                                                            >  3:389684/1‑145 (MQ=255)
    cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                          >  1:511791/1‑149 (MQ=255)
    cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                          >  3:264134/1‑149 (MQ=255)
    cTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                                                                                          >  4:120402/1‑149 (MQ=255)
      tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGt                                                                                                                                        <  2:99290/149‑1 (MQ=255)
      tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGt                                                                                                                                        <  1:74992/149‑1 (MQ=255)
        tttGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATAATGGTGa                                                                                                                                      <  2:313945/149‑1 (MQ=255)
         ttGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGAc                                                                                                                                     >  2:29245/1‑149 (MQ=255)
          tGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACa                                                                                                                                    >  3:289121/1‑149 (MQ=255)
            gCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACaaa                                                                                                                                  >  1:310005/1‑149 (MQ=255)
             ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACaaaa                                                                                                                                 >  1:521552/1‑149 (MQ=255)
              cTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                                                                                                <  3:143320/149‑1 (MQ=255)
                 aaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGAt                                                                                                                             >  4:344891/1‑149 (MQ=255)
                 aaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGAt                                                                                                                             >  2:40154/1‑149 (MQ=255)
                         ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGc                                                                                                                     >  2:133586/1‑149 (MQ=255)
                         ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGc                                                                                                                     <  1:53446/149‑1 (MQ=255)
                          tCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                                                                                                    <  2:261807/149‑1 (MQ=255)
                              aTTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACgg                                                                                                                <  4:434231/149‑1 (MQ=255)
                                 cAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAt                                                                                                             >  2:120514/1‑149 (MQ=255)
                                   aTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                                                                                           >  1:90551/1‑149 (MQ=255)
                                    tAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc                                                                                                          >  1:265797/1‑149 (MQ=255)
                                      aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                                                                                        <  4:197127/149‑1 (MQ=255)
                                      aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                                                                                        <  4:467905/149‑1 (MQ=255)
                                       aTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                                                                                       <  1:368117/149‑1 (MQ=255)
                                       aTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                                                                                       <  3:47795/149‑1 (MQ=255)
                                           cAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCt                                                                                                   <  4:264134/149‑1 (MQ=255)
                                              cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                                                                                <  3:545987/149‑1 (MQ=255)
                                               gTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                                                                                               <  1:120514/149‑1 (MQ=255)
                                                 cATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCata                                                                                             <  1:179455/149‑1 (MQ=255)
                                                        gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                      <  3:495760/149‑1 (MQ=255)
                                                        gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                      <  2:427908/149‑1 (MQ=255)
                                                        gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                      <  2:343244/149‑1 (MQ=255)
                                                        gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACatat                                                                                      <  1:285542/149‑1 (MQ=255)
                                                          cGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt                                                                                    <  3:344891/149‑1 (MQ=255)
                                                            gATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGa                                                                                  >  1:510676/1‑149 (MQ=255)
                                                            gATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGa                                                                                  >  4:479540/1‑149 (MQ=255)
                                                                gggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATc                                                                              <  2:521552/149‑1 (MQ=255)
                                                                 ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                             >  1:523565/1‑149 (MQ=255)
                                                                  gCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt                                                                            >  4:266094/1‑149 (MQ=255)
                                                                    gTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGc                                                                          <  2:255240/149‑1 (MQ=255)
                                                                             gcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAt                                                                 >  2:571462/1‑149 (MQ=255)
                                                                              cATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATa                                                                <  2:129851/149‑1 (MQ=255)
                                                                                   ggCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAagca                                                           >  4:514280/1‑149 (MQ=255)
                                                                                    gCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATaaa                                                              >  4:433396/1‑145 (MQ=255)
                                                                                    gCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATaaa                                                              <  3:433396/145‑1 (MQ=255)
                                                                                         tGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAgctgct                                                     >  4:66811/1‑149 (MQ=255)
                                                                                               aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                             >  4:115301/1‑119 (MQ=255)
                                                                                               aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                             >  4:177722/1‑119 (MQ=255)
                                                                                               aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                             <  3:115301/119‑1 (MQ=255)
                                                                                               aTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg                                                                             <  3:177722/119‑1 (MQ=255)
                                                                                                 cTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTc                                             <  2:101193/149‑1 (MQ=255)
                                                                                                  ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                            <  1:245198/149‑1 (MQ=255)
                                                                                                  ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                            >  3:17583/1‑149 (MQ=255)
                                                                                                  ttCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                            >  1:392675/1‑149 (MQ=255)
                                                                                                     aGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATCCATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGtt                                         <  4:243913/149‑1 (MQ=255)
                                                                                                     aGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGtt                                         <  4:17583/149‑1 (MQ=255)
                                                                                                      gTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGttt                                        <  1:45736/149‑1 (MQ=255)
                                                                                                      gTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGttt                                        >  1:363609/1‑149 (MQ=255)
                                                                                                       tATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                       <  1:395009/149‑1 (MQ=255)
                                                                                                        aTCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCg                                      >  1:106530/1‑149 (MQ=255)
                                                                                                           tGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCtt                                   >  3:123221/1‑149 (MQ=255)
                                                                                                            gACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCttt                                  <  1:127838/149‑1 (MQ=255)
                                                                                                             acacAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTg                                 >  1:63287/1‑149 (MQ=255)
                                                                                                                 aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                             >  3:534292/1‑149 (MQ=255)
                                                                                                                 aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                             <  3:529566/149‑1 (MQ=255)
                                                                                                                 aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                             >  1:282709/1‑149 (MQ=255)
                                                                                                                 aaaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCAt                             >  3:155080/1‑149 (MQ=255)
                                                                                                                  aaaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATa                            >  3:92697/1‑149 (MQ=255)
                                                                                                                   aaCTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATaa                           <  2:523565/149‑1 (MQ=255)
                                                                                                                     cTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACa                         <  1:492898/149‑1 (MQ=255)
                                                                                                                       aTCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAgg                       >  4:18900/1‑149 (MQ=255)
                                                                                                                            tttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCaaaa                  <  2:106530/149‑1 (MQ=255)
                                                                                                                            tttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCaaaa                  <  2:510676/149‑1 (MQ=255)
                                                                                                                             ttAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAAt                 >  1:126732/1‑149 (MQ=255)
                                                                                                                                 ccTTTCCGCTCGA‑CGGATATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATa             >  1:24436/1‑149 (MQ=255)
                                                                                                                                 ccTGTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATa             >  4:299912/1‑149 (MQ=255)
                                                                                                                                  cTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                       >  4:354680/1‑122 (MQ=255)
                                                                                                                                  cTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTc                                       <  3:354680/122‑1 (MQ=255)
                                                                                                                                   tttCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATaaa           >  1:401997/1‑149 (MQ=255)
                                                                                                                                    ttCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATaaaa          <  3:479540/149‑1 (MQ=255)
                                                                                                                                       cGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAAtt       <  2:63287/149‑1 (MQ=255)
                                                                                                                                         cTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAAc                                                 <  1:291310/105‑1 (MQ=255)
                                                                                                                                         cTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAAc                                                 >  2:291310/1‑105 (MQ=255)
                                                                                                                                            gA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATttgttg  >  4:183084/1‑149 (MQ=255)
                                                                                                                                            gA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATttgttg  >  2:363320/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTG  >  NC_000913/4000167‑4000455

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