Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,667,460 C→G 6.8% V41L (GTC→CTC) ‡ gadA ← glutamate decarboxylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,667,4600CG 6.8% 344.7 / 12.2 103V41L (GTC→CTC) ‡gadAglutamate decarboxylase A
Reads supporting (aligned to +/- strand):  ref base C (48/48);  new base G (3/4);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 104 total aligned reads displayed.

CGATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCA  >  NC_000913/3667325‑3667609
                                                                                                                                       |                                                                                                                                                     
cGATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCAttt                                                                                                                                         <  1:300110/150‑1 (MQ=255)
 gATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTc                                                                                                                                        >  1:1169587/1‑150 (MQ=255)
 gATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGAGACATCATCGCGCATTTc                                                                                                                                        >  1:1163771/1‑150 (MQ=255)
 gATCCAGGTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTc                                                                                                                                        <  1:1595483/150‑1 (MQ=255)
  aTCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGCGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGCGATGATGTCGCATTTca                                                                                                                                       <  1:754107/150‑2 (MQ=17)
   tCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGCGATGATGTCGCATTTcag                                                                                                                                      <  2:1339863/150‑3 (MQ=17)
    ccAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGCGATGATGTCGCATTTcaga                                                                                                                                     <  2:1461286/150‑4 (MQ=17)
    ccAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTg                                                                                                                                     >  2:1425081/1‑150 (MQ=255)
     cAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGc                                                                                                                                    >  2:1134394/1‑150 (MQ=255)
              ttGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAAt                                                                                                                           >  2:710564/1‑150 (MQ=255)
              ttGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAAt                                                                                                                           >  1:1255352/1‑150 (MQ=255)
                gATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCGGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCg                                                                                                                         >  2:1541840/1‑150 (MQ=255)
                gATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCg                                                                                                                         >  1:1600938/1‑150 (MQ=255)
                  tCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGtt                                                                                                                       >  2:130092/1‑150 (MQ=255)
                   cGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGttt                                                                                                                      <  1:234937/150‑1 (MQ=255)
                   cGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGttt                                                                                                                      <  1:690904/150‑1 (MQ=255)
                    gACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGtttt                                                                                                                     <  1:370139/150‑1 (MQ=255)
                     aCAAATCCATCAATTTATGGACGTGTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTATGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTg                                                                                                                    <  2:1422687/150‑1 (MQ=255)
                      cAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGa                                                                                                                   >  1:308444/1‑150 (MQ=255)
                       aaaTCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGAc                                                                                                                  <  1:872785/150‑1 (MQ=255)
                        aaTCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACt                                                                                                                 >  1:1589248/1‑150 (MQ=255)
                           ccATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCcg                                                                                                              >  1:1162376/1‑150 (MQ=255)
                           ccATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTAGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTCCAGCGGAAATCGTTTTGACTccc                                                                                                              >  1:1110338/1‑149 (MQ=255)
                            cATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCcgc                                                                                                             >  2:208872/1‑150 (MQ=255)
                             aTCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCcgcg                                                                                                            >  1:746203/1‑150 (MQ=255)
                              tCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGa                                                                                                           >  1:320792/1‑150 (MQ=255)
                                aaTTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGAta                                                                                                         >  1:328363/1‑150 (MQ=255)
                                aaTTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGAta                                                                                                         >  1:1336555/1‑150 (MQ=255)
                                 aTTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATGATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGAtag                                                                                                        <  2:383142/150‑1 (MQ=255)
                                    tATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCATGTTCTGTCGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCTGAAATCTTTTTGACTCCGCGAtagtag                                                                                                     <  2:504471/150‑1 (MQ=255)
                                    tATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGAtagtag                                                                                                     >  1:997031/1‑150 (MQ=255)
                                       ggACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGaaa                                                                                                  >  1:1600416/1‑150 (MQ=255)
                                        gACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAAt                                                                                                 >  1:1560027/1‑150 (MQ=255)
                                          cGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATgg                                                                                               <  1:1006263/150‑1 (MQ=255)
                                                  tcgtcCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTgcg                                                                                       <  2:112888/150‑1 (MQ=255)
                                                   cgtcCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTgcgc                                                                                      >  2:761392/1‑150 (MQ=255)
                                                     tcCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCGCTTgcgcaa                                                                                    >  2:23291/1‑148 (MQ=255)
                                                      cccAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCaa                                                                                   >  2:1005890/1‑150 (MQ=255)
                                                      cccAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCaa                                                                                   >  2:676660/1‑150 (MQ=255)
                                                        cAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGTCaaaa                                                                                 >  1:1309580/1‑150 (MQ=255)
                                                         aGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAAc                                                                                <  2:328363/150‑1 (MQ=255)
                                                           gTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGt                                                                              <  1:493817/150‑1 (MQ=255)
                                                            tCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTg                                                                             <  2:1083444/150‑1 (MQ=255)
                                                             cTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCgcgc                                                                                                                                             <  1:241041/85‑1 (MQ=255)
                                                             cTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCgcgc                                                                                                                                             >  2:241041/1‑85 (MQ=255)
                                                             cTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGa                                                                            >  1:720875/1‑150 (MQ=255)
                                                               ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAAt                                                                          >  2:71585/1‑150 (MQ=255)
                                                                 cAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGTAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCg                                                                        <  1:920414/150‑1 (MQ=255)
                                                                 cAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCg                                                                        <  1:1343405/150‑1 (MQ=255)
                                                                     aaGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGAGCCAAAACGTGAATCGagta                                                                    >  1:1327668/1‑150 (MQ=255)
                                                                      aGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGagtag                                                                   <  2:1486392/150‑1 (MQ=255)
                                                                       gTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGagtagt                                                                  <  1:1265187/150‑1 (MQ=255)
                                                                         ggCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTc                                                                <  1:267587/150‑1 (MQ=255)
                                                                         ggCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTAATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGATTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATAGAGTAGttg                                                                >  2:491068/1‑149 (MQ=255)
                                                                          gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCt                                                               >  2:830672/1‑150 (MQ=255)
                                                                          gCCAGGTTCTGACGAGCGGTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCt                                                               <  1:700176/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCTCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                           <  1:1497209/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                           <  1:208872/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                           <  2:1162376/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                           <  1:270107/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGCTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                           <  1:1134394/150‑1 (MQ=255)
                                                                               gTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCg                                                          <  1:1541840/150‑1 (MQ=255)
                                                                                  cTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGaa                                                       <  2:95073/150‑1 (MQ=255)
                                                                                  cTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGaa                                                       >  2:1329070/1‑150 (MQ=255)
                                                                                   tGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGaaa                                                      >  2:1002626/1‑150 (MQ=255)
                                                                                        aGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGt                                                 >  2:1513344/1‑150 (MQ=255)
                                                                                         gCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGtt                                                <  2:746203/150‑1 (MQ=255)
                                                                                          cGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTa                                               <  1:363150/150‑1 (MQ=255)
                                                                                          cGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTa                                               <  2:1600938/150‑1 (MQ=255)
                                                                                           gTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTaa                                              >  2:1228226/1‑150 (MQ=255)
                                                                                            ttGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAAc                                             <  2:1171611/150‑1 (MQ=255)
                                                                                             tGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACa                                            >  2:891407/1‑150 (MQ=255)
                                                                                              gCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAg                                           <  1:1066278/150‑1 (MQ=255)
                                                                                                 aTCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCtt                                        >  2:777729/1‑150 (MQ=255)
                                                                                                  tCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTc                                       >  2:615528/1‑150 (MQ=255)
                                                                                                   cAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCt                                      <  1:761392/150‑1 (MQ=255)
                                                                                                       atatAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTc                                  <  1:1005890/150‑1 (MQ=255)
                                                                                                        tataATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTcc                                 >  1:1030103/1‑150 (MQ=255)
                                                                                                          taATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCAt                               <  2:200508/150‑1 (MQ=255)
                                                                                                          taATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCAt                               >  1:825400/1‑150 (MQ=255)
                                                                                                              tcatcaTTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCg                           <  2:601359/150‑1 (MQ=255)
                                                                                                                atcatTGATAATCTGAAATGCGACATCAACGCGCATTTCGTCCAGCCGAAATCCTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGaa                         >  2:127373/1‑150 (MQ=255)
                                                                                                                      taTAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCtt                   <  1:1071929/149‑1 (MQ=255)
                                                                                                                      gATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCa                                <  2:17731/137‑1 (MQ=255)
                                                                                                                       aTAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTa                  <  2:320792/150‑1 (MQ=255)
                                                                                                                        tAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTaa                 <  1:1267732/150‑1 (MQ=255)
                                                                                                                        ctgcacgAAATGCGCGATGATGTCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTaa                 >  2:490142/7‑150 (MQ=25)
                                                                                                                          gcacgAAATGCGCGATGATGTCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAAc                        >  2:347067/5‑141 (MQ=25)
                                                                                                                          gcacgAAATGCGCGATGATGTCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAAc                        <  1:347067/137‑1 (MQ=25)
                                                                                                                           tCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAatt              >  2:705449/1‑150 (MQ=255)
                                                                                                                             tGAAATGCGACATCATCGCGCATTTAGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAattta            >  2:899590/1‑150 (MQ=255)
                                                                                                                             cgAAATGCGCGATGATGTCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAattta            >  1:1623865/2‑150 (MQ=25)
                                                                                                                                aaTGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAatttattt         <  1:777729/150‑1 (MQ=255)
                                                                                                                                aaTGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAatttattt         <  2:1327668/150‑1 (MQ=255)
                                                                                                                                 aTGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTg        <  2:166616/150‑1 (MQ=255)
                                                                                                                                 aTGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTg        <  2:1560027/150‑1 (MQ=255)
                                                                                                                                  tGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGa       >  2:913737/1‑150 (MQ=255)
                                                                                                                                     gACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGTCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAgg    <  2:1309580/150‑1 (MQ=255)
                                                                                                                                     gACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAgg    <  1:676660/150‑1 (MQ=255)
                                                                                                                                       catcatCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCa  >  1:798858/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
CGATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCA  >  NC_000913/3667325‑3667609

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: