Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,672,562 T→C L209P (CTG→CCG)  b1601 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,5620TC99.0% 313.3 / ‑4.7 103L209P (CTG→CCG) b1601putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (55/47);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 108 total aligned reads displayed.

GGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAATA  >  NC_000913/1672473‑1672651
                                                                                         |                                                                                         
ggTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGt                                                                                        <  2:785908/93‑1 (MQ=255)
  ttGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCg                                                                                      <  1:4013307/93‑1 (MQ=255)
   tGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGc                                                                                     >  1:2823876/1‑93 (MQ=255)
     aaGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCt                                                                                   >  2:633987/1‑93 (MQ=255)
      aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  >  1:3844159/1‑93 (MQ=255)
      aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  >  2:4234/1‑93 (MQ=255)
      aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  >  1:1845468/1‑93 (MQ=255)
       gggATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTgg                                                                                 >  1:3980201/1‑93 (MQ=255)
       gggATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTgg                                                                                 <  2:4550609/93‑1 (MQ=255)
       gggATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTgg                                                                                 <  2:2014539/93‑1 (MQ=255)
          aTGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCg                                                                              <  1:2087112/93‑1 (MQ=255)
           tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCCTTATCTGGTGCTGAAAACAGCCATCAGCATCATTACCGGCCCGGTCGCCTGGGCGa                                                                             >  2:1815295/1‑93 (MQ=255)
           tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             <  2:2573530/93‑1 (MQ=255)
           tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             <  1:2991328/93‑1 (MQ=255)
            ggcggcGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGAt                                                                            <  1:2346867/93‑1 (MQ=255)
             gcggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATg                                                                           >  1:506283/1‑93 (MQ=255)
               ggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCt                                                                         <  2:4205333/93‑1 (MQ=255)
                  gATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgc                                                                      >  1:1321792/1‑93 (MQ=255)
                  gATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgc                                                                      >  1:2860968/1‑93 (MQ=255)
                    ttCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgccg                                                                    >  2:3743993/1‑93 (MQ=255)
                    ttCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgccg                                                                    >  1:2383660/1‑93 (MQ=255)
                     tCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgccgc                                                                   <  2:3488814/93‑1 (MQ=255)
                     tCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTCCGCAACACAGCCATCAGCAGCATCACCGGCCCGGTCGCCTGGGCGATGCTcgccgc                                                                   <  2:4365434/93‑1 (MQ=255)
                        aCGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACt                                                                >  1:852747/1‑93 (MQ=255)
                           tGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGa                                                             >  1:394654/1‑93 (MQ=255)
                              gATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGt                                                          >  2:3452647/1‑93 (MQ=255)
                              gATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGt                                                          <  1:329219/93‑1 (MQ=255)
                                 tGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcg                                                       >  1:3109541/1‑93 (MQ=255)
                                  gACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgc                                                      >  1:2433401/1‑93 (MQ=255)
                                   aCAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgct                                                     >  1:2907692/1‑93 (MQ=255)
                                    cAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgctt                                                    <  2:1112384/93‑1 (MQ=255)
                                     aGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                                                   <  2:468690/93‑1 (MQ=255)
                                     aGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                                                   >  1:47081/1‑93 (MQ=255)
                                     aGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                                                   <  1:3248887/93‑1 (MQ=255)
                                       tgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcgc                                                 >  2:2421169/1‑93 (MQ=255)
                                       tgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcgc                                                 <  2:1561406/93‑1 (MQ=255)
                                       tgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcgc                                                 >  2:848754/1‑93 (MQ=255)
                                       tgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcgc                                                 <  1:313507/93‑1 (MQ=255)
                                        gtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTtcgcttcgct                                                >  1:3887003/1‑93 (MQ=255)
                                        gtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTtcgcttcgct                                                >  1:2905612/1‑93 (MQ=255)
                                           tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             >  2:3502695/1‑93 (MQ=255)
                                           tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             <  1:3144369/93‑1 (MQ=255)
                                           tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             <  1:1672701/93‑1 (MQ=255)
                                             tcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGt                                           >  2:1072803/1‑93 (MQ=255)
                                             tcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGt                                           >  2:2525144/1‑93 (MQ=255)
                                              cATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGtt                                          >  1:3095425/1‑92 (MQ=255)
                                              cATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTc                                          <  2:3536379/93‑1 (MQ=255)
                                               aTTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCt                                         >  2:3858959/1‑93 (MQ=255)
                                                ttATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg                                        >  2:544123/1‑93 (MQ=255)
                                                  aTCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg                                      <  1:3589399/93‑1 (MQ=255)
                                                  aTCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg                                      <  1:1384366/93‑1 (MQ=255)
                                                  aTCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg                                      >  1:1842928/1‑93 (MQ=255)
                                                  aTCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTggg                                      >  1:2574219/1‑93 (MQ=255)
                                                     tGGTGCTGAGAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAt                                   >  2:780201/1‑93 (MQ=255)
                                                      ggTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAtt                                  <  1:75200/93‑1 (MQ=255)
                                                      ggTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAtt                                  >  1:336176/1‑93 (MQ=255)
                                                        tGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGc                                <  1:1188880/93‑1 (MQ=255)
                                                        tGCTGAAAACAGCCACCACCATCATCACCGGCCCGGTCGCTTGGGCGATGCTCGCCGCGCTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGc                                <  2:1780973/93‑1 (MQ=255)
                                                         gCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGGTTGCt                               >  2:4013326/1‑93 (MQ=255)
                                                         gCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCt                               <  1:4115606/93‑1 (MQ=255)
                                                          cTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTg                              <  1:773792/93‑1 (MQ=255)
                                                          cTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTg                              <  1:3203648/93‑1 (MQ=255)
                                                          cTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTg                              >  1:682636/1‑93 (MQ=255)
                                                          cTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTg                              <  2:4295723/93‑1 (MQ=255)
                                                           tGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTgg                             <  2:1868583/93‑1 (MQ=255)
                                                            taaaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGATCGCCGCACTCGATGTTCGCTTAGCTTTTGTCTGGGGATTGCTGGc                            <  2:3800069/92‑1 (MQ=255)
                                                            gAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGc                            >  2:4144760/1‑93 (MQ=255)
                                                             aaaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGcc                           >  1:945445/1‑93 (MQ=255)
                                                             aaaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGcc                           <  1:843812/93‑1 (MQ=255)
                                                             aaaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGcc                           <  1:4298226/93‑1 (MQ=255)
                                                              aaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCt                          >  1:1357739/1‑93 (MQ=255)
                                                               aaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCtt                         <  2:1493770/93‑1 (MQ=255)
                                                                aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        >  1:4126440/1‑93 (MQ=255)
                                                                aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        <  2:3054379/93‑1 (MQ=255)
                                                                aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        <  1:1667416/93‑1 (MQ=255)
                                                                 cAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTg                       >  1:768473/1‑93 (MQ=255)
                                                                  aGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTgc                      >  2:2263603/1‑93 (MQ=255)
                                                                  aGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTgc                      >  1:3396268/1‑93 (MQ=255)
                                                                    ccATCAGCATCATCCCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTgcgc                    <  2:3326275/93‑1 (MQ=255)
                                                                     cATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCt                   >  1:217621/1‑93 (MQ=255)
                                                                     cATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCt                   <  1:3052422/93‑1 (MQ=255)
                                                                       tCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTa                 >  2:2168501/1‑93 (MQ=255)
                                                                        cAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTaa                <  1:2421169/93‑1 (MQ=255)
                                                                        cAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTaa                <  1:3839395/93‑1 (MQ=255)
                                                                        cAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTaa                >  1:3287876/1‑93 (MQ=255)
                                                                        cAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTaa                <  1:4134103/93‑1 (MQ=255)
                                                                         aGCATCATCACCGGGCCGGTCGCCTGGGCGATGCTCGCCGCCCTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAt               <  1:439257/93‑1 (MQ=255)
                                                                         aGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAt               <  1:1185226/93‑1 (MQ=255)
                                                                          gCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAtt              <  1:3509310/93‑1 (MQ=255)
                                                                           catcatCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTa             <  2:3702211/93‑1 (MQ=255)
                                                                           catcatCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTa             <  1:1344550/93‑1 (MQ=255)
                                                                             tcatcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACa           >  1:4031476/1‑93 (MQ=255)
                                                                             tcatcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACa           >  2:927437/1‑93 (MQ=255)
                                                                               atcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATc         <  2:1103934/93‑1 (MQ=255)
                                                                                tcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATcc        >  2:1635745/1‑93 (MQ=255)
                                                                                  aCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTACTTACATcccc      >  2:3911559/1‑92 (MQ=255)
                                                                                  aCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCg      <  2:1872247/93‑1 (MQ=255)
                                                                                   ccGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGa     >  1:3651895/1‑93 (MQ=255)
                                                                                   ccGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGa     >  1:1843292/1‑93 (MQ=255)
                                                                                   ccGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGa     >  1:1824953/1‑93 (MQ=255)
                                                                                         |                                                                                         
GGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAATA  >  NC_000913/1672473‑1672651

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