Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,038,858 G→C 5.0% P296A (CCG→GCG)  xerD ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,038,8580GC 5.0% 272.6 / 4.3 80P296A (CCG→GCG) xerDsite‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base C (2/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

CGCTGCTTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTG  >  NC_000913/3038711‑3039002
                                                                                                                                                   |                                                                                                                                                
cGCTGCTTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGt                                                                                                                                                >  2:585072/1‑150 (MQ=255)
    gctTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATg                                                                                                                                            >  1:1427520/1‑150 (MQ=255)
    gctTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATg                                                                                                                                            <  1:1414640/150‑1 (MQ=255)
    gctTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATg                                                                                                                                            <  2:113884/150‑1 (MQ=255)
     ctTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGc                                                                                                                                           >  1:897231/1‑150 (MQ=255)
     ctTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGc                                                                                                                                           >  2:230670/1‑150 (MQ=255)
       tttGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACCCGCGGGCGTGATgtct                                                                                                                                         <  2:277632/150‑4 (MQ=255)
           aTGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACCCGCGGGCGTGATgtctgaaa                                                                                                                                     <  1:1337717/150‑8 (MQ=255)
            tGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTtgat                                                                                                                                    <  2:27149/150‑1 (MQ=255)
             gCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTtgatg                                                                                                                                   >  2:1373427/1‑150 (MQ=255)
              cccATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTtgatga                                                                                                                                  <  1:1353728/150‑1 (MQ=255)
              cccATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTtgatga                                                                                                                                  >  2:1207568/1‑150 (MQ=255)
                cATTTTGGCTAACGTTTGTTGAATTGCCTCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAg                                                                                                                                >  2:1030151/1‑150 (MQ=255)
                 aTTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGt                                                                                                                               <  1:1430486/150‑1 (MQ=255)
                     tGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGcc                                                                                                                           <  1:883718/150‑1 (MQ=255)
                     tGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGcc                                                                                                                           <  1:1207568/150‑1 (MQ=255)
                       gCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCggg                                                                                                                                                 >  1:963340/1‑126 (MQ=255)
                       gCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCggg                                                                                                                                                 <  2:963340/126‑1 (MQ=255)
                           aCGTTTGTTGAATTGCCGCGGCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGAc                                                                                                                     <  2:942613/150‑1 (MQ=255)
                                   tgAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTa                                                                                                             >  1:1400053/1‑150 (MQ=255)
                                     aaTTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTATCAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGTTGAAGTTGCCGCAGACGCTCGGTAGc                                                                                                           <  1:934346/150‑1 (MQ=255)
                                      aTTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCg                                                                                                          <  1:133597/150‑1 (MQ=255)
                                      aTTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCg                                                                                                          >  1:870528/1‑150 (MQ=255)
                                      aTTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCg                                                                                                          >  2:1317251/1‑150 (MQ=255)
                                             cgTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCg                                                                                                   <  2:505204/150‑1 (MQ=255)
                                             cgTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCg                                                                                                   <  1:37033/150‑1 (MQ=255)
                                                 atcaGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGtata                                                                                               >  1:1376255/1‑150 (MQ=255)
                                                     gCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAAtt                                                                                           >  1:135213/1‑150 (MQ=255)
                                                      ccTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAAttt                                                                                          >  2:8881/1‑150 (MQ=255)
                                                       cTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAg                                                                                                                       >  1:1114703/1‑120 (MQ=255)
                                                       cTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAg                                                                                                                       <  2:1114703/120‑1 (MQ=255)
                                                        tGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGc                                                                                        <  1:641403/150‑1 (MQ=255)
                                                           gCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCgtg                                                                                     <  2:1400053/150‑1 (MQ=255)
                                                               aggcTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGtggtgg                                                                                 <  1:611769/147‑1 (MQ=255)
                                                               aGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGTTGATTATGTTTATGAAGTTGCCGCATAATCTCGTTTGTTAAATGCGTATAAATTTGAgttgtgt                                                                                 <  2:516334/150‑2 (MQ=255)
                                                               aGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGtggtgg                                                                                 >  2:1055040/1‑150 (MQ=255)
                                                                  cTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgaga                                                                              >  2:1099004/1‑150 (MQ=255)
                                                                  cTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCAGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgaga                                                                              >  1:529677/1‑150 (MQ=255)
                                                                   tGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgagag                                                                             >  2:303974/1‑150 (MQ=255)
                                                                   tGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCCTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgagag                                                                             >  2:723731/1‑150 (MQ=255)
                                                                    gAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgagaga                                                                            >  2:107945/1‑150 (MQ=255)
                                                                    gAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGgagaga                                                                            >  1:1325127/1‑150 (MQ=255)
                                                                        aCGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGc                                                                        <  1:104205/150‑1 (MQ=255)
                                                                        aCGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGc                                                                        <  1:1373427/150‑1 (MQ=255)
                                                                         cgccgcTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCt                                                                       >  1:96802/1‑150 (MQ=255)
                                                                            cgcTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCtgtg                                                                    <  1:8881/150‑1 (MQ=255)
                                                                             gcTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTgg                                                                   <  1:342133/150‑1 (MQ=255)
                                                                              cTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGc                                                                  <  1:200735/150‑1 (MQ=255)
                                                                              cTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGc                                                                  >  2:696921/1‑150 (MQ=255)
                                                                                aaCAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGccc                                                                >  1:1276819/1‑150 (MQ=255)
                                                                                aaCAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGccc                                                                <  1:585072/150‑1 (MQ=255)
                                                                                  cAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAg                                                              <  2:1048753/150‑1 (MQ=255)
                                                                                   aaaGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGt                                                             >  1:419145/1‑150 (MQ=255)
                                                                                      gTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGc                                                          >  2:828294/1‑150 (MQ=255)
                                                                                          acaacaTAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCt                                                      <  1:1317251/150‑1 (MQ=255)
                                                                                             acaTAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGca                                                   <  1:362582/150‑1 (MQ=255)
                                                                                               aTAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGcacc                                                 >  2:1165407/1‑150 (MQ=255)
                                                                                                tAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGccc                                                                >  1:641092/1‑134 (MQ=255)
                                                                                                tAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGccc                                                                <  2:641092/134‑1 (MQ=255)
                                                                                                  aaaCCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCAcg                                              >  2:1205327/1‑150 (MQ=255)
                                                                                                    aCCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCAcgcg                                            <  2:96802/150‑1 (MQ=255)
                                                                                                           ttcATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATcc                                     >  1:416820/1‑150 (MQ=255)
                                                                                                             cATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGc                                   <  1:107945/150‑1 (MQ=255)
                                                                                                               tAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCAc                                 >  1:346058/1‑150 (MQ=255)
                                                                                                                    cTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGa                            >  1:858567/1‑150 (MQ=255)
                                                                                                                    cTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGa                            <  2:419145/150‑1 (MQ=255)
                                                                                                                     ttCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGAt                           <  1:255860/150‑1 (MQ=255)
                                                                                                                       cccGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGAttt                         <  1:1099004/150‑1 (MQ=255)
                                                                                                                         cGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTtaa                       <  1:1205327/150‑1 (MQ=255)
                                                                                                                             cTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATaaa                   >  1:1197885/1‑150 (MQ=255)
                                                                                                                                 tCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAg               <  2:960794/150‑1 (MQ=255)
                                                                                                                                  acagcATCACCCGCGGGCGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGt              >  1:1349590/5‑150 (MQ=255)
                                                                                                                                    agcATCACCCGCGGGCGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTgg            >  2:502575/3‑150 (MQ=255)
                                                                                                                                       aTCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCaa         >  2:1088958/1‑150 (MQ=255)
                                                                                                                                       aTCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACACTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCaa         >  1:415147/1‑150 (MQ=255)
                                                                                                                                         cACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGc                                                                        <  1:1165454/85‑1 (MQ=255)
                                                                                                                                         cACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGc                                                                        >  2:1165454/1‑85 (MQ=255)
                                                                                                                                          aCGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGTACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAg      <  1:230670/150‑1 (MQ=255)
                                                                                                                                          aCGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAg      >  2:81483/1‑150 (MQ=255)
                                                                                                                                            gCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCg    >  1:1068361/1‑150 (MQ=255)
                                                                                                                                             cccGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGt   >  1:867445/1‑150 (MQ=255)
                                                                                                                                              ccGCGGGTTATGCTGTTGATGTATTTGCCGCAGATTCTCGGTAGCGACATGCGTATAAATTTGAGTTTTGGAGAGATCGCTGTGTCCCAGTATCATCTTCACCACGCGTAAATCCGCACGATTATTTATTAAATTAGTGGCAAATGCGTg  <  2:513842/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGCTGCTTTTGATGCCCATTTTGGCTAACGTTTGTTGAATTGCCGCGTCATCAGCCTGAGCAAAGCCTGAAAACGCCGCTAACAAAGTAAACAACATAAAACCTTTCTTCATAAATCTTCCCGTTCTTTTCAGACATCACGCCCGCGGGTGATGCTGTTGATGAAGTTGCCGCAGACGCTCGGTAGCGACATGCGTATAAATTTGCGTGGTGGAGAGATCGCTGTGGCCCAGTAGCATCTGCACCACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTG  >  NC_000913/3038711‑3039002

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: