Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,667,467 G→T 6.8% R38R (CGC→CGA gadA ← glutamate decarboxylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,667,4670GT 6.8% 346.2 / 9.5 103R38R (CGC→CGAgadAglutamate decarboxylase A
Reads supporting (aligned to +/- strand):  ref base G (49/47);  new base T (3/4);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 105 total aligned reads displayed.

CGATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAAT  >  NC_000913/3667325‑3667611
                                                                                                                                              |                                                                                                                                                
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:690904/150‑1 (MQ=255)
                   cGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGttt                                                                                                                        <  1:234937/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)
                           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)
                                                                  aGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGa                                                                         >  1:92217/1‑150 (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)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCTCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                             <  1:1497209/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                             <  2:1162376/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                             <  1:270107/150‑1 (MQ=255)
                                                                              ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGc                                                             <  1:208872/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)
                                                                                                              tcatcaTTGATAATCTGAAATGCGACATCATCGCGCATTTGGTGCTGCTGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTTCGCCAAAACGTTAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCAtttag                             <  2:186297/150‑3 (MQ=255)
                                                                                                              tcatcaTTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCg                             <  2:601359/150‑1 (MQ=255)
                                                                                                                      taTAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCtt                     <  1:1071929/149‑1 (MQ=255)
                                                                                                                      gATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCa                                  <  2:17731/137‑1 (MQ=255)
                                                                                                                      gATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCa                                  >  1:17731/1‑137 (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           <  2:1327668/150‑1 (MQ=255)
                                                                                                                                aaTGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAatttattt           <  1:777729/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)
                                                                                                                                         tcatcGCGCATTTCGTGCAGCGGAAATCGTTTTGCCTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCGTTAAATTTATTTGAAGGCAAt  >  2:1212594/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGATCCAGTTTTTATTGATCGACAAATCCATCAATTTATGGACGTTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCATCGCGCATTTCGTGCAGCGGAAATCGTTTTGACTCCGCGATAGTAGAAATGGCCTTTGCGCCAAAACGTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAAT  >  NC_000913/3667325‑3667611

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: