Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,672,529 T→C V198A (GTG→GCG)  b1601 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,5290TC99.5% 628.5 / ‑4.4 182V198A (GTG→GCG) b1601putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (97/84);  total (97/85)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

Only 100 of 191 total aligned reads displayed.

TGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAT  >  NC_000913/1672439‑1672618
                                                                                          |                                                                                         
tGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGgcgc                                                                                         >  1:4835132/1‑93 (MQ=255)
 gCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCt                                                                                        <  2:5667368/93‑1 (MQ=255)
 gCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCt                                                                                        >  2:2267653/1‑93 (MQ=255)
 gCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCt                                                                                        >  2:6214811/1‑93 (MQ=255)
     ggAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGaaa                                                                                    >  2:3391534/1‑93 (MQ=255)
     ggAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGaaa                                                                                    >  2:6203102/1‑93 (MQ=255)
       aaaaTTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGaaaaa                                                                                  >  1:5443594/1‑92 (MQ=255)
        aaaTTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACa                                                                                 >  2:5839478/1‑93 (MQ=255)
          aTTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGc                                                                               >  1:1442693/1‑93 (MQ=255)
          aTTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGc                                                                               >  2:341387/1‑93 (MQ=255)
          aTTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGc                                                                               <  2:5796574/93‑1 (MQ=255)
             tCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCAt                                                                            <  1:3081476/93‑1 (MQ=255)
             tCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCAt                                                                            >  2:2814362/1‑93 (MQ=255)
              cAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATc                                                                           >  1:2993698/1‑93 (MQ=255)
               aGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCa                                                                          <  2:4871614/93‑1 (MQ=255)
                gCAAATGATGGCGCGTCGGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAg                                                                         <  2:1748270/93‑1 (MQ=255)
                gCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAg                                                                         <  2:5597657/93‑1 (MQ=255)
                gCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAg                                                                         <  2:3827977/93‑1 (MQ=255)
                 cAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGc                                                                        <  2:3334680/93‑1 (MQ=255)
                    atgatgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGcatc                                                                     >  2:4792548/1‑93 (MQ=255)
                    atgatgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGcatc                                                                     >  1:5545459/1‑93 (MQ=255)
                     tgatgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGcatca                                                                    <  2:3822490/93‑1 (MQ=255)
                      gatgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGcatcat                                                                   <  2:4410018/93‑1 (MQ=255)
                      gatgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGcatcat                                                                   <  2:2991124/93‑1 (MQ=255)
                       atgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCatcatc                                                                  <  2:677207/93‑1 (MQ=255)
                        tgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCAtcatca                                                                 >  2:5075977/1‑93 (MQ=255)
                        tgGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCACATTATCTGGCGCTGAAAACAGCCATCAGCAtcatca                                                                 <  2:211201/93‑1 (MQ=255)
                         ggCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCAc                                                                >  2:3990042/1‑93 (MQ=255)
                          gcgcGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCAcc                                                               <  1:3500030/93‑1 (MQ=255)
                           cgcgTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCg                                                              <  2:4641523/93‑1 (MQ=255)
                            gcgTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCgg                                                             <  2:421302/93‑1 (MQ=255)
                             cgTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGc                                                            >  2:6265398/1‑93 (MQ=255)
                             cgTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGc                                                            <  1:5171486/93‑1 (MQ=255)
                             cgTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGc                                                            <  1:3092108/93‑1 (MQ=255)
                              ttCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGcc                                                           <  2:5283001/92‑1 (MQ=255)
                              gTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGcc                                                           <  1:887171/93‑1 (MQ=255)
                               tCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCt                                                          <  1:4895724/93‑1 (MQ=255)
                               tCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCt                                                          >  1:221936/1‑93 (MQ=255)
                                 cGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTgg                                                        >  2:817726/1‑93 (MQ=255)
                                   gTTGAAGGGATGGCGGCGATTCAACGGGGGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTc                                                      <  1:5333962/93‑1 (MQ=255)
                                    ttGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCg                                                     >  1:4783175/1‑93 (MQ=255)
                                     tGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGc                                                    >  1:3775318/1‑93 (MQ=255)
                                      gAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGcc                                                   >  2:2936400/1‑93 (MQ=255)
                                      gAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGcc                                                   <  1:322397/93‑1 (MQ=255)
                                      gAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGcc                                                   <  2:2392193/93‑1 (MQ=255)
                                       aaGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCt                                                  >  1:2710802/1‑93 (MQ=255)
                                        aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTg                                                 <  1:4794309/93‑1 (MQ=255)
                                        aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTg                                                 <  1:6031571/93‑1 (MQ=255)
                                        aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTg                                                 <  1:6071805/93‑1 (MQ=255)
                                        aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTg                                                 <  2:5943172/93‑1 (MQ=255)
                                          ggATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTggg                                               <  2:2082371/93‑1 (MQ=255)
                                            aTGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCg                                             >  1:4318135/1‑93 (MQ=255)
                                               gcggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATg                                          >  2:1614416/1‑93 (MQ=255)
                                               gcggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATg                                          >  2:2021916/1‑93 (MQ=255)
                                                cggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGc                                         >  1:1095458/1‑93 (MQ=255)
                                                 ggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCt                                        <  1:3735217/93‑1 (MQ=255)
                                                 ggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCt                                        <  1:599514/93‑1 (MQ=255)
                                                   cgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTcg                                      >  2:3799851/1‑93 (MQ=255)
                                                     aTTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTcgcc                                    >  1:408116/1‑93 (MQ=255)
                                                     aTTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTcgcc                                    <  1:1168704/93‑1 (MQ=255)
                                                      ttCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTcgccg                                   >  1:3951092/1‑93 (MQ=255)
                                                      ttCAACGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTcgccg                                   >  2:3905993/1‑93 (MQ=255)
                                                          aCGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACt                               >  2:1544142/1‑93 (MQ=255)
                                                          aCGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACt                               <  2:213890/93‑1 (MQ=255)
                                                          aCGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACt                               <  1:1579709/93‑1 (MQ=255)
                                                          aCGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACt                               >  1:6206873/1‑93 (MQ=255)
                                                          aCGTGCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACt                               >  1:5925561/1‑93 (MQ=255)
                                                              gCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGAt                           <  2:2275615/93‑1 (MQ=255)
                                                              gCGATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGAt                           <  1:200869/93‑1 (MQ=255)
                                                                gATTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGt                         <  1:4591725/93‑1 (MQ=255)
                                                                 aTTGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGtt                        <  2:3677811/93‑1 (MQ=255)
                                                                  ttGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                       <  2:709620/93‑1 (MQ=255)
                                                                  ttGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                       <  1:10622/93‑1 (MQ=255)
                                                                  ttGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                       >  2:1130051/1‑93 (MQ=255)
                                                                  ttGACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                       >  1:1941290/1‑93 (MQ=255)
                                                                    gACAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgc                     >  2:5200886/1‑93 (MQ=255)
                                                                      cAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgctt                   <  1:5317000/93‑1 (MQ=255)
                                                                      cAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgctt                   <  2:5905606/93‑1 (MQ=255)
                                                                      cAGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgctt                   >  1:3391661/1‑93 (MQ=255)
                                                                       aGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                  <  1:3584467/93‑1 (MQ=255)
                                                                       aGTGTTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                  <  1:2311503/93‑1 (MQ=255)
                                                                        gtgtTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcg                 >  1:5404911/1‑93 (MQ=255)
                                                                        gtgtTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcg                 >  1:1286661/1‑93 (MQ=255)
                                                                          gtTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTtcgcttcgct               >  2:4295092/1‑93 (MQ=255)
                                                                          gtTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTtcgcttcgct               >  2:1704825/1‑93 (MQ=255)
                                                                          gtTTCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTtcgcttcgct               >  1:622242/1‑93 (MQ=255)
                                                                           tttCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTcgcttcgctt              <  2:3569823/93‑1 (MQ=255)
                                                                            ttCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCttt             >  2:2446896/1‑93 (MQ=255)
                                                                            ttCTCATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCttt             <  2:4955157/93‑1 (MQ=255)
                                                                             tctcATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt            <  1:69077/93‑1 (MQ=255)
                                                                                cATTATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTc         >  1:3290938/1‑93 (MQ=255)
                                                                                  ttATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg       >  2:873177/1‑93 (MQ=255)
                                                                                  ttATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg       <  1:965788/93‑1 (MQ=255)
                                                                                  ttATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg       <  1:283710/93‑1 (MQ=255)
                                                                                  ttATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg       <  1:6145601/93‑1 (MQ=255)
                                                                                   tATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgg      <  1:1364037/93‑1 (MQ=255)
                                                                                   tATCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgg      >  1:595996/1‑93 (MQ=255)
                                                                                    aTCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg     <  1:2884914/93‑1 (MQ=255)
                                                                                    aTCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg     >  2:4485126/1‑93 (MQ=255)
                                                                                     tCTGGCGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgggg    >  2:5487970/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAT  >  NC_000913/1672439‑1672618

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: