Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,509,139 A→G Y40C (TAC→TGC) 
T19A (ACG→GCG) 
b2390 →
b2391 →
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,1390AG99.2% 824.7 / ‑2.4 237Y40C (TAC→TGC) 
T19A (ACG→GCG) 
b2390
b2391
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (117/118);  total (119/118)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-02

Only 100 of 239 total aligned reads displayed.

CCCTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTACGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGTCATG  >  NC_000913/2509049‑2509228
                                                                                          |                                                                                         
cccTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCg                                                                                         <  2:4386942/93‑1 (MQ=255)
cccTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCg                                                                                         >  2:5207769/1‑93 (MQ=255)
  cTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCTATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGt                                                                                       <  1:5241579/93‑1 (MQ=255)
  cTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGt                                                                                       >  2:1357196/1‑93 (MQ=255)
  cTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGt                                                                                       <  2:3315345/93‑1 (MQ=255)
  cTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGt                                                                                       >  2:6011201/1‑93 (MQ=255)
    gATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCTATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTaa                                                                                     <  2:4518262/93‑1 (MQ=255)
    gATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTaa                                                                                     <  2:3337707/93‑1 (MQ=255)
    gATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTaa                                                                                     >  1:3406749/1‑93 (MQ=255)
     aTGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAAt                                                                                    >  2:1135127/1‑93 (MQ=255)
        gCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTAc                                                                                 <  2:1287152/93‑1 (MQ=255)
        gCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTAc                                                                                 <  2:6394622/93‑1 (MQ=255)
          aTTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTg                                                                               >  1:2894729/1‑93 (MQ=255)
           tttACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTgg                                                                              <  1:5078151/93‑1 (MQ=255)
            ttACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTggg                                                                             <  2:5776002/93‑1 (MQ=255)
             tACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGa                                                                            <  1:375128/93‑1 (MQ=255)
             tACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGa                                                                            >  2:3427835/1‑93 (MQ=255)
              aCCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGAc                                                                           >  1:106066/1‑93 (MQ=255)
              aCCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGAc                                                                           <  2:354585/93‑1 (MQ=255)
               ccccGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACg                                                                          >  1:4608876/1‑93 (MQ=255)
                cccGCTTGCAGCACACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACgg                                                                         <  1:5689650/93‑1 (MQ=255)
                cccGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACgg                                                                         >  1:5324328/1‑93 (MQ=255)
                  cGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtg                                                                       >  2:506644/1‑93 (MQ=255)
                   gCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtgg                                                                      <  2:3017194/93‑1 (MQ=255)
                   gCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtgg                                                                      >  1:150162/1‑93 (MQ=255)
                    cTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtggt                                                                     >  2:2953808/1‑93 (MQ=255)
                    cTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtggt                                                                     >  2:4710101/1‑93 (MQ=255)
                    cTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtggt                                                                     <  2:755490/93‑1 (MQ=255)
                    cTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACggtggt                                                                     <  1:4462348/93‑1 (MQ=255)
                     ttGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTc                                                                    <  1:1770247/93‑1 (MQ=255)
                     ttGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTc                                                                    >  1:6000173/1‑93 (MQ=255)
                     ttGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTc                                                                    <  1:3612340/93‑1 (MQ=255)
                     ttGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTc                                                                    <  2:4858102/93‑1 (MQ=255)
                       gcagcaAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCAc                                                                  <  2:5888229/93‑1 (MQ=255)
                       gcagcaAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCAc                                                                  <  1:308420/93‑1 (MQ=255)
                        cagcaAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACt                                                                 <  2:1790286/93‑1 (MQ=255)
                          gcaAACGCAGGTGAAATCGCCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTgg                                                               <  2:3234239/93‑1 (MQ=255)
                           caAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgc                                                              <  2:2131621/93‑1 (MQ=255)
                           caAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgc                                                              <  1:1669352/93‑1 (MQ=255)
                            aaaCGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgcg                                                             >  1:2441495/1‑93 (MQ=255)
                             aaCGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgcgc                                                            >  2:4935585/1‑93 (MQ=255)
                             aaCGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgcgc                                                            >  2:4119821/1‑93 (MQ=255)
                             aaCGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGgcgc                                                            <  2:19869/93‑1 (MQ=255)
                                gCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGAc                                                         <  2:1192024/93‑1 (MQ=255)
                                gCAGGTGAAATCACCCTACTGCAATCAATCAAATTACAAATTGGCGACCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGAc                                                         <  1:4896139/93‑1 (MQ=255)
                                  aGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCg                                                       >  2:662430/1‑93 (MQ=255)
                                  aGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCg                                                       <  1:4639112/93‑1 (MQ=255)
                                  aGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCg                                                       <  1:2828262/93‑1 (MQ=255)
                                    gTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTg                                                     >  2:3926811/1‑93 (MQ=255)
                                      gAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGac                                                   >  1:107073/1‑93 (MQ=255)
                                      gAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGac                                                   <  1:4957144/93‑1 (MQ=255)
                                         aTCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGactac                                                <  1:4254959/93‑1 (MQ=255)
                                          tCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGactact                                               <  1:1027625/93‑1 (MQ=255)
                                           cACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTg                                              >  1:368529/1‑93 (MQ=255)
                                            aCCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTgg                                             <  1:5388731/93‑1 (MQ=255)
                                             cccTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGc                                            >  2:3192251/1‑93 (MQ=255)
                                              ccTACTGCCATCAATAAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCa                                           <  2:2386407/93‑1 (MQ=255)
                                                tACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATc                                         <  2:1508105/93‑1 (MQ=255)
                                                tACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATc                                         <  1:485117/93‑1 (MQ=255)
                                                tACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATc                                         >  2:2655505/1‑93 (MQ=255)
                                                 aCTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCg                                        <  1:2645802/93‑1 (MQ=255)
                                                 aCTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCg                                        <  1:1378950/93‑1 (MQ=255)
                                                 aCTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCg                                        >  2:2483288/1‑93 (MQ=255)
                                                   tGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCa                                      <  1:5110702/93‑1 (MQ=255)
                                                    gCCATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCaa                                     >  1:1592843/1‑93 (MQ=255)
                                                      cATCAATCAAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAAtt                                   >  1:5123806/1‑93 (MQ=255)
                                                       atcaatcaAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTa                                  >  2:5337/1‑93 (MQ=255)
                                                       atcaatcaAATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTa                                  >  1:2339301/1‑93 (MQ=255)
                                                        tcaatcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTAt                                 >  2:227778/1‑93 (MQ=255)
                                                          aatcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGa                               >  1:2199097/1‑93 (MQ=255)
                                                           atcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAg                              >  1:4084929/1‑93 (MQ=255)
                                                            tcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGt                             <  1:3460635/93‑1 (MQ=255)
                                                            tcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGt                             <  1:4523295/93‑1 (MQ=255)
                                                            tcaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGt                             <  1:639955/93‑1 (MQ=255)
                                                            acaaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGt                             <  2:5257408/92‑1 (MQ=255)
                                                             caaATTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTg                            >  2:5041735/1‑93 (MQ=255)
                                                              aaaTTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTgg                           <  2:3576860/93‑1 (MQ=255)
                                                              aaaTTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTgg                           >  2:3991568/1‑93 (MQ=255)
                                                               aaTTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGgc                          <  2:30947/93‑1 (MQ=255)
                                                               aaTTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGgc                          >  2:3825326/1‑93 (MQ=255)
                                                               aaTTACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGgc                          >  1:1990090/1‑93 (MQ=255)
                                                                  tACAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCa                       >  1:930086/1‑93 (MQ=255)
                                                                    cAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCaaa                     <  1:3186723/93‑1 (MQ=255)
                                                                    cAAATTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCaaa                     >  1:6183809/1‑93 (MQ=255)
                                                                     aaaTTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCaaaa                    >  1:6325675/1‑93 (MQ=255)
                                                                      aaTTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCaaaaa                   >  2:6253865/1‑93 (MQ=255)
                                                                      aaTTGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCaaaaa                   >  2:6409879/1‑93 (MQ=255)
                                                                        ttGGCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAAcc                 >  2:5586360/1‑93 (MQ=255)
                                                                          ggCGATCGCGATCATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGt               <  2:101091/93‑1 (MQ=255)
                                                                           gcgatcgcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGtt              >  1:5279129/1‑93 (MQ=255)
                                                                           gcgatcgcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGtt              <  1:3545317/93‑1 (MQ=255)
                                                                              atcgcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGTGTGGCGCAAAAACCGTtggt           >  1:4999553/1‑93 (MQ=255)
                                                                              atcgcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTtggt           <  2:1917007/93‑1 (MQ=255)
                                                                               tcgcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTtggtg          <  2:2759408/93‑1 (MQ=255)
                                                                                 gcgatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGc        >  1:3861075/1‑93 (MQ=255)
                                                                                   gatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGt      >  2:2204892/1‑93 (MQ=255)
                                                                                   gatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGt      <  1:5500495/93‑1 (MQ=255)
                                                                                   gatcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGt      >  2:2436202/1‑93 (MQ=255)
                                                                                     tcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGTCa    >  1:1233293/1‑93 (MQ=255)
                                                                                     tcATTGCGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGTCa    >  1:3362498/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCCTGATGGCATTTACCCCGCTTGCAGCAAACGCAGGTGAAATCACCCTACTGCCATCAATCAAATTACAAATTGGCGATCGCGATCATTACGGTAATTACTGGGACGGTGGTCACTGGCGCGACCGTGACTACTGGCATCGCAATTATGAGTGGCGCAAAAACCGTTGGTGGCGTCATG  >  NC_000913/2509049‑2509228

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