Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,151,093 A→G C20R (TGC→CGC)  b2073 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,0930AG99.6% 748.9 / ‑4.8 230C20R (TGC→CGC) b2073hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (110/119);  minor base T (0/1);  total (110/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 234 total aligned reads displayed.

CATCACGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCAACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAGGTC  >  NC_000913/2151003‑2151182
                                                                                          |                                                                                         
catcaCGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGAc                                                                                         <  1:486303/93‑1 (MQ=255)
 atcaCGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACg                                                                                        >  2:5627692/1‑93 (MQ=255)
   caCGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTg                                                                                      >  2:645586/1‑93 (MQ=255)
    aCGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTgg                                                                                     >  2:727879/1‑93 (MQ=255)
     cGAAAGGTAACCACTCCGGCTTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAATCATCCAGCAGCAAAATGCAGGGGCGACGTGGt                                                                                    <  1:4665226/93‑1 (MQ=255)
       aaaGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTc                                                                                  <  2:1643744/93‑1 (MQ=255)
        aaGGTACCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAACCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTcc                                                                                 <  1:2375480/93‑1 (MQ=255)
        aaGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTcc                                                                                 <  1:875207/93‑1 (MQ=255)
        aaGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCTTCATAGAACCAGAAAAATCCAGCAGCAAAAAGCAGGGGCGACGTGGTTCCg                                                                                >  2:1438035/1‑93 (MQ=255)
         aGGTAACCAATCCGGCGTTAAGTTCCTTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCg                                                                                >  1:5254963/1‑93 (MQ=255)
          ggTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCgg                                                                               <  1:5461311/93‑1 (MQ=255)
            tAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCggct                                                                             <  1:3468090/93‑3 (MQ=255)
            tAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGt                                                                             <  1:2671998/93‑1 (MQ=255)
            tAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGt                                                                             <  2:1167947/93‑1 (MQ=255)
            tAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGt                                                                             <  1:3127868/93‑1 (MQ=255)
             aaCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGtt                                                                            >  2:5634944/1‑93 (MQ=255)
             aaCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGtt                                                                            <  2:6522817/93‑1 (MQ=255)
             aaCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGtt                                                                            >  1:3110033/1‑93 (MQ=255)
              aCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTa                                                                           >  1:1490704/1‑93 (MQ=255)
              aCCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTa                                                                           >  2:5496840/1‑93 (MQ=255)
               ccAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTAc                                                                          <  1:6406864/93‑1 (MQ=255)
               ccAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTAc                                                                          <  2:4750960/93‑1 (MQ=255)
                cAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACt                                                                         <  1:3221362/93‑1 (MQ=255)
                 aaTCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTg                                                                        >  1:6646950/1‑93 (MQ=255)
                 aaTCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTg                                                                        >  2:3850003/1‑93 (MQ=255)
                  aTCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTgg                                                                       >  1:899024/1‑93 (MQ=255)
                   tCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGc                                                                      <  2:3392721/93‑1 (MQ=255)
                    ccGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCa                                                                     >  1:3996218/1‑93 (MQ=255)
                    ccGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCa                                                                     >  2:492609/1‑93 (MQ=255)
                     cGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCaa                                                                    <  2:3247983/93‑1 (MQ=255)
                        cGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAAt                                                                 >  1:991469/1‑93 (MQ=255)
                         gTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATc                                                                <  1:1690238/93‑1 (MQ=255)
                          ttAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCg                                                               >  2:5272638/1‑93 (MQ=255)
                          ttAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGAAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCg                                                               <  2:3566045/93‑1 (MQ=255)
                           tAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGc                                                              <  2:3147661/93‑1 (MQ=255)
                            aaGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCt                                                             >  2:4470710/1‑93 (MQ=255)
                              gTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCtgg                                                           >  1:6250086/1‑93 (MQ=255)
                               ttCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCtggt                                                          <  2:5827047/93‑1 (MQ=255)
                                tCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCtggtg                                                         >  2:4847955/1‑93 (MQ=255)
                                 cATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCtggtgg                                                        >  1:6582879/1‑93 (MQ=255)
                                  aTTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGc                                                       <  1:2781011/93‑1 (MQ=255)
                                  aTTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGc                                                       >  1:5750634/1‑93 (MQ=255)
                                   ttAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCa                                                      <  1:5894260/93‑1 (MQ=255)
                                   ttAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCa                                                      >  2:1679823/1‑93 (MQ=255)
                                   ttAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCa                                                      >  1:4761631/1‑93 (MQ=255)
                                   ttAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCa                                                      >  2:4945712/1‑93 (MQ=255)
                                    tAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCaa                                                     >  2:2670895/1‑93 (MQ=255)
                                     aaTAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCaaa                                                    <  1:2536779/93‑1 (MQ=255)
                                     aaTAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCaaa                                                    <  1:2913460/93‑1 (MQ=255)
                                      aTAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCATGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAAc                                                   <  2:5376163/93‑1 (MQ=255)
                                      aTAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAAc                                                   <  1:4275735/93‑1 (MQ=255)
                                        aGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGt                                                 >  2:4681216/1‑93 (MQ=255)
                                            cTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGAtt                                             <  2:3501431/93‑1 (MQ=255)
                                            cTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGAtt                                             <  2:394529/93‑1 (MQ=255)
                                             tACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGAttt                                            <  2:5892282/93‑1 (MQ=255)
                                             tACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGAttt                                            <  1:448036/93‑1 (MQ=255)
                                             tACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGAttt                                            >  1:5558732/1‑93 (MQ=255)
                                               ccGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGt                                          >  1:6501607/1‑93 (MQ=255)
                                               ccGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGt                                          <  1:2169459/93‑1 (MQ=255)
                                                cGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGtt                                         <  1:3906243/93‑1 (MQ=255)
                                                 gTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTc                                        >  1:5063451/1‑93 (MQ=255)
                                                 gTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTc                                        <  2:3359415/93‑1 (MQ=255)
                                                  ttCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTct                                       >  1:4730669/1‑92 (MQ=255)
                                                  ttCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCg                                       <  2:1274618/93‑1 (MQ=255)
                                                   tCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGc                                      <  2:6271149/93‑1 (MQ=255)
                                                    cATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCt                                     <  2:5206866/93‑1 (MQ=255)
                                                     aTAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTc                                    >  1:2435210/1‑93 (MQ=255)
                                                       aGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCAt                                  <  2:363218/93‑1 (MQ=255)
                                                        gAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATa                                 >  1:3795290/1‑93 (MQ=255)
                                                        gAACCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATa                                 >  1:2934175/1‑93 (MQ=255)
                                                         aaCCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATaa                                >  2:5087290/1‑93 (MQ=255)
                                                          aCCAGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAg                               <  2:3135657/93‑1 (MQ=255)
                                                             aGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGaaa                            <  1:2093389/93‑1 (MQ=255)
                                                             aGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGaaa                            <  2:2322097/93‑1 (MQ=255)
                                                             aGAAACATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGCGATTTGTTCGCTCATAAGaaa                            >  1:3516619/1‑93 (MQ=255)
                                                                aaCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATcc                         <  2:6073886/93‑1 (MQ=255)
                                                                aaCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATcc                         >  2:1351645/1‑93 (MQ=255)
                                                                aaCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATcc                         >  1:4150463/1‑93 (MQ=255)
                                                                aaCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATcc                         >  2:1593394/1‑93 (MQ=255)
                                                                aaCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATcc                         <  1:3313150/93‑1 (MQ=255)
                                                                 aCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCt                        <  1:5501997/93‑1 (MQ=255)
                                                                 aCATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCt                        <  2:216791/93‑1 (MQ=255)
                                                                  cATCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCtt                       <  1:3606179/93‑1 (MQ=255)
                                                                    tCCAGCAGCAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTc                     >  2:1789175/1‑93 (MQ=255)
                                                                       agcagcAAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTa                  <  2:3544654/93‑1 (MQ=255)
                                                                        gcagcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTag                 >  2:1972439/1‑93 (MQ=255)
                                                                        gcagcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTag                 >  2:4137078/1‑93 (MQ=255)
                                                                         cagcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTaga                <  1:166510/93‑1 (MQ=255)
                                                                           gcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAgaga              <  1:5359252/93‑1 (MQ=255)
                                                                           gcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAgaga              <  1:1944368/93‑1 (MQ=255)
                                                                           gcaAAATGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAgaga              >  1:1227413/1‑93 (MQ=255)
                                                                            caAAATGCAGGGGCTACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGaa             <  1:1340961/93‑1 (MQ=255)
                                                                               aaTGCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTa          <  2:5053837/93‑1 (MQ=255)
                                                                                  gCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCg       <  2:544370/93‑1 (MQ=255)
                                                                                  gCAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCg       >  1:6317400/1‑93 (MQ=255)
                                                                                   cAGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGa      <  2:674104/93‑1 (MQ=255)
                                                                                    aGGGGCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAg     <  1:5187663/93‑1 (MQ=255)
                                                                                     ggggCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAgg    >  2:2864361/1‑93 (MQ=255)
                                                                                     ggggCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAgg    >  1:1068363/1‑93 (MQ=255)
                                                                                       ggCGACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAGGTc  <  1:743647/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATCACGAAAGGTAACCAATCCGGCGTTAAGTTCATTAATAGGTCTACCGTTCATAGAACCAGAAACATCCAGCAGCAAAATGCAGGGGCAACGTGGTTCCGGGTTACTGGCAAAATCGCTGGTGGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAGGTC  >  NC_000913/2151003‑2151182

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: