Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 5.9% intergenic (+34/‑164) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0510CG 5.9% 233.3 / 4.5 82intergenic (+34/‑164)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base G (2/3);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGC  >  NC_000913/3780904‑3781200
                                                                                                                                                   |                                                                                                                                                      
aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAg                                                                                                                                                      <  1:590716/150‑1 (MQ=255)
aCGCTGACTGGCGCGAAATCGATCAGCGAAAATACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                      <  2:399568/150‑1 (MQ=255)
    tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgag                                                                                                                                                  <  1:756257/150‑1 (MQ=255)
     gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                                 >  1:922828/1‑150 (MQ=255)
      ccTGGCGCGAAAACGATCAGCGAAATTACGCTAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGaga‑g                                                                                                                                               <  2:127687/149‑1 (MQ=255)
        tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCGGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑ggg                                                                                                                                             <  2:808086/150‑1 (MQ=255)
         ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGTTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCCTAGGAGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGA‑gggg                                                                                                                                            >  2:291062/1‑150 (MQ=255)
         ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑ggcg                                                                                                                                            >  2:22703/1‑148 (MQ=255)
           cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCTCTCCCGCAGGGAGA‑GGGGAc                                                                                                                                          >  2:1243213/1‑150 (MQ=255)
              gAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGA‑gggttaggg                                                                                                                                       <  1:1042834/150‑7 (MQ=255)
               aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGtt                                                                                                                                      >  1:1140216/1‑150 (MQ=255)
               aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGtt                                                                                                                                      >  2:102340/1‑150 (MQ=255)
                   cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGgc                                                                                                                                  >  2:1110024/1‑150 (MQ=255)
                   cGATCAGCGAAATTACGCAAGATTCGCGGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGgc                                                                                                                                  <  1:1127712/150‑1 (MQ=255)
                     aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGgcgc                                                                                                                                <  1:102340/150‑1 (MQ=255)
                      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATTGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCt                                                                                                                               <  1:956623/150‑1 (MQ=255)
                           gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACCCCCTAACCCTCTCCCTCAGGGCGA‑GGGGACCGTTCGGCGCtgtatg                                                                                                                          >  2:293336/1‑150 (MQ=255)
                              aTTACGCAAGATTCGCTGGTGCAGCGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTAc                                                                                                                       <  2:935531/150‑1 (MQ=255)
                                  cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGG‑GGGGACCGGTCGGCGCTGGATGTACTccc                                                                                                                   >  1:49036/1‑150 (MQ=255)
                                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGTGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCg                                                                                                                                    >  1:915669/1‑132 (MQ=255)
                                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGTGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCg                                                                                                                                    <  2:915669/132‑1 (MQ=255)
                                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCAGTATGTACTCCCt                                                                                                                  >  1:1252285/1‑150 (MQ=255)
                                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCACAGGGAGA‑GGGGACCGGTCGGCGCTGTATGTACTgccgc                                                                                                                 >  2:1651225/1‑148 (MQ=255)
                                      aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCAc                                                                                                               <  1:547811/150‑1 (MQ=255)
                                          tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                                                           <  1:1633867/150‑1 (MQ=255)
                                          tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                                                           <  2:1349098/150‑1 (MQ=255)
                                           cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                          >  2:41033/1‑150 (MQ=255)
                                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTCAGGGAGA‑GGGGACCGTTCGGCGCGGTATGTACTCCCGCACTCTGaaa                                                                                                        >  2:329874/1‑150 (MQ=255)
                                              tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGTCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                       <  1:105842/150‑1 (MQ=255)
                                              tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                       <  2:186891/150‑1 (MQ=255)
                                               ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                      <  2:514886/150‑1 (MQ=255)
                                                     ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACACCCTCACTCTGCAACGCCACCg                                                                                                >  2:317760/1‑150 (MQ=255)
                                                     ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCTCAGGGAGA‑GGGGACCGTTAGGCGCTGTATGTACTACCTCACTCTGAAACGACACCg                                                                                                >  2:744662/1‑150 (MQ=255)
                                                     ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCTCACCCGCAGGGGGA‑GGGGGCCGTTCGGGGCTGTATGTACTCCCTCCCTCTTAAACgcccccc                                                                                                >  2:1704495/1‑144 (MQ=255)
                                                             gTAAAGCGGTGCCTGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGGCGTGTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                        <  1:1705572/150‑1 (MQ=255)
                                                                aaGAGTTGCCTTCGGCACTGGCTACCATTGCTAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTTAAACGACACCTCACTCtttttt                                                                                     <  1:987900/150‑1 (MQ=255)
                                                                aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttttt                                                                                     <  2:1001118/150‑1 (MQ=255)
                                                                aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCcttttt                                                                                     <  1:1446444/150‑1 (MQ=255)
                                                                 agagTTGCCTGCGGATCTTGCTCCCATTGATAAAGTGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTCAGGGAGA‑GGGGACCGTTCGGAGCTGTATGGACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                    <  2:170564/150‑1 (MQ=255)
                                                                   agTTGCCTGCGGCACTTGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCATGGAGA‑GGGGACCTTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                  <  1:975397/150‑1 (MQ=255)
                                                                    gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                 <  1:1010026/150‑1 (MQ=255)
                                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                <  1:18094/150‑1 (MQ=255)
                                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                <  2:922828/150‑1 (MQ=255)
                                                                       gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGA‑GGGGACCGTTCGGCGCTGTAAGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCAccca                                                                              >  2:289878/1‑149 (MQ=255)
                                                                         cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                                            <  1:1294798/150‑1 (MQ=255)
                                                                          tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGACGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGTAGA‑CGGTACCTATCGGCGCTGTATGTACTCCCACACAAATAAACGCCACCGCACTCTTTATTTCTCCCTCGc                                                                           >  1:283824/1‑150 (MQ=255)
                                                                                    gCTCCCATGGCGAAAGTGAATGAGGCATAGTAGTTTGCGCCCCCTCAACCTTACCCTCTCCCTCAGCGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCTCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                 <  2:280889/150‑1 (MQ=255)
                                                                                    gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                 <  2:473541/150‑1 (MQ=255)
                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                <  1:249584/150‑1 (MQ=255)
                                                                                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCTCCCCCTAACCCACTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATTTACTACATCAATATGAAACGACACCGCCCACATTTTTTATACCACGCCCCTCCGGGgaga                                                              >  2:948254/1‑150 (MQ=255)
                                                                                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                              <  1:1136973/150‑1 (MQ=255)
                                                                                        ccATGGCGAACGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                             <  2:735171/150‑1 (MQ=255)
                                                                                        ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                             <  1:620065/150‑1 (MQ=255)
                                                                                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTACGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                           <  2:1081238/150‑1 (MQ=255)
                                                                                            ggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGG‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgggag                                                         <  1:1607845/150‑3 (MQ=255)
                                                                                             gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                        <  2:817015/150‑1 (MQ=255)
                                                                                              cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                       >  2:190091/1‑150 (MQ=255)
                                                                                              cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGGTCGGCGCGGGATGTACTCCCGCCCGCGGAAACGCCCCCGCCCGCGTTTTTTCGCCCGCGCCCCGCCGGGGGGGGGGGCgg                                                       >  1:1155639/1‑150 (MQ=255)
                                                                                              cGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                       <  1:1211960/150‑1 (MQ=255)
                                                                                               gAAAGCGAATGCGGCATAGACGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGA‑GGGGACCCTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                      <  2:507334/150‑1 (MQ=255)
                                                                                                 aaGGGAATGCGGCATAGGCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                    <  1:642393/150‑1 (MQ=255)
                                                                                                  agcgAATGCGGCATAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTCTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                   <  2:817707/147‑1 (MQ=255)
                                                                                                  acggAATGCCGCATAGTCGTTTGCCCCCCCACAACCAAAACCACAACCTCAGGGACA‑CGGGACCGTTAGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTATTTTTTTCTCCCACCCCCCGCCGGGGAGAGGGCCGGGGTg                                                   >  2:692980/3‑150 (MQ=255)
                                                                                                     gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                <  1:771661/150‑1 (MQ=255)
                                                                                                        ttcgtaaTAGTCGTTTGAGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                             <  2:1641994/144‑1 (MQ=255)
                                                                                                        tGCGGCATAGTCGTTTGCGCCCCCACCCCCAAACCCACGACCACAGGGAGA‑GGGGGCCGTTAGGCGCAGTATGTTCTCCCACACACGGAAACGACCCCCCCCACGTTTTTTATACCACCCCCCAACCGGGAGAGGGGCCGGGTGTggggg                                             >  2:596144/1‑149 (MQ=255)
                                                                                                           ggCATAGGCGTGAGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                          <  2:312372/150‑1 (MQ=255)
                                                                                                                    gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                 >  2:1128062/1‑150 (MQ=255)
                                                                                                                          gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGcgccc                           >  2:664110/1‑146 (MQ=255)
                                                                                                                          gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTAGCCCGCCCCCGGCGGGGGTGGGGGCGGGGGGGGGGGTAAAAGCCCCCCCTGcgccc                           >  2:13894/1‑146 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCTCCCTCAGGGAGACGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGTGGAAAAGGACGCACTGCTCACCAAt                       >  1:663554/1‑150 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGGTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCAGCTCCCGAATTTTCt                  >  1:1646719/1‑150 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                  >  1:1353249/1‑150 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTGTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCGCTGCTCACCAATTTTCt                  >  2:133635/1‑150 (MQ=255)
                                                                                                                                    ccTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCTCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAGCAATTTTCTg                 >  1:545359/1‑150 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTAACCAATTTTCTGCTa              >  2:118322/1‑150 (MQ=255)
                                                                                                                                         cccTCTCCCTGAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                          <  2:242303/88‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTGAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                          >  1:242303/1‑88 (MQ=255)
                                                                                                                                          ccTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc           >  2:1480454/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGCAACGCCCCCGCCGTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCCGGCTGAGGGGACAAGGCCGCCGTGCTCTCCTATTTTCTGCTATCCt          >  1:452807/1‑150 (MQ=255)
                                                                                                                                               cccTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc      >  1:156131/1‑150 (MQ=255)
                                                                                                                                                 cTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc    >  2:52972/1‑150 (MQ=255)
                                                                                                                                                   cAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc  >  2:1343818/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGA‑GGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGC  >  NC_000913/3780904‑3781200

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