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.2% 754.3 / NA 238V137A (GTG→GCG) b2442putative prophage integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (126/110);  minor base A (2/0);  total (128/110)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289
                                                                                          |                                                                                         
 ccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAggcggc                                                                                        <  1:1506057/93‑1 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCCAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       >  2:1100144/1‑93 (MQ=255)
   ttATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGcggcgg                                                                                      >  1:3175234/1‑93 (MQ=255)
     aTTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGaa                                                                                    <  2:4020986/93‑1 (MQ=255)
      ttGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAg                                                                                   <  2:4102639/93‑1 (MQ=255)
       tGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCTGCGGAAGa                                                                                  <  2:1211776/93‑1 (MQ=255)
       tGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGa                                                                                  >  2:1994152/1‑93 (MQ=255)
       tGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGa                                                                                  <  2:5548947/93‑1 (MQ=255)
          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               <  1:87313/93‑1 (MQ=255)
           aaTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTgg                                                                              >  2:6003541/1‑93 (MQ=255)
           aaTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTgg                                                                              >  1:2525206/1‑93 (MQ=255)
            aTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGt                                                                             <  2:1828146/93‑1 (MQ=255)
            aTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGt                                                                             <  1:4938136/93‑1 (MQ=255)
             tGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGtt                                                                            >  1:4610109/1‑93 (MQ=255)
                 aTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCt                                                                        >  2:2392248/1‑93 (MQ=255)
                  ttAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTa                                                                       <  2:4860538/93‑1 (MQ=255)
                   tAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAc                                                                      <  1:3439727/93‑1 (MQ=255)
                     aCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCa                                                                    >  2:2022106/1‑93 (MQ=255)
                      cccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   <  1:5226388/93‑1 (MQ=255)
                        cAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAggg                                                                 >  2:5867305/1‑93 (MQ=255)
                        cAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAggg                                                                 >  1:38405/1‑93 (MQ=255)
                         aGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGCTTTAAGAAGGCGGCGGCAGACTGGTTTGCTACCAGGGt                                                                >  1:2989952/1‑93 (MQ=255)
                         aGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGt                                                                >  1:5428296/1‑93 (MQ=255)
                         aGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGt                                                                >  1:181992/1‑93 (MQ=255)
                          gTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGtt                                                               >  2:4583176/1‑93 (MQ=255)
                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTCTCAGGGCTgg                                                             >  2:3751194/1‑93 (MQ=255)
                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             <  2:2143566/93‑1 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            <  2:5312593/93‑1 (MQ=255)
                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAg                                                           <  1:1547555/93‑1 (MQ=255)
                               aCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAgg                                                          >  2:431332/1‑93 (MQ=255)
                               aCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAgg                                                          >  1:596116/1‑93 (MQ=255)
                               aCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAgg                                                          >  1:1313112/1‑93 (MQ=255)
                                cAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAggg                                                         <  1:2424027/93‑1 (MQ=255)
                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGATGACTGGTTTGCTACCAGGGTTGGAGggt                                                        <  2:5336481/93‑1 (MQ=255)
                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggt                                                        <  1:3618545/93‑1 (MQ=255)
                                   gCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtgg                                                      <  1:5130637/93‑1 (MQ=255)
                                     ccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTc                                                    >  1:3081842/1‑93 (MQ=255)
                                     ccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTc                                                    >  2:3807705/1‑93 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCg                                                   >  2:127793/1‑93 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCg                                                   >  2:5154810/1‑93 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCg                                                   >  1:3983160/1‑93 (MQ=255)
                                        cAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGa                                                 >  1:1017064/1‑93 (MQ=255)
                                         aaGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGaa                                                <  1:4568773/93‑1 (MQ=255)
                                         aaGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGaa                                                <  2:4644944/93‑1 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  2:4077319/1‑93 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  2:1630576/1‑93 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  2:1618307/1‑93 (MQ=255)
                                           ggAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCGGGGTTGGAGGGTGGTCGGAATc                                              >  1:2022187/1‑93 (MQ=255)
                                           ggAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATc                                              >  2:2614830/1‑93 (MQ=255)
                                           ggAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATc                                              >  1:2653143/1‑93 (MQ=255)
                                            gAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATcc                                             <  1:4880000/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           <  1:4233434/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           <  2:5457705/93‑1 (MQ=255)
                                               aaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTAt                                          <  2:2913888/93‑1 (MQ=255)
                                               aaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTAt                                          >  1:545145/1‑93 (MQ=255)
                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGc                                        >  1:1081781/1‑93 (MQ=255)
                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    >  2:5099579/1‑93 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   >  1:6081624/1‑93 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   <  1:5303980/93‑1 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   >  1:3783890/1‑93 (MQ=255)
                                                       tGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACa                                  <  2:4484895/93‑1 (MQ=255)
                                                        gATGCGCTCAACATGGCGGTAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAg                                 >  2:6133269/1‑93 (MQ=255)
                                                         aTGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGGCGGAATCCTATGCAAAACAgg                                >  2:4551410/1‑93 (MQ=255)
                                                          tGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGt                               <  1:1734730/93‑1 (MQ=255)
                                                            cgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACATGTCa                             <  1:2370182/93‑1 (MQ=255)
                                                            cgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCa                             <  1:662658/93‑1 (MQ=255)
                                                            cgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCa                             <  1:4011667/93‑1 (MQ=255)
                                                            cgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGCTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCa                             >  2:2321509/1‑93 (MQ=255)
                                                             gcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAg                            >  1:134023/1‑93 (MQ=255)
                                                             gcTCAACATGGCGGCAGGATTTAAGAAGAAGGCGGGAGACTGGTTTGCTACCAGGGTTGGGGGGTGGGCGGCATCCTATGCAAAACAGGGCAg                            >  2:5619458/1‑93 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           >  2:832404/1‑93 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           >  2:3605642/1‑93 (MQ=255)
                                                               tCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGAt                          <  2:38405/93‑1 (MQ=255)
                                                                cAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATc                         <  1:2540576/93‑1 (MQ=255)
                                                                 aaCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCg                        >  1:5560423/1‑93 (MQ=255)
                                                                  aCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCgg                       >  1:2292592/1‑93 (MQ=255)
                                                                  aCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCgg                       >  2:2473711/1‑93 (MQ=255)
                                                                   cATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTGTGCAAAACAGGTCAGATCGGc                      <  2:4008629/93‑1 (MQ=255)
                                                                   cATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGc                      <  2:5047343/93‑1 (MQ=255)
                                                                   cATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGc                      >  2:1454623/1‑93 (MQ=255)
                                                                      ggcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACt                   >  1:4366455/1‑93 (MQ=255)
                                                                        cggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAg                 >  2:1044267/1‑93 (MQ=255)
                                                                         ggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGa                <  2:334996/93‑1 (MQ=255)
                                                                         ggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGa                <  1:5723803/93‑1 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  1:5488541/1‑93 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               <  1:6059287/93‑1 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  2:1511893/1‑93 (MQ=255)
                                                                           cAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaa              >  2:3296251/1‑93 (MQ=255)
                                                                             ggATTTAAGACGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaacaa            >  1:2224554/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGCGGTCCGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  1:5823773/1‑93 (MQ=255)
                                                                               aTTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTGGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGa          <  1:4443098/93‑1 (MQ=255)
                                                                               aTTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGa          <  2:3384080/93‑1 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  1:274304/93‑1 (MQ=255)
                                                                                  tAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGt       >  2:1334566/1‑93 (MQ=255)
                                                                                   aagaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTc      >  1:3482855/1‑93 (MQ=255)
                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     <  2:4370748/93‑1 (MQ=255)
                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     <  2:5894421/93‑1 (MQ=255)
                                                                                      aagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTAt   >  2:4149036/1‑93 (MQ=255)
                                                                                       agGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATc  >  2:51184/1‑93 (MQ=255)
                                                                                       agGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATc  <  2:3241625/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289

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