Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,999 A→T 8.9% L670* (TTA→TAA)  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,9990AT 8.9% 267.3 / 13.4 90L670* (TTA→TAA) torZtrimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base T (5/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

CAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCG  >  NC_000913/1954851‑1955143
                                                                                                                                                    |                                                                                                                                                
cAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTaa                                                                                                                                                 >  1:1744778/1‑150 (MQ=255)
 aGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAAc                                                                                                                                                >  2:409477/1‑150 (MQ=255)
     cGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCt                                                                                                                                            <  1:190697/150‑1 (MQ=255)
        cATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTtct                                                                                                                                         <  1:1774767/150‑1 (MQ=255)
        cATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTtct                                                                                                                                         <  2:505655/150‑1 (MQ=255)
            cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCtcgtc                                                                                                                                     >  2:1757477/1‑150 (MQ=255)
            cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      <  1:1585273/149‑4 (MQ=255)
            cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCAGTTACAGCTtctgac                                                                                                                                      >  2:1585273/1‑146 (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:623893/1‑150 (MQ=255)
                             cagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCAtt                                                                                                                    >  1:1000394/1‑150 (MQ=255)
                               gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                  <  1:1140456/150‑1 (MQ=255)
                               gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                  <  2:1556017/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:563861/150‑1 (MQ=255)
                                                                                               cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTg                                                  <  2:623893/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)
                                                                                                                                                    |                                                                                                                                                
CAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCG  >  NC_000913/1954851‑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: ‑