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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,0740TC 6.5% 370.3 / 11.5 123intergenic (‑71/+43)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base T (46/69);  new base C (2/6);  total (48/75)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 125 total aligned reads displayed.

GCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACATAGC  >  NC_000913/644925‑645223
                                                                                                                                                     |                                                                                                                                                     
gCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCtt                                                                                                                                                       >  2:1178549/1‑150 (MQ=255)
  aaCGCTAAGGCGTTGGCAGAAGAGAAGGGGAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAc                                                                                                                                                     <  2:1661871/150‑1 (MQ=255)
   aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                                    <  1:651168/150‑1 (MQ=255)
   aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                                    >  1:1587585/1‑150 (MQ=255)
    cGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                                   >  1:969546/1‑150 (MQ=255)
    cGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                                   >  1:1635319/1‑150 (MQ=255)
      cTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTgg                                                                                                                                                 >  1:1394495/1‑150 (MQ=255)
        aaGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcg                                                                                                                                               <  2:1873832/150‑1 (MQ=255)
         aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                              <  1:132403/150‑1 (MQ=255)
         aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                              >  2:654337/1‑150 (MQ=255)
          ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcggg                                                                                                                                             <  1:244166/150‑1 (MQ=255)
          ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcggg                                                                                                                                             <  2:1268236/150‑1 (MQ=255)
          cgCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGCGCCACGCGGCGGTACGCCAGACTGCTTGCATACTGGGTGGGGAACGCATACCTCCACTGAATACTATCTATTAACTCCTAAAACAGCCCGCCAGGTGGGCTTACCTggcggg                                                                                                                                             <  1:1169190/149‑1 (MQ=255)
             gTTGGCAGACGACAAGGGAAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATCCTTACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                          <  1:1473755/150‑1 (MQ=255)
                ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgat                                                                                                                                       <  1:1225523/150‑1 (MQ=255)
                ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAAAGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgat                                                                                                                                       >  1:1376297/1‑150 (MQ=255)
                  cAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgggctttttta                                                                                                                                     <  2:771810/150‑8 (MQ=255)
                  cAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgggctttttta                                                                                                                                     <  1:1105386/150‑8 (MQ=255)
                   agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                                    <  2:996977/150‑1 (MQ=255)
                   agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                                    >  2:895070/1‑150 (MQ=255)
                   agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATAAACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                                    <  2:1797365/150‑1 (MQ=255)
                          aagGGCAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGACTGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCa                                                                                                                             <  1:1699164/150‑1 (MQ=255)
                            gggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGCCTTACCTGGCGGGCGTGATGATTTATTCacc                                                                                                                           >  2:728571/1‑148 (MQ=255)
                             ggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCg                                                                                                                          <  2:1587585/150‑1 (MQ=255)
                              gAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGt                                                                                                                         <  2:1656975/150‑1 (MQ=255)
                                  cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                                     <  1:504724/150‑1 (MQ=255)
                                   agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                                    <  2:456384/150‑1 (MQ=255)
                                    gagaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAATAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTAATCGTTTGGCg                                                                                                                   >  2:1819011/1‑150 (MQ=255)
                                    gagaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCg                                                                                                                   >  1:662588/1‑150 (MQ=255)
                                      gaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGaa                                                                                                                 >  2:396687/1‑150 (MQ=255)
                                       aaaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAAc                                                                                                                <  1:541251/150‑1 (MQ=255)
                                        aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                               <  2:1320850/150‑1 (MQ=255)
                                        aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                               <  1:1081825/150‑1 (MQ=255)
                                        aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                               <  2:1636312/150‑1 (MQ=255)
                                          ccGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTa                                                                                                             <  1:1823361/150‑1 (MQ=255)
                                            gcgAGCAACGGGGGGTTACGCCAGAATGCTTTGATAATGGGGGTGGATCTCATACATACACTGAATACTATCTATTAACTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                           <  2:1618672/150‑1 (MQ=255)
                                            gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATACTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                           <  2:1532380/150‑1 (MQ=255)
                                            gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                           <  1:1463974/150‑1 (MQ=255)
                                            gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACAAGCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                           <  2:1256773/150‑1 (MQ=255)
                                             cgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACTCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTa                                                                                                          <  1:242765/150‑1 (MQ=255)
                                               aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGTCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                        >  1:1536470/1‑150 (MQ=255)
                                               aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                        >  2:1037392/1‑150 (MQ=255)
                                               aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                        >  2:1123897/1‑150 (MQ=255)
                                               aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                        >  2:191788/1‑150 (MQ=255)
                                                 cAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACCTAATAGtta                                                                                                      >  2:1551236/1‑149 (MQ=255)
                                                  caCGCGGAGTTACGCCAGAATGCTTTCATAATTGGTGTGGAACTCATACATACACTGAATACTAGCTCTTAAATCATAAAAAAGCCCGCCCGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTCGTTcc                                                                                                     <  2:916827/149‑1 (MQ=255)
                                                    cgcgGCGTTACGCCAGAATGCTTTCATAATGGGTGTAGAACTCATACATACACTGAATACTATCTTTTATATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCTGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcac                                                                                                   >  2:1907126/1‑150 (MQ=255)
                                                     gcggcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcaca                                                                                                  <  2:970404/150‑1 (MQ=255)
                                                        gcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATgg                                                                                               <  1:1043162/150‑1 (MQ=255)
                                                            tACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCgag                                                                                           >  1:1652064/1‑150 (MQ=255)
                                                             aCGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCCCATGGCgaga                                                                                          >  2:1630365/1‑150 (MQ=255)
                                                                  aGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGGGTGATGATTTATTCAGAGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATAgg                                                                                     >  2:64463/1‑150 (MQ=255)
                                                                  aGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCgg                                                                                     >  1:649985/1‑150 (MQ=255)
                                                                  aGAATGCTTTCATAATGGGTGTGGAACTCATACATAAACTGATAACTATCTATTAAATCATAAAAAAGCCCGACAGGTGGGCTTAACTGGCGGGAGTGATGATTTATTCAGCGTTTGGCGAACTTATTAGTTCCACATGGCGAGAATctg                                                                                     >  1:668903/1‑148 (MQ=255)
                                                                      tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                                 <  1:1713824/150‑1 (MQ=255)
                                                                      tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                                 <  2:842273/150‑1 (MQ=255)
                                                                      tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                                 <  1:584566/150‑1 (MQ=255)
                                                                      tGCTTTCATAATGGGGGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                                 <  1:814174/150‑1 (MQ=255)
                                                                        cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                               <  1:370995/150‑1 (MQ=255)
                                                                        cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                               >  1:785859/1‑150 (MQ=255)
                                                                        cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                               >  1:1750469/1‑150 (MQ=255)
                                                                        cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCAAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                               >  1:536587/1‑150 (MQ=255)
                                                                             aTAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATAATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTATGGCGAACGTATTAGTTCCACATGGCGAGAATCTGCCAGcaaacca                                                                          >  1:1203756/1‑150 (MQ=255)
                                                                                aTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACa                                                                       <  2:1631364/150‑1 (MQ=255)
                                                                                 tGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAg                                                                      <  1:222016/150‑1 (MQ=255)
                                                                                      gtgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGGGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACa                                                                 <  1:1665042/150‑1 (MQ=255)
                                                                                        ggAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATa                                                               >  2:882845/1‑150 (MQ=255)
                                                                                        ggAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCAGGCGGGCGTGATGATTTATTCAGCGTTTGTCGAACGTATTAGTTCCACAAGGCGAGAATCGGCCAGCCAACCAACAGCAACATa                                                               >  1:677565/1‑150 (MQ=255)
                                                                                            cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAgaga                                                           <  2:1414574/150‑1 (MQ=255)
                                                                                              catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                         <  1:337072/150‑1 (MQ=255)
                                                                                              catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                         <  2:87222/150‑1 (MQ=255)
                                                                                              catacataCACTGAATACTATCTAATAAATCATAAAAAAGTCCGCCCGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                         <  2:939971/150‑1 (MQ=255)
                                                                                               atacatacACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGt                                                        <  2:1776726/150‑1 (MQ=255)
                                                                                                       acTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcc                                                >  2:1579715/1‑150 (MQ=255)
                                                                                                          gAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCg                                             <  2:407750/150‑1 (MQ=255)
                                                                                                                  tctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgc                                     >  2:1853174/1‑150 (MQ=255)
                                                                                                                  tctatAAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgc                                     <  2:395455/150‑1 (MQ=255)
                                                                                                                   ctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAATAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                                    >  1:970614/1‑150 (MQ=255)
                                                                                                                    tatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCACCCAACAGCAACATAGACATGTAAATCACCCCGATGATTgcgc                                   >  2:1632584/1‑150 (MQ=255)
                                                                                                                    tatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGGTTGGCGAACGGATTAGTTCCACAGGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAAGCCCCCCGAAGATTgcgc                                   >  1:947704/1‑150 (MQ=255)
                                                                                                                          aTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGAc                             >  1:422116/1‑150 (MQ=255)
                                                                                                                            cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGc                           >  2:1035681/1‑150 (MQ=255)
                                                                                                                             aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCATCATAGAGATGTAAATCCCCCCGATGATTGCGCCAAGACGcc                          >  1:363308/1‑150 (MQ=255)
                                                                                                                             aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                          <  2:645840/150‑1 (MQ=255)
                                                                                                                             aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                          >  1:1069540/1‑150 (MQ=255)
                                                                                                                             aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                          <  2:1687684/150‑1 (MQ=255)
                                                                                                                             aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                          >  2:1600997/1‑150 (MQ=255)
                                                                                                                              aaaaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCa                         <  2:277947/149‑1 (MQ=255)
                                                                                                                                aaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGt                       <  1:1206448/150‑1 (MQ=255)
                                                                                                                                   aaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAAACGGCCAGCCAACCAACAGCAAAATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGCCGCCAGTAAt                    <  1:941559/150‑1 (MQ=255)
                                                                                                                                    cgCCCGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTAATCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGATTCGGCCAGCCACCCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGCCGACAGTAAta                   >  1:1450271/2‑149 (MQ=255)
                                                                                                                                      cccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt                 >  2:607520/1‑150 (MQ=255)
                                                                                                                                       ccGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCttt                >  1:825018/1‑150 (MQ=255)
                                                                                                                                        cGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtttt               >  1:627947/1‑150 (MQ=255)
                                                                                                                                         gCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTg              <  2:144683/150‑1 (MQ=255)
                                                                                                                                         gCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTg              <  2:1074900/150‑1 (MQ=255)
                                                                                                                                              gTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTc         <  2:1536470/150‑1 (MQ=255)
                                                                                                                                                  gCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACAt     <  2:873993/150‑1 (MQ=255)
                                                                                                                                                   cTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACATa    <  2:662588/150‑1 (MQ=255)
                                                                                                                                                     tACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACATAGc  >  1:1306206/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCACATAGC  >  NC_000913/644925‑645223

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