Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,510,339 T→C V221A (GTG→GCG)  b1440 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3390TC97.9% 266.4 / ‑2.9 94V221A (GTG→GCG) b1440putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (47/44);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGTG  >  NC_000913/1510250‑1510428
                                                                                         |                                                                                         
tGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCt                                                                                        <  2:4569027/93‑1 (MQ=255)
tGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCAGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCt                                                                                        >  1:3202791/1‑93 (MQ=255)
    gACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaaa                                                                                    <  1:2755630/93‑1 (MQ=255)
     aCGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAg                                                                                   <  2:2411275/93‑1 (MQ=255)
      cGCTGCGTGGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGt                                                                                  <  1:949974/93‑1 (MQ=255)
         tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                               >  2:3220126/1‑93 (MQ=255)
         tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                               >  1:2167140/1‑93 (MQ=255)
          gCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctg                                                                              >  2:753554/1‑93 (MQ=255)
           ggttgttCGTTAAAGCCGCTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctgc                                                                             <  2:1175398/92‑1 (MQ=255)
             ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                           >  2:1874496/1‑93 (MQ=255)
             ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                           <  2:3283001/93‑1 (MQ=255)
             ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                           >  1:3237057/1‑93 (MQ=255)
               gttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgcgc                                                                         >  1:529966/1‑93 (MQ=255)
               gttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgcgc                                                                         <  1:3717247/93‑1 (MQ=255)
                ttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCt                                                                        >  2:2442452/1‑93 (MQ=255)
                 tCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTc                                                                       >  1:686997/1‑93 (MQ=255)
                    ttAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAac                                                                    >  1:3275874/1‑93 (MQ=255)
                      aaaGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAacac                                                                  >  2:2547364/1‑93 (MQ=255)
                       acgCCTCTCAGCCGAAGTAGGGCATCAAAAGTCCGTATCAACTCACCGAAGAACAGTACAAGGCGGCGCTACAAGTGCTGCGCGCTCAAcaca                                                                 <  2:4682897/91‑1 (MQ=255)
                        aGCCCCTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAg                                                                <  2:47821/93‑1 (MQ=255)
                        aGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAg                                                                >  1:1482880/1‑93 (MQ=255)
                        aGCCACTCAGCCGCAGGTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAg                                                                <  1:714039/93‑1 (MQ=255)
                        aGCCACGCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAg                                                                <  2:1385909/93‑1 (MQ=255)
                         gCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGt                                                               <  2:1886689/93‑1 (MQ=255)
                          ccACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtt                                                              >  2:4115630/1‑93 (MQ=255)
                          ccACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtt                                                              <  1:2465290/93‑1 (MQ=255)
                          ccAATCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtt                                                              <  2:2903542/93‑1 (MQ=255)
                           cACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGttt                                                             >  1:1622749/1‑93 (MQ=255)
                            aCTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtgtg                                                            >  2:758252/1‑93 (MQ=255)
                             cTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGa                                                           <  1:2168881/93‑1 (MQ=255)
                             cTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGa                                                           <  2:3035560/93‑1 (MQ=255)
                                 gCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCa                                                       >  2:548001/1‑93 (MQ=255)
                                 gCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCa                                                       >  2:3661777/1‑93 (MQ=255)
                                   cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                     <  2:4058265/93‑1 (MQ=255)
                                    gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                    >  1:3969530/1‑93 (MQ=255)
                                     cAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCCTCGc                                                   >  2:1496677/1‑93 (MQ=255)
                                       gTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTa                                                 <  1:3654046/93‑1 (MQ=255)
                                       gTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTa                                                 <  1:3654045/93‑1 (MQ=255)
                                        ttGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                                <  1:3839456/93‑1 (MQ=255)
                                         tGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACt                                               >  2:2722207/1‑93 (MQ=255)
                                          gggCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTg                                              >  1:442031/1‑93 (MQ=255)
                                          gggCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTg                                              <  1:1959994/93‑1 (MQ=255)
                                           gccATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTgg                                             <  1:361366/91‑1 (MQ=255)
                                             cATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCa                                           <  2:3965169/93‑1 (MQ=255)
                                              aTCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGATTGATCCATCGCTACTGGCAt                                          >  2:2557733/1‑93 (MQ=255)
                                               tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                         >  2:2651516/1‑93 (MQ=255)
                                               tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                         >  2:1198600/1‑93 (MQ=255)
                                                cAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGa                                        <  1:4083793/93‑1 (MQ=255)
                                                  gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACtgacatgaca                                      >  1:4542873/1‑93 (MQ=255)
                                                  gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGaca                                      >  1:4108571/1‑93 (MQ=255)
                                                    gATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACt                                    >  2:203120/1‑93 (MQ=255)
                                                    gATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACt                                    >  1:1148825/1‑93 (MQ=255)
                                                     aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                   <  1:995310/93‑1 (MQ=255)
                                                     aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                   >  1:2815081/1‑93 (MQ=255)
                                                     aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                   >  1:1750467/1‑93 (MQ=255)
                                                     aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                   >  1:3521790/1‑93 (MQ=255)
                                                       ccGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAcc                                 <  1:1630090/93‑1 (MQ=255)
                                                        cGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCg                                <  1:1297360/93‑1 (MQ=255)
                                                         ggaTCAACTCACCGAAGAACAGTACCAGGCGGAGCTGAAAGTGCTGCGCGCTCAACACAGCTTGATCGATCGCTACTGGCATGACACTAcccg                               >  2:1260489/3‑91 (MQ=255)
                                                         gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                               >  2:1566415/1‑93 (MQ=255)
                                                         gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                               >  1:2458247/1‑93 (MQ=255)
                                                         gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                               >  1:3361320/1‑93 (MQ=255)
                                                         gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                               <  1:4604044/93‑1 (MQ=255)
                                                          tATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTg                              >  1:2460596/1‑93 (MQ=255)
                                                          tATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTg                              >  1:3406991/1‑93 (MQ=255)
                                                            tCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCa                            >  1:3571468/1‑93 (MQ=255)
                                                             cAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaa                           >  2:1506922/1‑93 (MQ=255)
                                                             cAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaa                           <  1:584571/93‑1 (MQ=255)
                                                              aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                          <  1:4587682/93‑1 (MQ=255)
                                                              aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                          <  2:4308437/93‑1 (MQ=255)
                                                              aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                          <  2:3147677/93‑1 (MQ=255)
                                                              aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                          <  1:52998/93‑1 (MQ=255)
                                                               aCTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                         >  2:3547370/1‑93 (MQ=255)
                                                               aCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCCAAt                         >  2:2266274/1‑93 (MQ=255)
                                                               aCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                         <  2:35059/93‑1 (MQ=255)
                                                               aCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                         <  1:4150823/93‑1 (MQ=255)
                                                               aCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                         >  2:3673858/1‑93 (MQ=255)
                                                                 tCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTCCCGTCCCAATGa                       >  2:3156170/1‑93 (MQ=255)
                                                                 tCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGa                       >  2:1667550/1‑93 (MQ=255)
                                                                   aCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGGAACTACCCTGCAAAtaggg                     >  1:1127798/1‑89 (MQ=255)
                                                                   aCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGc                     <  2:954875/93‑1 (MQ=255)
                                                                      gaagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAt                  <  1:674766/93‑1 (MQ=255)
                                                                      gaagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAt                  >  2:4036124/1‑93 (MQ=255)
                                                                         gaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTc               <  1:876104/93‑1 (MQ=255)
                                                                         gaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTc               >  2:3831306/1‑93 (MQ=255)
                                                                         gaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTc               <  1:2322608/93‑1 (MQ=255)
                                                                           aCAGTACCCGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaa             >  2:4541737/1‑93 (MQ=255)
                                                                           aCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaa             >  1:931413/1‑93 (MQ=255)
                                                                            cAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaa            <  2:4143077/93‑1 (MQ=255)
                                                                             aGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaa           <  2:994574/93‑1 (MQ=255)
                                                                              gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCcaaaa          >  1:4384843/1‑93 (MQ=255)
                                                                              gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa          <  1:4316648/93‑1 (MQ=255)
                                                                               tACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAAc         <  1:3155799/93‑1 (MQ=255)
                                                                                 ccAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGa       <  2:363364/93‑1 (MQ=255)
                                                                                    ggcggcGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAggg    <  2:345940/93‑1 (MQ=255)
                                                                                     gcggcgCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGgt   >  2:3815951/1‑93 (MQ=255)
                                                                                      cGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGgtg  >  2:1014304/1‑93 (MQ=255)
                                                                                         |                                                                                         
TGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGTG  >  NC_000913/1510250‑1510428

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