Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,310,728 T→C 78.6% I85V (ATT→GTT)  alsA ← D‑allose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,310,7280TC78.6% 412.7 / 103.2 205I85V (ATT→GTT) alsAD‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (74/85);  total (98/107)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 216 total aligned reads displayed.

CGATAATATTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCG  >  NC_000913/4310629‑4310825
                                                                                                   |                                                                                                 
cGATAATATTGACGCCACAGATTTTTTTCGTCAGATGACGACCACTATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGataaat                                                                                                   >  2:3367566/1‑100 (MQ=255)
   taataTTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAaacaat                                                                                                <  2:2677564/100‑1 (MQ=255)
   taataTTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTCTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAaacaat                                                                                                >  2:20391/1‑100 (MQ=255)
    aataTTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATc                                                                                               <  1:3879484/100‑1 (MQ=255)
     atatTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATcc                                                                                              <  1:2214425/100‑1 (MQ=255)
      tatTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATccc                                                                                             <  2:1866663/100‑1 (MQ=255)
            cGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATAcc                                                                                       <  1:1607249/100‑1 (MQ=255)
             gCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCg                                                                                      <  2:2147094/100‑1 (MQ=255)
               cacaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGCTTCCCTGATAAATAATCCCAATAccaa                                                                                     >  1:3826128/1‑97 (MQ=255)
               cacaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAg                                                                                    >  1:4083030/1‑100 (MQ=255)
               cacaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAg                                                                                    >  1:315797/1‑100 (MQ=255)
                acaGATTTTTTTTTTTTTTTGTCGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCgtgt                                                                                   >  2:1049882/1‑100 (MQ=255)
                acaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGt                                                                                   <  2:1134243/100‑1 (MQ=255)
                  aGATTTTTTTCGTCAGATGACGACCAATATATATATTTTCCATTACGGTTAATTCATCAATAACGCTGAGTACCTGACAAACAATCCCGATACCAAGtt                                                                                  <  2:3031216/99‑1 (MQ=255)
                  aGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTtg                                                                                 >  1:1040955/1‑100 (MQ=255)
                  aGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTtg                                                                                 <  1:999083/100‑1 (MQ=255)
                   gATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTtgt                                                                                >  1:3253236/1‑100 (MQ=255)
                     tttttttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTtgtg                                                                               <  2:579045/99‑1 (MQ=255)
                      ttttttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATACCGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGcc                                                                             <  1:2921700/100‑1 (MQ=255)
                      ttttttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGcc                                                                             <  1:1575588/100‑1 (MQ=255)
                         tttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCt                                                                          <  2:2338869/100‑1 (MQ=255)
                         tttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACACTGAGTTCTTGATAAACAATCCCGATACCGAGTTGTGCCGc                                                                           >  2:4021166/1‑99 (MQ=255)
                             gTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAAtt                                                                      <  1:3547538/100‑1 (MQ=255)
                             gTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAAAAATCCCGATACAGAGTTATGCCGCTAAtt                                                                      >  1:2033225/1‑100 (MQ=255)
                              tCAGATGACGACCAATGTATTTATTTTACAGTACGGTTAATTCATAAATACCACTGAGTTGGTGATAAACAATCCCGATGCCGAGTTGTGACGAGAAttt                                                                     <  2:201925/100‑1 (MQ=255)
                                aGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTAt                                                                   <  1:1730071/100‑1 (MQ=255)
                                aGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTAt                                                                   <  2:795674/100‑1 (MQ=255)
                                 gATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATg                                                                  <  2:163815/100‑1 (MQ=255)
                                 gATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATg                                                                  >  2:1349650/1‑100 (MQ=255)
                                 gATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATg                                                                  <  2:1589873/100‑1 (MQ=255)
                                   tGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGAt                                                                >  1:1775687/1‑100 (MQ=255)
                                    gacgacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATc                                                               <  2:2764109/100‑1 (MQ=255)
                                     acgacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATcc                                                              <  2:1296970/100‑1 (MQ=255)
                                     acgacCAATATATAAATTTTCCAGTACAGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATcc                                                              <  2:2544495/100‑1 (MQ=255)
                                      cgacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCa                                                             <  1:1761875/100‑1 (MQ=255)
                                      cgacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAAGCCCGATACCGATTTGTGCCGCTAATTTATGATCCa                                                             >  2:1802543/1‑100 (MQ=255)
                                        acCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGc                                                           >  2:448182/1‑100 (MQ=255)
                                         ccAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCt                                                          <  2:1820951/100‑1 (MQ=255)
                                         ccAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCTCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCt                                                          >  2:1911206/1‑100 (MQ=255)
                                         ccAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCt                                                          <  2:694174/100‑1 (MQ=255)
                                          cAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATGCCGAGTTGTGCCGCTAATTTATGATCCAGCtt                                                         >  2:3570171/1‑100 (MQ=255)
                                            atatatAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGt                                                       <  2:1953583/100‑1 (MQ=255)
                                              atataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTta                                                     >  2:3130383/1‑100 (MQ=255)
                                               tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                    >  1:1944160/1‑100 (MQ=255)
                                               tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                    <  2:240330/100‑1 (MQ=255)
                                               tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                    >  1:2302129/1‑100 (MQ=255)
                                               tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                    >  1:1169540/1‑100 (MQ=255)
                                               tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                    >  1:3528532/1‑100 (MQ=255)
                                                ataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtata                                                   <  2:928435/100‑1 (MQ=255)
                                                ataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtata                                                   >  2:116556/1‑100 (MQ=255)
                                                 taAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAg                                                  >  1:2867484/1‑100 (MQ=255)
                                                  aaaTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGc                                                 <  2:4064661/100‑1 (MQ=255)
                                                  aaaTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGc                                                 <  1:57751/100‑1 (MQ=255)
                                                         ccAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGtt                                          >  1:1990250/1‑100 (MQ=255)
                                                          cAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGtta                                         <  1:708095/100‑1 (MQ=255)
                                                          cAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGtta                                         >  2:515389/1‑100 (MQ=255)
                                                            gTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttatt                                       <  2:3670229/100‑1 (MQ=255)
                                                            gTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttatt                                       >  1:2535668/1‑100 (MQ=255)
                                                              aCGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaa                                     >  2:2114204/1‑100 (MQ=255)
                                                              aCGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaa                                     <  1:781160/100‑1 (MQ=255)
                                                               cGGTTAATNCATCAATAACGCTGCGTNCATGNTANACAATCCCGATACCGAGTTGTGGCGCTATTTTATGATCCAGCTTGTTATAGCTAATGTNATtcat                                    >  1:3379483/1‑100 (MQ=255)
                                                                ggTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatg                                   >  2:242833/1‑100 (MQ=255)
                                                                 gTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatgg                                  >  1:1372027/1‑100 (MQ=255)
                                                                  ttAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggt                                 <  2:1144533/100‑1 (MQ=255)
                                                                   tAATTCATCAATAACGCTGAGTTCCTGATAAACATTCCCGATACCGAGTTGTGCAGCTAATTTAAGATCCAGCTTGTTATAGCTAATGTTATtaatggta                                <  2:2271443/100‑1 (MQ=255)
                                                                    aaTTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTACAGCTAATGTTATtaatgctaa                               >  2:15186/1‑100 (MQ=255)
                                                                      ttCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatg                             >  1:813337/1‑100 (MQ=255)
                                                                      ttCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatg                             <  2:383099/100‑1 (MQ=255)
                                                                      ttCATCAANAACGCTGAGTTCCTGATNAATANTCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtagtngtaatg                             >  1:3423322/1‑100 (MQ=255)
                                                                        catcaATAACGCTGAGTTCCTGATACATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatgg                            >  2:841061/1‑99 (MQ=255)
                                                                        catcaATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTaatggtaatggt                           <  1:106184/100‑1 (MQ=255)
                                                                        catcaATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCAGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTCTtaatggtaatgg                            >  2:3327957/1‑99 (MQ=255)
                                                                         atcaATACCGCAGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTg                          >  2:2188541/1‑100 (MQ=255)
                                                                         atcaATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTg                          >  1:2597325/1‑100 (MQ=255)
                                                                          tcaATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGc                         >  1:420191/1‑100 (MQ=255)
                                                                          tcaATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTCTAGCTAATGTTATTAATGGTAATGGTGc                         >  2:3038518/1‑100 (MQ=255)
                                                                           caATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGcc                        >  2:2788372/1‑100 (MQ=255)
                                                                           caATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATAATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGcc                        >  2:4170744/1‑100 (MQ=255)
                                                                            aaTAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCt                       >  2:3053893/1‑100 (MQ=255)
                                                                             aTAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCtt                      >  1:3616247/1‑100 (MQ=255)
                                                                              tAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCttt                     <  2:1852653/100‑1 (MQ=255)
                                                                                aCGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTgg                   <  2:3905146/100‑1 (MQ=255)
                                                                                   cTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTc                 <  1:1189592/99‑1 (MQ=255)
                                                                                   cTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCg                <  2:2361755/100‑1 (MQ=255)
                                                                                   cTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCg                <  1:1668827/100‑1 (MQ=255)
                                                                                     gAGTTCCTGATAAACAATCCCGATACCGAGTTGTGTCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc              <  2:226067/100‑1 (MQ=255)
                                                                                     gAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc              <  2:3264389/100‑1 (MQ=255)
                                                                                     gAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc              <  1:2059851/100‑1 (MQ=255)
                                                                                       gTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTc            >  1:2451249/1‑100 (MQ=255)
                                                                                          ccTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg         <  1:1105963/100‑1 (MQ=255)
                                                                                          ccTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg         <  2:2973496/100‑1 (MQ=255)
                                                                                           cTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt        <  1:3947572/100‑1 (MQ=255)
                                                                                           cTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt        >  2:3342245/1‑100 (MQ=255)
                                                                                           cTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt        >  2:3155744/1‑100 (MQ=255)
                                                                                           cTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt        <  1:2707386/100‑1 (MQ=255)
                                                                                            tGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTa       <  2:740147/100‑1 (MQ=255)
                                                                                             gATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTAt      >  1:461640/1‑100 (MQ=255)
                                                                                             gATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTAt      >  2:2578065/1‑100 (MQ=255)
                                                                                               taaataaTCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTAtt     >  2:3220665/1‑99 (MQ=255)
                                                                                                 aataatCCCGATACCGTGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCg  <  1:1022740/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CGATAATATTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCG  >  NC_000913/4310629‑4310825

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