Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,620,2190AG61.5% 61.2 / 100.2 96M1133V (ATG→GTG) rhsBRhsB
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base G (25/34);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 106 total aligned reads displayed.

GTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCA  >  NC_000913/3620198‑3620311
                     |                                                                                            
 tGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgc                      >  2:756610/1‑93 (MQ=17)
 tGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgc                      <  2:4295647/93‑1 (MQ=17)
 tGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgc                      <  2:1363828/93‑1 (MQ=17)
 tGGAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCCGCCCCTGGCGCt                     >  2:3557265/1‑93 (MQ=11)
 gggAGCAGATAAAAACCCAGGTGGCCCGGGTTCACACGCCGCCGCAAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgc                      <  2:4446367/92‑1 (MQ=11)
  ggAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCt                     >  1:2178669/1‑93 (MQ=18)
  ggAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCt                     <  1:2483896/93‑1 (MQ=33)
   gAGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtc                    >  1:1109062/1‑93 (MQ=21)
    aGCAGATGAAAACCCAGGGGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtca                   <  2:448555/93‑1 (MQ=11)
    aGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtca                   <  1:840966/93‑1 (MQ=255)
    aGCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtca                   <  2:4143217/93‑1 (MQ=25)
     gCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcat                  >  1:4482942/1‑93 (MQ=255)
     gCAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcat                  <  1:2999999/93‑1 (MQ=37)
     gCAGATGAAAAACCAGGTGGACCCGGTGTACACGACGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcat                  >  1:3362157/1‑93 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTCCACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:2743789/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:3154723/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 >  1:3662488/1‑93 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:3725687/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:11169/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 >  1:676018/1‑93 (MQ=37)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  2:552841/93‑1 (MQ=255)
      cAGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatc                 <  1:1079212/93‑1 (MQ=255)
       ataTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  1:4236245/91‑1 (MQ=255)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  2:2728058/93‑1 (MQ=17)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  2:3789910/93‑1 (MQ=17)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  1:3652556/93‑1 (MQ=17)
       aGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGACCTGCCGCTGGCGCtcatca                >  1:3393420/1‑93 (MQ=12)
       aGATGCAAAACAAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                >  1:3324298/1‑93 (MQ=9)
       aGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  2:1702373/93‑1 (MQ=255)
       aGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                >  2:2802763/1‑93 (MQ=255)
       aGATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCtcatca                <  1:3057162/93‑1 (MQ=255)
        gATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:2727987/1‑93 (MQ=18)
        gATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:146731/1‑93 (MQ=18)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:2491969/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  1:3500276/1‑93 (MQ=37)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               >  2:4145375/1‑93 (MQ=255)
        gATGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAg               <  1:4065239/93‑1 (MQ=255)
         aTGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              >  2:1189112/1‑93 (MQ=18)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  2:4186420/93‑1 (MQ=255)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              <  1:3277476/93‑1 (MQ=255)
         aTGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGc              >  1:2485046/1‑93 (MQ=255)
          tgcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  1:4456971/90‑1 (MQ=255)
          tgcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             >  1:899784/4‑93 (MQ=18)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:2512762/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             >  1:1730888/1‑93 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:2165563/93‑1 (MQ=37)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  2:3452484/93‑1 (MQ=255)
          tGAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCa             <  1:1332391/93‑1 (MQ=255)
           gcacaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAACCCGCCTGTACCACTGCGCCCATCGCGGCCTGCCGCTGGCTCTCATCAGCAc            <  2:1829287/89‑1 (MQ=11)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  2:4421416/91‑1 (MQ=25)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  1:2192994/91‑1 (MQ=25)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  2:2161794/3‑93 (MQ=21)
           gcaaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            >  2:4265977/3‑93 (MQ=21)
           gAAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  1:3101861/93‑1 (MQ=255)
           aaaaaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCAc            <  1:1396394/92‑1 (MQ=255)
            caaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTCCCACTGCGCCCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           <  1:636830/92‑1 (MQ=255)
            aaaaaCCAGGTGGACCCGTTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg           <  2:4277448/93‑1 (MQ=255)
             aaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACgg          >  1:3072360/1‑93 (MQ=37)
              aaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         >  2:3345110/1‑93 (MQ=37)
              aaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         >  1:2889552/1‑93 (MQ=37)
              aaaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa         <  1:115107/93‑1 (MQ=255)
              aaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAACACGGa         >  1:3085780/1‑93 (MQ=37)
               aaCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        <  2:1009061/93‑1 (MQ=255)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        >  1:2641603/1‑93 (MQ=37)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        >  2:4539965/1‑93 (MQ=37)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        >  2:2633600/1‑93 (MQ=37)
               aaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGaa        >  2:3987076/1‑93 (MQ=37)
                aCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       <  2:2312056/93‑1 (MQ=255)
                aCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       <  1:4079330/93‑1 (MQ=255)
                aCCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  2:2721715/1‑93 (MQ=37)
                aCCAGATGGACCGGGTGTACACGCCGCCGCAAAAAATCCCCCTGTACCACGGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       <  1:3311436/93‑1 (MQ=14)
                aCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       >  2:1344245/1‑93 (MQ=37)
                aCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg       <  1:1766118/93‑1 (MQ=39)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCGCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:1947022/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:3293473/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      >  1:4187299/1‑93 (MQ=37)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:3981085/93‑1 (MQ=255)
                 ccAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:2383853/93‑1 (MQ=255)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:3991587/93‑1 (MQ=37)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:4228266/93‑1 (MQ=39)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  1:3760188/93‑1 (MQ=39)
                 ccAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgg      <  2:1650791/93‑1 (MQ=37)
                  cAGTTGGACCGGGTGTACACGCGGGCGCGAAAAACCCACCTGTACCACGGCGGCCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:3139314/93‑1 (MQ=14)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  1:3144337/1‑93 (MQ=37)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  2:3163121/1‑93 (MQ=37)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:1045565/93‑1 (MQ=255)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  2:2301529/1‑93 (MQ=37)
                  cAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:4207704/93‑1 (MQ=255)
                  cAGATGGACCCGGTGTACACGCGGGCGCGAAAAATCCCCCTGTACCACTGCGACCACCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:2063224/93‑1 (MQ=255)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  1:3624563/93‑1 (MQ=39)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     <  2:597997/93‑1 (MQ=38)
                  cAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAggg     >  2:137964/1‑93 (MQ=37)
                   aGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgggg    >  2:4483826/1‑93 (MQ=37)
                   aGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAgggg    >  2:3193464/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   >  1:4072265/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   >  1:2877955/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   >  1:4074375/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   >  2:1094612/1‑93 (MQ=37)
                    gATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGc   <  1:2069361/93‑1 (MQ=39)
                     aTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCa  >  1:3124432/1‑93 (MQ=37)
                     |                                                                                            
GTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCA  >  NC_000913/3620198‑3620311

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