Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,996 T→A 9.2% Q671L (CAG→CTG)  torZ ← trimethylamine N‑oxide reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,9960TA 9.2% 270.8 / 11.5 93Q671L (CAG→CTG) torZtrimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base A (6/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCG  >  NC_000913/1954847‑1955143
                                                                                                                                                     |                                                                                                                                                   
gCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGTGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCt                                                                                                                                                     >  2:971875/1‑150 (MQ=255)
 cACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTg                                                                                                                                                    >  1:829528/1‑150 (MQ=255)
 cACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTCCCGGATGTGCGGTCAGAAGCTg                                                                                                                                                    >  1:1381455/1‑150 (MQ=255)
    cAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTaa                                                                                                                                                 >  1:1744778/1‑150 (MQ=255)
     aGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAAc                                                                                                                                                >  2:409477/1‑150 (MQ=255)
         cGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCt                                                                                                                                            <  1:190697/150‑1 (MQ=255)
            cATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTtct                                                                                                                                         <  2:505655/150‑1 (MQ=255)
            cATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTtct                                                                                                                                         <  1:1774767/150‑1 (MQ=255)
                cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCtcgtc                                                                                                                                     >  2:1757477/1‑150 (MQ=255)
                cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      >  2:1585273/1‑146 (MQ=255)
                cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      <  1:1585273/149‑4 (MQ=255)
                 ccAATACCCACCCGTGCAGCATCTTCGGCGTGAATAGTGATTGGTTCACTATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGGGCGGTCAGAAGCAGTTACAGCTtctgaccg                                                                                                                                    <  2:1027948/149‑6 (MQ=255)
                  cAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTcgg                                                                                                                                   >  2:1582908/1‑150 (MQ=255)
                   aaTACCAAAGCGTGCATCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTcggc                                                                                                                                  <  2:799882/150‑1 (MQ=255)
                   aaTACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTcggc                                                                                                                                  <  2:1402665/150‑1 (MQ=255)
                     tACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTcggcgg                                                                                                                                <  2:1824054/150‑1 (MQ=255)
                         aaaGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGTGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTAcc                                                                                                                            <  1:971875/150‑1 (MQ=255)
                         aaaGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCTTCGGCGGTAcc                                                                                                                            <  2:848974/150‑1 (MQ=255)
                           aGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCtt                                                                                                                          >  1:1865204/1‑150 (MQ=255)
                             cGTGCAGCATCTTCGGTGTGAATAGTGCTTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCTGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTc                                                                                                                        <  2:1503673/150‑1 (MQ=255)
                             cGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGTCTGTGTAAACGGTGTGGCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTc                                                                                                                        <  2:1443610/150‑1 (MQ=255)
                                gcagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCAt                                                                                                                     >  2:612599/1‑150 (MQ=255)
                                 cagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCAtt                                                                                                                    >  1:1000394/1‑150 (MQ=255)
                                 cagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCAtt                                                                                                                    >  1:623893/1‑150 (MQ=255)
                                   gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                  <  2:1556017/150‑1 (MQ=255)
                                   gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                  <  1:1140456/150‑1 (MQ=255)
                                   gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                  >  1:563861/1‑150 (MQ=255)
                                      tCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtcatca                                                                                                               <  1:1582908/150‑1 (MQ=255)
                                      tCTTCGGTGGGAATAGTGATTGGGTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtcatca                                                                                                               <  1:1798905/150‑1 (MQ=255)
                                       cTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAg                                                                                                              >  1:704794/1‑150 (MQ=255)
                                       cTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGAATATTTTTTACGCAGATCCGCATAGTTAAGCTGCCTGTCTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAAGTGCTTCTCGTCGGCGGTAAGAATCCATTcagcag                                                                                                              >  2:1321582/1‑150 (MQ=255)
                                        ttCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAgg                                                                                                             >  2:1907074/1‑150 (MQ=255)
                                         tCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGa                                                                                                            >  2:905540/1‑150 (MQ=255)
                                           ggTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGagc                                                                                                          >  2:188273/1‑150 (MQ=255)
                                             tgtgAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGagcca                                                                                                        >  2:804919/1‑150 (MQ=255)
                                               tgAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGagccagc                                                                                                      >  2:1812326/1‑150 (MQ=255)
                                                 aaTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCACTCATCAGGAgccagcca                                                                                                    <  1:767706/150‑1 (MQ=255)
                                                  aTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGccagccag                                                                                                   <  1:1587200/150‑1 (MQ=255)
                                                  aTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGccagccag                                                                                                   <  1:97139/150‑1 (MQ=255)
                                                       gATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTgg                                                                                              >  1:667070/1‑150 (MQ=255)
                                                        aTTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTggg                                                                                             <  2:829528/150‑1 (MQ=255)
                                                        aTTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTggg                                                                                             <  2:1006568/150‑1 (MQ=255)
                                                          tGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGgtg                                                                                           <  2:717497/150‑1 (MQ=255)
                                                          tGGTTCACGATCGGCGATCGCATATTTTTTACGCAGATCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGgtg                                                                                           >  2:849340/1‑150 (MQ=255)
                                                            gTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGtgtg                                                                                         <  1:152541/150‑1 (MQ=255)
                                                                 cGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCggg                                                                                    >  2:1762077/1‑150 (MQ=255)
                                                                  gATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGc                                                                                   <  1:615570/150‑1 (MQ=255)
                                                                   aTCGGCGCTCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCa                                                                                  >  2:1155759/1‑150 (MQ=255)
                                                                       gCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTCCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATcc                                                                              <  2:1381455/150‑1 (MQ=255)
                                                                          ttCGCATATTTTTTACGCAGTTCCGCATAGTTATGCTGACTGTGTATACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTta                                                                           <  2:1079318/149‑1 (MQ=255)
                                                                            cGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTtata                                                                         <  1:803448/150‑1 (MQ=255)
                                                                            cGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGGAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTAACTCCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTtata                                                                         >  2:843587/1‑150 (MQ=255)
                                                                               atatTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGGCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGcc                                                                      <  2:1022952/150‑1 (MQ=255)
                                                                                  ttttttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCaaa                                                                   <  1:804919/150‑1 (MQ=255)
                                                                                       tACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAAttttt                                                              <  1:1757477/150‑1 (MQ=255)
                                                                                        aCGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTc                                                             <  2:1186739/150‑1 (MQ=255)
                                                                                         cGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGTTGTGCCGGATGTGCGGTGAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTcc                                                            <  2:1210991/150‑1 (MQ=255)
                                                                                         cGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTcc                                                            >  1:1504189/1‑150 (MQ=255)
                                                                                         cGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCCGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTcc                                                            >  2:1183837/1‑150 (MQ=255)
                                                                                         cGCAGTTCCGCATAGTTAAGCTGACTGTGGAAACGGTGGGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTcc                                                            <  2:1275311/150‑1 (MQ=255)
                                                                                            aGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATGGCCAAATTTTTCCAAc                                                         <  1:905540/150‑1 (MQ=255)
                                                                                               tCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGtt                                                      >  1:28344/1‑150 (MQ=255)
                                                                                                  gCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGttttt                                                   >  1:1794095/1‑150 (MQ=255)
                                                                                                   cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTg                                                  <  2:623893/150‑1 (MQ=255)
                                                                                                   cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTg                                                  <  2:563861/150‑1 (MQ=255)
                                                                                                   cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      >  2:1146666/1‑63 (MQ=37)
                                                                                                   cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      <  1:1146666/66‑4 (MQ=37)
                                                                                                       gTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATa                                              >  1:1315340/1‑150 (MQ=255)
                                                                                                         tAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGa                                            >  2:1598133/1‑150 (MQ=255)
                                                                                                           aGCTGACTGTGTAAACGGTGTTCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGAtt                                          >  1:1795855/1‑150 (MQ=255)
                                                                                                           aGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGAtt                                          >  1:1367781/1‑150 (MQ=255)
                                                                                                            gCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGAttt                                         >  2:1583393/1‑150 (MQ=255)
                                                                                                              tGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCACGTTGGGTGTGCCGGGCAATCCTTATAGCAAAAATTTTCAAACGTTTTTGAATAGATTTCa                                       >  1:1263557/1‑150 (MQ=255)
                                                                                                                aCTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAt                                     <  2:621654/150‑1 (MQ=255)
                                                                                                                   gtgtAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                                  >  1:1908487/1‑150 (MQ=255)
                                                                                                                     gtAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGc                                <  2:1334905/150‑1 (MQ=255)
                                                                                                                     gtAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGc                                <  2:1794095/150‑1 (MQ=255)
                                                                                                                      tAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGcc                               >  1:1519633/1‑150 (MQ=255)
                                                                                                                       aaaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGGCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGGGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCg                              <  1:1321582/150‑1 (MQ=255)
                                                                                                                        aaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGc                             >  1:1246602/1‑150 (MQ=255)
                                                                                                                        aaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTGGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGc                             <  1:1812326/150‑1 (MQ=255)
                                                                                                                           ggTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGTCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTg                          <  1:1199690/150‑1 (MQ=255)
                                                                                                                           ggTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTg                          <  2:667070/150‑1 (MQ=255)
                                                                                                                            gtgtGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                                  >  2:281515/1‑141 (MQ=255)
                                                                                                                            gtgtGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                                  <  1:281515/141‑1 (MQ=255)
                                                                                                                             tgtgCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGc                        >  1:762587/1‑150 (MQ=255)
                                                                                                                              gtgCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCg                       >  2:403450/1‑150 (MQ=255)
                                                                                                                                 ccGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTAc                    <  1:1762077/150‑1 (MQ=255)
                                                                                                                                        tgCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCCTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAgcgc             >  1:1197231/1‑150 (MQ=255)
                                                                                                                                          cGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAt           <  2:1193197/150‑1 (MQ=255)
                                                                                                                                           gacgaGAAGCAGTTACAGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAtt          >  2:1535055/5‑150 (MQ=255)
                                                                                                                                            acgaGAAGCAGTTACAGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAttt         >  1:1614509/4‑150 (MQ=255)
                                                                                                                                             cgaGAAGCAGTTACAGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAtttt        >  1:239840/3‑150 (MQ=255)
                                                                                                                                                   gCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCg  >  2:426519/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCG  >  NC_000913/1954847‑1955143

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: