Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,074 T→C 6.2% intergenic (‑71/+43) rna ← / ← citT RNase I/citrate:succinate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,0740TC 6.2% 240.3 / 5.2 80intergenic (‑71/+43)rna/citTRNase I/citrate:succinate antiporter
Reads supporting (aligned to +/- strand):  ref base T (27/48);  new base C (2/3);  total (29/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACAT  >  NC_000913/644926‑645220
                                                                                                                                                    |                                                                                                                                                  
cAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTa                                                                                                                                                   >  2:894923/1‑150 (MQ=255)
  aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                                 <  2:404686/150‑1 (MQ=255)
     cTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTgg                                                                                                                                              >  2:408397/1‑150 (MQ=255)
       aaGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcg                                                                                                                                            <  1:1483198/150‑1 (MQ=255)
        aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                           <  2:1188915/150‑1 (MQ=255)
        aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                           >  2:609171/1‑150 (MQ=255)
        aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgg                                                                                                                                           <  2:1657094/150‑1 (MQ=255)
         ggCGTTGACAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcggg                                                                                                                                          <  1:1167011/150‑1 (MQ=255)
          gCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgggc                                                                                                                                         <  2:1024669/150‑1 (MQ=255)
              tGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtga                                                                                                                                     <  2:600938/150‑1 (MQ=255)
               ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgat                                                                                                                                    <  1:1084404/150‑1 (MQ=255)
               ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgat                                                                                                                                    <  2:145068/150‑1 (MQ=255)
                  agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                                 >  2:111507/1‑150 (MQ=255)
                      gagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAt                                                                                                                             <  1:627743/150‑1 (MQ=255)
                      gagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAt                                                                                                                             <  1:1078192/150‑1 (MQ=255)
                        gaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                           <  2:1138487/150‑1 (MQ=255)
                         aagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCa                                                                                                                          <  2:737796/150‑1 (MQ=255)
                            ggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCg                                                                                                                       <  2:733483/150‑1 (MQ=255)
                             gAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGt                                                                                                                      <  1:1415223/150‑1 (MQ=255)
                                 cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCATTTGGGCTTACCTGGCTGTCGTGATTATTTATTCAGCGTTTgg                                                                                                                  <  2:149521/150‑1 (MQ=255)
                                  agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                                 >  2:712338/1‑150 (MQ=255)
                                  agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                                 >  1:1055831/1‑150 (MQ=255)
                                    agaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAg                                                                                                                         >  1:988157/1‑140 (MQ=255)
                                    agaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAg                                                                                                                         <  2:988157/140‑1 (MQ=255)
                                    agaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGa                                                                                                               <  1:1495314/150‑1 (MQ=255)
                                       aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                            <  2:1248846/150‑1 (MQ=255)
                                        aCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGt                                                                                                           <  2:414606/150‑1 (MQ=255)
                                         ccGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTa                                                                                                          >  1:436940/1‑150 (MQ=255)
                                          cgcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                         <  1:916766/150‑1 (MQ=255)
                                          cgcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                         <  1:922566/150‑1 (MQ=255)
                                             gAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAg                                                                                                      >  1:597997/1‑150 (MQ=255)
                                               gCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGtt                                                                                                    >  2:649538/1‑150 (MQ=255)
                                                 aaCGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTcc                                                                                                  <  2:469307/150‑1 (MQ=255)
                                                 aaCGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTcc                                                                                                  <  2:138475/150‑1 (MQ=255)
                                                 aaCGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACAATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTcc                                                                                                  <  2:1498532/150‑1 (MQ=255)
                                                    gcggcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcaca                                                                                               >  2:1209478/1‑150 (MQ=255)
                                                         gTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCg                                                                                          <  2:449645/150‑1 (MQ=255)
                                                         gTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCg                                                                                          <  2:549804/150‑1 (MQ=255)
                                                         gTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCg                                                                                          <  2:555466/150‑1 (MQ=255)
                                                            aCGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCgaga                                                                                       >  1:1036244/1‑150 (MQ=255)
                                                             cGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGaa                                                                                      >  1:717001/1‑150 (MQ=255)
                                                               ccAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATc                                                                                    <  2:1634444/150‑1 (MQ=255)
                                                                 aGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCgg                                                                                  >  2:453107/1‑150 (MQ=255)
                                                                  gAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgc                                                                                 <  2:186108/150‑1 (MQ=255)
                                                                  gAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgc                                                                                 <  1:894923/150‑1 (MQ=255)
                                                                     tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                              <  2:473022/150‑1 (MQ=255)
                                                                           cATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaacc                                                                        <  2:1370243/150‑1 (MQ=255)
                                                                            aTAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaacca                                                                       <  1:980125/150‑1 (MQ=255)
                                                                              aaTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCcaaccaac                                                                     <  1:798555/150‑1 (MQ=255)
                                                                                  ggTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCGAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAAGCAACAGCa                                                                 <  2:717001/150‑1 (MQ=255)
                                                                                  ggTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCa                                                                 <  1:111507/150‑1 (MQ=255)
                                                                                   gtgtGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCaa                                                                >  2:1524974/1‑150 (MQ=255)
                                                                                    tgtgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGGAAc                                                               <  1:453107/150‑1 (MQ=255)
                                                                                    tgtgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAAc                                                               <  2:436940/150‑1 (MQ=255)
                                                                                     gtgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACa                                                              <  1:1586248/150‑1 (MQ=255)
                                                                                      tgtaaCTCATACATACACTGAATACTATCTATTAAATCATAAAAACGCCCGCCAGGTGGGCTTACGTGGCGGGCGTGATGATTTATTCAGCGTTTGGCTAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACCGCAACAt                                                             <  2:1603896/147‑1 (MQ=255)
                                                                                          aCTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATagag                                                         <  2:1486967/150‑1 (MQ=255)
                                                                                           cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGAGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAgaga                                                        >  2:735002/1‑150 (MQ=255)
                                                                                             catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                      >  2:454540/1‑150 (MQ=255)
                                                                                             catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                      >  2:492599/1‑150 (MQ=255)
                                                                                                   tacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATc                                                <  1:169507/150‑1 (MQ=255)
                                                                                                     cacaGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAc                                              >  1:746528/1‑150 (MQ=255)
                                                                                                      aatGAAAACAATCTATTAAATCATAAAAAAGCCCGCCAGCTGGCCTTACATGCCGGGCCAGATGATTTATTAAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcc                                             >  2:543727/3‑150 (MQ=255)
                                                                                                       cTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAccc                                            >  2:967917/1‑150 (MQ=255)
                                                                                                        tGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcccc                                           >  2:973621/1‑150 (MQ=255)
                                                                                                             aCTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGa                                      <  1:713102/150‑1 (MQ=255)
                                                                                                              ctatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGAt                                     >  2:1116815/1‑150 (MQ=255)
                                                                                                                   tatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcgc                                <  1:147022/150‑1 (MQ=255)
                                                                                                                      tAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCaa                             >  1:1457612/1‑150 (MQ=255)
                                                                                                                      tAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCaa                             >  2:202554/1‑150 (MQ=255)
                                                                                                                           cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATCGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGc                        >  2:870255/1‑150 (MQ=255)
                                                                                                                            aTAAAAAAGCCCGCCAGGTGGGCTTACCGGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                       <  2:597997/150‑1 (MQ=255)
                                                                                                                                 cacgCCCGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTaa                  >  1:392629/4‑150 (MQ=255)
                                                                                                                                     cccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt              >  1:1017910/1‑150 (MQ=255)
                                                                                                                                     cccGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt              >  2:1454395/1‑150 (MQ=255)
                                                                                                                                         ccAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGa          <  1:1524974/150‑1 (MQ=255)
                                                                                                                                             gTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTc      <  1:454540/150‑1 (MQ=255)
                                                                                                                                               gggCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTcac    <  1:164599/150‑1 (MQ=255)
                                                                                                                                               gggCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTcac    <  1:376998/150‑1 (MQ=255)
                                                                                                                                                 gCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACAt  <  1:649538/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACAT  >  NC_000913/644926‑645220

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