Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,557,200 T→C V137A (GTG→GCG)  b2442 → putative prophage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,557,2000TC99.5% 646.5 / ‑5.0 209V137A (GTG→GCG) b2442putative prophage integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (101/107);  total (101/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 223 total aligned reads displayed.

CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289
                                                                                          |                                                                                         
cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         >  1:1362513/1‑93 (MQ=255)
cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         >  1:384090/1‑93 (MQ=255)
 ccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAggcggc                                                                                        >  2:7140652/1‑93 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       >  2:4134599/1‑93 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       >  1:1187813/1‑93 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       >  1:3770994/1‑93 (MQ=255)
     aTTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGaa                                                                                    >  1:586345/1‑93 (MQ=255)
       tGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGa                                                                                  <  1:1632983/93‑1 (MQ=255)
        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTCAGAAGGCGGCGGAAGAc                                                                                 >  1:7258572/1‑93 (MQ=255)
        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGAc                                                                                 >  2:2708458/1‑93 (MQ=255)
        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGAc                                                                                 <  2:5128055/93‑1 (MQ=255)
         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACt                                                                                <  1:5853381/93‑1 (MQ=255)
         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                <  2:1483712/93‑1 (MQ=255)
         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                >  2:2364360/1‑93 (MQ=255)
          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               <  2:608167/93‑1 (MQ=255)
             tGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGtt                                                                            <  2:3433421/93‑1 (MQ=255)
             tGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGtt                                                                            >  2:4418447/1‑93 (MQ=255)
              ggAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGttt                                                                           <  2:6371822/93‑1 (MQ=255)
              ggAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGttt                                                                           <  1:2676503/93‑1 (MQ=255)
                  ttAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTa                                                                       >  1:4789032/1‑93 (MQ=255)
                  ttAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTa                                                                       <  1:986897/93‑1 (MQ=255)
                   tAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAc                                                                      <  1:6353292/93‑1 (MQ=255)
                   tAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAc                                                                      <  1:6178867/93‑1 (MQ=255)
                    aaCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAcc                                                                     >  1:2500352/1‑93 (MQ=255)
                     aCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCa                                                                    >  1:4141510/1‑93 (MQ=255)
                       ccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAgg                                                                  <  2:3858097/93‑1 (MQ=255)
                           tGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTg                                                              <  2:1433747/93‑1 (MQ=255)
                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             >  2:7603774/1‑93 (MQ=255)
                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             >  1:5911179/1‑93 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            <  1:3717245/93‑1 (MQ=255)
                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAg                                                           <  1:239067/93‑1 (MQ=255)
                               aCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGCTTGCTACCAGGGGTGGAgg                                                          >  1:4325559/1‑93 (MQ=255)
                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggt                                                        <  1:4628094/93‑1 (MQ=255)
                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGGTTGCTACCAGGGTTGGAGggt                                                        >  2:1806342/1‑93 (MQ=255)
                                    cccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtggt                                                     <  1:4777601/93‑1 (MQ=255)
                                     ccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTc                                                    <  2:2565618/93‑1 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCg                                                   >  1:204839/1‑93 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCg                                                   <  2:456331/93‑1 (MQ=255)
                                       gCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCgg                                                  >  1:3367356/1‑93 (MQ=255)
                                       gCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCgg                                                  >  2:4983585/1‑93 (MQ=255)
                                       gCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCgg                                                  >  2:4037802/1‑93 (MQ=255)
                                       gCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCgg                                                  >  2:1265649/1‑93 (MQ=255)
                                       gCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCgg                                                  <  2:1331196/93‑1 (MQ=255)
                                         aaGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGaa                                                >  2:4943034/1‑93 (MQ=255)
                                         aaGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGaa                                                <  2:7758252/93‑1 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  1:1085232/1‑93 (MQ=255)
                                            gAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATcc                                             <  1:3657023/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           <  1:6325016/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           <  2:3636869/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           >  1:3251342/1‑93 (MQ=255)
                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGc                                        >  2:3183607/1‑93 (MQ=255)
                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    >  1:7427784/1‑93 (MQ=255)
                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    >  1:5158607/1‑93 (MQ=255)
                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    <  1:4037808/93‑1 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGTAGGTTTTAAGAAGGCGGCGGAAGACTGGATTTCTTCCTGGGTTGTAGGGTGGTCGGAAGTCTTTGCAATAc                                   >  2:7037784/1‑93 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   >  2:5028791/1‑93 (MQ=255)
                                                       tGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACa                                  >  1:7043398/1‑93 (MQ=255)
                                                       tGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACa                                  >  1:2119117/1‑93 (MQ=255)
                                                        gATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAg                                 >  1:5393022/1‑93 (MQ=255)
                                                         aTGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAgg                                <  1:6298378/93‑1 (MQ=255)
                                                         aTGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAgg                                <  1:1870006/93‑1 (MQ=255)
                                                          tGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGt                               <  2:283171/93‑1 (MQ=255)
                                                          tGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGt                               >  2:3761665/1‑93 (MQ=255)
                                                             gcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTccg                            >  2:7430316/1‑91 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           <  1:5596889/93‑1 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGCGGTCGGAATCCTATGCAAAACAGGTCAGa                           >  1:2763377/1‑93 (MQ=255)
                                                               tCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGAt                          >  2:5524801/1‑93 (MQ=255)
                                                                cAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATc                         <  1:2300779/93‑1 (MQ=255)
                                                                cAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATc                         >  2:4480181/1‑93 (MQ=255)
                                                                cAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATc                         <  1:3707272/93‑1 (MQ=255)
                                                                    aTGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCa                     >  1:823758/1‑93 (MQ=255)
                                                                     tGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCAc                    >  1:5282846/1‑93 (MQ=255)
                                                                       gcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  >  1:3115065/1‑93 (MQ=255)
                                                                       gcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  <  1:1800796/93‑1 (MQ=255)
                                                                       gcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  <  1:264460/93‑1 (MQ=255)
                                                                        cggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAg                 >  2:3817842/1‑93 (MQ=255)
                                                                         ggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGa                <  2:3688266/93‑1 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  1:2260657/1‑93 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  1:1436360/1‑93 (MQ=255)
                                                                           cAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaa              >  2:4606449/1‑93 (MQ=255)
                                                                           cAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaa              >  2:948279/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  2:7674930/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            <  2:2455235/93‑1 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  2:4770481/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTAAGATCGGCACTAGaaaaa            >  2:7212265/1‑93 (MQ=255)
                                                                              gATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAg           >  2:3978149/1‑93 (MQ=255)
                                                                              gATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAg           >  1:4993135/1‑93 (MQ=255)
                                                                               gtttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGa          <  1:6823585/92‑1 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  1:6457179/93‑1 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  1:5466183/93‑1 (MQ=255)
                                                                                ttgaagaagTGGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACATGTCAGATCGGCACTAGAAAAAGAc         <  2:4716373/90‑1 (MQ=255)
                                                                                ttaaagaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  2:6928292/90‑1 (MQ=255)
                                                                                 ttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACg        >  2:593345/1‑93 (MQ=255)
                                                                                 ttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACg        <  1:3954092/93‑1 (MQ=255)
                                                                                 ttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACg        <  1:1161703/93‑1 (MQ=255)
                                                                                  tAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGt       <  1:241773/93‑1 (MQ=255)
                                                                                    agaagGCGGCGGAAGACTGGTCTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     <  1:7203386/93‑1 (MQ=255)
                                                                                     gaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTa    <  1:3455297/93‑1 (MQ=255)
                                                                                      aagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTAt   >  1:5239684/1‑93 (MQ=255)
                                                                                       agGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATc  <  1:6709226/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: