Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,069 G→A 6.4% intergenic (‑66/+48) 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,0690GA 6.4% 396.5 / 9.6 125intergenic (‑66/+48)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base G (50/67);  new base A (2/6);  total (52/73)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 126 total aligned reads displayed.

CCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTC  >  NC_000913/644922‑645216
                                                                                                                                                   |                                                                                                                                                   
 cTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                  <  2:232393/150‑1 (MQ=255)
 cTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                  >  1:147422/1‑150 (MQ=255)
 cTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                  >  1:689541/1‑150 (MQ=255)
 cTGCAACGATAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACGGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                  <  2:1415384/150‑1 (MQ=255)
   gCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCtt                                                                                                                                                >  2:1178549/1‑150 (MQ=255)
     aaCGCTAAGGCGTTGGCAGAAGAGAAGGGGAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAc                                                                                                                                              <  2:1661871/150‑1 (MQ=255)
      aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                             <  2:1547510/150‑1 (MQ=255)
      aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                             >  1:1587585/1‑150 (MQ=255)
      aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                             <  1:651168/150‑1 (MQ=255)
       cGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                            >  1:1635319/1‑150 (MQ=255)
        gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGCCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTg                                                                                                                                           >  2:911169/1‑150 (MQ=255)
         cTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTgg                                                                                                                                          >  1:1394495/1‑150 (MQ=255)
            aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                       <  1:132403/150‑1 (MQ=255)
             cgCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGCGCCACGCGGCGGTACGCCAGACTGCTTGCATACTGGGTGGGGAACGCATACCTCCACTGAATACTATCTATTAACTCCTAAAACAGCCCGCCAGGTGGGCTTACCTggcggg                                                                                                                                      <  1:1169190/149‑1 (MQ=255)
                gTTGGCAGACGACAAGGGAAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATCCTTACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                   <  1:1473755/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:895070/1‑150 (MQ=255)
                      agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                             <  2:996977/150‑1 (MQ=255)
                      agaagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATCCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCAAAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                             <  1:125491/150‑1 (MQ=255)
                           agaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAtt                                                                                                                        <  1:1858887/150‑1 (MQ=255)
                            gaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCACAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                       <  1:1696773/150‑1 (MQ=255)
                             aagGGCAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGACTGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCa                                                                                                                      <  1:1699164/150‑1 (MQ=255)
                              agGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATTCTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAg                                                                                                                     <  2:659107/150‑1 (MQ=255)
                               gggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGCCTTACCTGGCGGGCGTGATGATTTATTCacc                                                                                                                    >  2:728571/1‑148 (MQ=255)
                                ggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCg                                                                                                                   <  2:1587585/150‑1 (MQ=255)
                                    gCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTg                                                                                                               >  1:1159541/1‑150 (MQ=255)
                                     cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                              <  1:504724/150‑1 (MQ=255)
                                      agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                             <  2:456384/150‑1 (MQ=255)
                                       gagaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAATAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTAATCGTTTGGCg                                                                                                            >  2:1819011/1‑150 (MQ=255)
                                         gaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGaa                                                                                                          >  2:396687/1‑150 (MQ=255)
                                           aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                        <  2:1636312/150‑1 (MQ=255)
                                           aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                        <  2:1320850/150‑1 (MQ=255)
                                           aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACg                                                                                                        <  1:1081825/150‑1 (MQ=255)
                                               gcgAGCAACGGGGGGTTACGCCAGAATGCTTTGATAATGGGGGTGGATCTCATACATACACTGAATACTATCTATTAACTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                    <  2:1618672/150‑1 (MQ=255)
                                               gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATACTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                    <  2:1532380/150‑1 (MQ=255)
                                               gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAtt                                                                                                    <  2:1613204/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)
                                                  aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                 >  2:191788/1‑150 (MQ=255)
                                                  aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                 >  2:1123897/1‑150 (MQ=255)
                                                  aGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGt                                                                                                 >  2:1037392/1‑150 (MQ=255)
                                                    cAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACCTAATAGtta                                                                                               >  2:1551236/1‑149 (MQ=255)
                                                       cgcgGCGTTACGCCAGAATGCTTTCATAATGGGTGTAGAACTCATACATACACTGAATACTATCTTTTATATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCTGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcac                                                                                            >  2:1907126/1‑150 (MQ=255)
                                                        gcggcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcaca                                                                                           <  2:970404/150‑1 (MQ=255)
                                                           gcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATgg                                                                                        <  1:52073/150‑1 (MQ=255)
                                                           gcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATgg                                                                                        <  1:1043162/150‑1 (MQ=255)
                                                              ttACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCga                                                                                     <  2:186634/150‑1 (MQ=255)
                                                               tACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCgag                                                                                    >  1:1652064/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)
                                                                       aaTGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgcc                                                                            >  2:581009/1‑150 (MQ=255)
                                                                         tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                          <  1:1713824/150‑1 (MQ=255)
                                                                         tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                          <  2:842273/150‑1 (MQ=255)
                                                                         tGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                          <  1:584566/150‑1 (MQ=255)
                                                                           cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                        >  1:1750469/1‑150 (MQ=255)
                                                                           cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                        >  1:785859/1‑150 (MQ=255)
                                                                           cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                        <  1:370995/150‑1 (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)
                                                                                          tgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACAt                                                         <  2:1635319/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                                                  <  2:87222/150‑1 (MQ=255)
                                                                                                 catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                  <  1:337072/150‑1 (MQ=255)
                                                                                                 catacataCACTGAATACTATCTAATAAATCATAAAAAAGTCCGCCCGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                  <  2:939971/150‑1 (MQ=255)
                                                                                                  atacatacACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGt                                                 <  2:1776726/150‑1 (MQ=255)
                                                                                                   tacatacaCTGATTACTTTCTATTAATTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTa                                                <  1:360862/150‑1 (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)
                                                                                                                        atTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGcc                           <  1:9158/150‑1 (MQ=255)
                                                                                                                             aTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGAc                      >  1:422116/1‑150 (MQ=255)
                                                                                                                               cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGc                    >  2:1714156/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:1563257/150‑1 (MQ=255)
                                                                                                                                           cGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtttt        >  1:627947/1‑150 (MQ=255)
                                                                                                                                            gCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTg       <  2:1139273/150‑1 (MQ=255)
                                                                                                                                             ccAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGa      >  2:1375600/1‑150 (MQ=255)
                                                                                                                                                 gTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTc  <  2:1536470/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTC  >  NC_000913/644922‑645216

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: