Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,955,002 T→A 9.0% Q669L (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,955,0020TA 9.0% 266.1 / 13.1 89Q669L (CAG→CTG) torZtrimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (37/44);  new base A (5/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCG  >  NC_000913/1954856‑1955151
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                              >  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:623893/1‑150 (MQ=255)
                        cagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCAtt                                                                                                                            >  1:1000394/1‑150 (MQ=255)
                          gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                          <  1:1140456/150‑1 (MQ=255)
                          gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                          >  1:563861/1‑150 (MQ=255)
                          gcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATtca                                                                                                                          <  2:1556017/150‑1 (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                                          <  1:281515/141‑1 (MQ=255)
                                                                                                                   gtgtGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                                          >  2:281515/1‑141 (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)
                                                                                                                                                  tGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGgcg  >  1:163230/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCG  >  NC_000913/1954856‑1955151

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: ‑