Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,620,2190AG60.0% 50.9 / 93.2 85M1133V (ATG→GTG) rhsBRhsB
Reads supporting (aligned to +/- strand):  ref base A (18/16);  new base G (23/28);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCA  >  NC_000913/3620198‑3620311
                     |                                                                                            
gTGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcg                       >  2:3120469/1‑93 (MQ=14)
 tGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgc                      >  1:5673153/1‑93 (MQ=17)
   gAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtc                    >  1:1109495/1‑93 (MQ=21)
   gAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtc                    <  2:5937377/93‑1 (MQ=21)
   gAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtc                    >  2:5110020/1‑93 (MQ=21)
   gAGCAGATGAAAAACCAGGTGGACCCGGTGGACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtc                    >  1:4524567/1‑93 (MQ=255)
    aGCAGATGAAAAACCAGGTGGCCCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtca                   <  2:3595368/93‑1 (MQ=255)
    aGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtca                   <  1:4222171/93‑1 (MQ=25)
     gCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcat                  <  2:5002896/93‑1 (MQ=37)
      cAGATGAAAACCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:5451492/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:2379705/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:2708329/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  2:5917725/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:3054734/93‑1 (MQ=255)
       aGATGCAAACCCAGATGGCCCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  2:5345596/93‑1 (MQ=255)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  1:4926990/93‑1 (MQ=17)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  2:4922460/93‑1 (MQ=17)
        gATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:1201879/1‑93 (MQ=18)
        gATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:3059147/1‑93 (MQ=18)
        gATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:5122527/1‑93 (MQ=18)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:3101350/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:5032289/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:421572/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCACTGGCGCTCATCAg               >  2:5799897/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGACGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:2743037/1‑93 (MQ=255)
         aTGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              >  2:5315639/1‑93 (MQ=18)
         aTGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  1:6246646/93‑1 (MQ=21)
         aTGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  1:1668513/93‑1 (MQ=18)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  1:178066/93‑1 (MQ=255)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  2:2541993/93‑1 (MQ=255)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              >  1:3009568/1‑93 (MQ=255)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  1:4742338/93‑1 (MQ=255)
          tgcaaacCCAGATGGCCCCGTTGTACACGCCGGCGCGAAAAATCCCCCTGTACCGCTGCGCCCAGCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:5234025/90‑1 (MQ=11)
          tgcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  1:4602706/90‑1 (MQ=16)
          tgcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             >  1:3831833/4‑93 (MQ=18)
          tgcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:1303627/90‑1 (MQ=25)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:4357282/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:3783267/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  1:3689921/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  1:2075887/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             >  1:4719001/1‑93 (MQ=37)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  2:1320412/91‑1 (MQ=25)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  2:3498633/3‑93 (MQ=21)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  2:3482276/93‑1 (MQ=255)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  1:1502999/1‑93 (MQ=255)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  1:4472566/93‑1 (MQ=255)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  1:77626/1‑93 (MQ=255)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  1:3759446/1‑93 (MQ=255)
            caaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           >  2:269583/2‑93 (MQ=25)
            caaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           <  2:2346065/92‑1 (MQ=37)
            aaaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           <  2:5547963/93‑1 (MQ=255)
            aaaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           >  1:2846212/1‑93 (MQ=255)
            aaaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           >  2:3906666/1‑93 (MQ=255)
            aaaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           >  1:5004762/1‑93 (MQ=255)
             aaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACgg          <  1:5272939/93‑1 (MQ=255)
             aaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGGGCTCATCAGCACgg          >  2:4743249/1‑93 (MQ=37)
             aaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACgg          >  1:781616/1‑93 (MQ=37)
              aaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         <  1:1508268/93‑1 (MQ=255)
              aaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         >  2:5049895/1‑93 (MQ=37)
              aaaCCAGATGGACCCGGTGTACACGCCGGCGCAAAAAATCCACCTGTGCCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         <  2:2243938/93‑1 (MQ=255)
               aaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        <  1:3592756/93‑1 (MQ=255)
               aaCCAGGTGGACCCCGTGTGCACGCGGGCGCAAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        <  1:5108479/93‑1 (MQ=14)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        <  2:2081997/93‑1 (MQ=39)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        >  2:3308450/1‑93 (MQ=37)
                aCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  1:2788878/1‑93 (MQ=255)
                aCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  2:579572/1‑93 (MQ=255)
                aCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  1:4419551/1‑93 (MQ=37)
                aCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       <  2:4837139/93‑1 (MQ=38)
                aCCAGATGGACCCGGCGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  1:3578831/1‑93 (MQ=37)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAATCCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:5642385/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:3962793/93‑1 (MQ=37)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:4110457/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:458834/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTCTGCACGCCGGCGCTAAAATTCCACCTGTCCCACTGCGACCATCGCGTCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:3884824/93‑1 (MQ=11)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:518671/93‑1 (MQ=39)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:570883/93‑1 (MQ=38)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      >  1:6079716/1‑93 (MQ=37)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      >  1:4787002/1‑93 (MQ=37)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:5450286/93‑1 (MQ=39)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGCCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:4326637/93‑1 (MQ=37)
                  cAGGTGGGCCCGGTGTACCCGCCGGCGCGAAAAATCCGCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:2691581/93‑1 (MQ=11)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  2:5007161/1‑93 (MQ=37)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  1:1487353/93‑1 (MQ=255)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:5257562/93‑1 (MQ=255)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  1:3789920/1‑93 (MQ=37)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  1:1260989/93‑1 (MQ=255)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  1:3340662/93‑1 (MQ=255)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  1:2921762/1‑93 (MQ=37)
                  cAGATGGACCCGGTGTCCACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  1:2353514/93‑1 (MQ=255)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  1:4195996/1‑93 (MQ=37)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:5529619/93‑1 (MQ=38)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  1:2803263/1‑93 (MQ=37)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:3513317/93‑1 (MQ=37)
                   aGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGCCCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgggg    <  2:1004830/93‑1 (MQ=255)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   >  2:3108319/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   <  2:1189730/93‑1 (MQ=39)
                     aTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCa  >  2:2345126/1‑93 (MQ=37)
                     aTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCa  >  2:6117173/1‑93 (MQ=37)
                     |                                                                                            
GTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCA  >  NC_000913/3620198‑3620311

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: