Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,441,062 G→A 10.2% S205F (TCC→TTC) ‡ rpsD ← 30S ribosomal subunit protein S4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,441,0620GA10.2% 538.6 / 38.7 176S205F (TCC→TTC) ‡rpsD30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base G (80/78);  new base A (6/12);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 177 total aligned reads displayed.

GAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440915‑3441209
                                                                                                                                                   |                                                                                                                                                   
gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  1:1007357/149‑1 (MQ=255)
gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  2:1018964/149‑1 (MQ=255)
gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  1:542372/149‑1 (MQ=255)
gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  4:217133/149‑1 (MQ=255)
gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  4:929290/149‑1 (MQ=255)
    aTGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTaa                                                                                                                                                <  4:114327/149‑1 (MQ=255)
       gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGc                                                                                                                                             <  4:1042801/149‑1 (MQ=255)
        ccAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCt                                                                                                                                            <  2:1095176/149‑1 (MQ=255)
         cAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtt                                                                                                                                           <  4:438786/149‑2 (MQ=255)
               cACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGAt                                                                                                                                     <  4:433097/149‑1 (MQ=255)
               cACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGAt                                                                                                                                     <  1:787779/149‑1 (MQ=255)
                aCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGTATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt                                                                                                                                            <  1:36069/141‑1 (MQ=255)
                aCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt                                                                                                                                            >  2:36069/1‑141 (MQ=255)
                 cGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCa                                                                                                                                   <  2:160852/149‑1 (MQ=255)
                  gCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAg                                                                                                                                  <  1:1050014/149‑1 (MQ=255)
                      tAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgtgt                                                                                                                              >  2:109080/1‑149 (MQ=255)
                        aaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtta                                                                                                                                          <  4:271671/135‑3 (MQ=255)
                             tCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCt                                                                                                                                            <  4:1382650/128‑1 (MQ=255)
                               aaGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATg                                                                                                                     >  4:122216/1‑149 (MQ=255)
                                aGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGt                                                                                                                    >  2:1335804/1‑149 (MQ=255)
                                   gTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTaaa                                                                                                                                               <  1:1254352/119‑1 (MQ=255)
                                   gTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTaaa                                                                                                                                               >  2:1254352/1‑119 (MQ=255)
                                    tCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGc                                                                                                                <  2:246799/149‑1 (MQ=255)
                                        cTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagac                                                                                                            >  1:1384111/1‑149 (MQ=255)
                                        cTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagac                                                                                                            <  4:981091/149‑1 (MQ=255)
                                         ttGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagaca                                                                                                           >  1:294887/1‑149 (MQ=255)
                                            gcgtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGAt                                                                                                        <  2:1384111/149‑1 (MQ=255)
                                              gtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCATAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTAGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCa                                                                                                      <  4:904954/149‑1 (MQ=255)
                                                      aaCTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTc                                                                                              <  2:1147760/149‑1 (MQ=255)
                                                         tCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCgg                                                                                           >  1:1010859/1‑149 (MQ=255)
                                                         tCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtta                                                                                                                                          >  4:1364296/1‑100 (MQ=255)
                                                         tCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtta                                                                                                                                          <  3:1364296/102‑3 (MQ=255)
                                                              tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTAc                                                                                      >  4:524522/1‑149 (MQ=255)
                                                               gCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  1:1239978/86‑1 (MQ=255)
                                                                 tCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcaga                                                                                                             >  1:1077103/1‑123 (MQ=255)
                                                                    atatTAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    >  1:467086/1‑81 (MQ=255)
                                                                    atatTAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                    <  2:467086/81‑1 (MQ=255)
                                                                         aaCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           >  3:874138/1‑149 (MQ=255)
                                                                          aCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          >  3:713182/1‑149 (MQ=255)
                                                                             aGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCa                                                                                                      <  1:1453678/118‑1 (MQ=255)
                                                                             aGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCtt                                                                       >  2:941046/1‑149 (MQ=255)
                                                                              ggCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTc                                                                      <  3:1319608/149‑1 (MQ=255)
                                                                              ggCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTc                                                                      <  3:122216/149‑1 (MQ=255)
                                                                                     ttttAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                        <  2:1052544/76‑1 (MQ=255)
                                                                                        tAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCa                                                            <  3:524522/149‑1 (MQ=255)
                                                                                           aaaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTc                                                                                                                   <  2:457849/91‑1 (MQ=255)
                                                                                           aaaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           <  1:1094221/131‑1 (MQ=255)
                                                                                           aaaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGcc                                                             >  1:1307016/1‑145 (MQ=255)
                                                                                            aaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagaca                                                                                                           >  1:1105986/1‑98 (MQ=255)
                                                                                              ctctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGc                                                                                    >  3:1466734/1‑119 (MQ=255)
                                                                                               tctGTGACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGATTAAACGTACCTTCCATCTTGCCAGCATCaa                                                     >  1:415551/1‑149 (MQ=255)
                                                                                               tctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCaa                                                     >  1:493911/1‑149 (MQ=255)
                                                                                               tctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCaa                                                     <  4:713182/149‑1 (MQ=255)
                                                                                                   tCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCt                                                                                                                                            >  1:287083/1‑58 (MQ=37)
                                                                                                   tCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCt                                                                                                                                            <  2:287083/58‑1 (MQ=37)
                                                                                                      caGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccag                                              <  3:920266/149‑1 (MQ=255)
                                                                                                      caGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGAAAGATCa                                                                                                      >  3:904311/1‑93 (MQ=255)
                                                                                                       aGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTc                                                                                              >  2:883565/1‑100 (MQ=255)
                                                                                                           cccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTTGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          <  4:850033/116‑1 (MQ=255)
                                                                                                           cccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          >  4:32952/1‑116 (MQ=255)
                                                                                                           cccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGc                                                          <  3:1265379/132‑1 (MQ=255)
                                                                                                           cccTGCATTATGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          <  3:32952/116‑1 (MQ=255)
                                                                                                             cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTc                                                                                                                   >  3:262425/1‑73 (MQ=255)
                                                                                                             cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTaaa                                                                               >  2:619317/1‑109 (MQ=255)
                                                                                                             cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           >  1:188863/1‑113 (MQ=255)
                                                                                                             cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGtt                                       <  1:933601/149‑1 (MQ=255)
                                                                                                              tGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGc                                                                                          >  2:293416/1‑97 (MQ=255)
                                                                                                                  ttGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGt                                        >  2:749725/1‑143 (MQ=255)
                                                                                                                     tgtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGGCAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCTGCCAGGTTGGCTTTTc                               <  4:623953/149‑1 (MQ=255)
                                                                                                                      gtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCa                              <  2:415551/149‑1 (MQ=255)
                                                                                                                       tCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAAt                                                                                                                      >  3:1122959/1‑60 (MQ=255)
                                                                                                                       tCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCAc                             >  1:1354471/1‑149 (MQ=255)
                                                                                                                       tCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCAc                             >  1:930424/1‑149 (MQ=255)
                                                                                                                         ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          >  4:1159507/1‑102 (MQ=255)
                                                                                                                         ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          <  3:1159507/102‑1 (MQ=255)
                                                                                                                         ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                           <  2:493911/149‑1 (MQ=255)
                                                                                                                         ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                           <  2:1405094/149‑1 (MQ=255)
                                                                                                                          tctctTTGGTGCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           <  2:1276680/100‑1 (MQ=255)
                                                                                                                           ctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCGTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctg                         <  3:996519/149‑1 (MQ=255)
                                                                                                                            tctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCt                                   <  4:1101897/138‑1 (MQ=255)
                                                                                                                            tctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCt                                   >  3:1101897/1‑138 (MQ=255)
                                                                                                                            tctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgc                        >  3:1275971/1‑149 (MQ=255)
                                                                                                                              tttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTc                      >  1:972361/1‑149 (MQ=255)
                                                                                                                              tttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTc                      >  2:1130647/1‑149 (MQ=255)
                                                                                                                              tttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCAAGGTTGGCTTTTCACGCTGCTc                      >  1:1253967/1‑149 (MQ=255)
                                                                                                                                   tACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTAc                                                                                      <  1:881358/80‑1 (MQ=255)
                                                                                                                                   tACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTAc                                                                                      >  2:881358/1‑80 (MQ=255)
                                                                                                                                      tAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcc                                                <  3:1376796/115‑1 (MQ=255)
                                                                                                                                        aGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGca                                                                                                               >  3:488805/1‑50 (MQ=255)
                                                                                                                                        aGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc            >  2:1050365/1‑149 (MQ=255)
                                                                                                                                         gCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCa           >  1:1429/1‑149 (MQ=255)
                                                                                                                                          cTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagc                   >  4:1294782/1‑140 (MQ=255)
                                                                                                                                          cTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagc                   <  1:701628/140‑1 (MQ=255)
                                                                                                                                          cTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCag          >  4:1226326/1‑149 (MQ=255)
                                                                                                                                           tttACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCaga         <  2:930424/149‑1 (MQ=255)
                                                                                                                                            ttACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCagag        >  4:1055163/1‑149 (MQ=255)
                                                                                                                                                ttGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTCCCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGtt                                       <  4:1386444/114‑1 (MQ=255)
                                                                                                                                                ttGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGc    >  1:1174123/1‑149 (MQ=255)
                                                                                                                                                ttGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGACAGCATCAACTTCCAGCCAGGtt                                       >  3:1386444/1‑114 (MQ=255)
                                                                                                                                                  ggAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCtt  >  3:580647/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440915‑3441209

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