Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 329,512 T→A 5.6% V17E (GTG→GAG)  betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,5120TA 5.6% 373.9 / 8.8 135V17E (GTG→GAG) betTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base T (69/55);  new base A (2/6);  total (71/64)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 140 total aligned reads displayed.

ATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATT  >  NC_000913/329364‑329660
                                                                                                                                                    |                                                                                                                                                    
aTCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCggtggt                                                                                                                                                      <  2:267513/149‑1 (MQ=255)
  cTATTAAAGCGGTTAATGAGTGCACGTGCAATATAAAATGTCTCTTCACTGCGCCGAATTAGAGTTTACATAGTACCAATACCAGCGGGGCTACCGGATGCCAGACCTTTCACACAGCAGGGACCAGGACACCAGCCATCCGGTGgtgt                                                                                                                                                    <  2:876049/149‑1 (MQ=255)
   tATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGtt                                                                                                                                                   >  1:1426338/1‑149 (MQ=255)
    aTTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGttt                                                                                                                                                  >  1:1415126/1‑149 (MQ=255)
    aTTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGAACTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGttt                                                                                                                                                  >  1:1388443/1‑149 (MQ=255)
     ttAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCAAACAGCAGGGAAAAGTACAAAATAAATCAGGTGGTGtttt                                                                                                                                                 >  1:264976/1‑149 (MQ=255)
      tAAAGCGGTTATTGATTGGAAGTTCAATATAATATGTGTCTTAATTGTTACGATTTTGATTTTAAATATTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAATGGACAAAATCAATCCGGTGGTGTTTTa                                                                                                                                                <  1:62729/149‑1 (MQ=255)
        aaGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACATTAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTaca                                                                                                                                              >  1:39485/1‑149 (MQ=255)
        aaGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTaca                                                                                                                                              <  1:142997/149‑1 (MQ=255)
          gCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACAcc                                                                                                                                            >  2:817838/1‑149 (MQ=255)
           cGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTATACACCt                                                                                                                                           >  1:1036184/1‑149 (MQ=255)
             gTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTcc                                                                                                                                         >  2:1289638/1‑149 (MQ=255)
              ttATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTccg                                                                                                                                        >  2:1032599/1‑149 (MQ=255)
                 ttgattgGACGTGCCATATATCAGGTGTCTTTATTGTTACGACTTTGATTGTAAATCGGACCACTAACAGTGTGGCTACTGGATGCCAGCCCTTTCCCACTGCAGGGACAAGGACAAAATCAATCCGGTGGGGTGTTACCCCTccgccg                                                                                                                                     <  2:1417928/149‑1 (MQ=255)
                    attgGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGAc                                                                                                                                  <  1:691947/149‑1 (MQ=255)
                    attgGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGAc                                                                                                                                  >  1:1476480/1‑149 (MQ=255)
                       ggATGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGa                                                                                                                               >  2:7190/1‑149 (MQ=255)
                       ggACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGa                                                                                                                               >  2:533490/1‑149 (MQ=255)
                        gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGTAAAAGGCCAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                              >  2:1549542/1‑149 (MQ=255)
                        gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                              >  1:1336042/1‑149 (MQ=255)
                        gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                              <  1:429430/149‑1 (MQ=255)
                        gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                              >  2:476036/1‑149 (MQ=255)
                        gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACATTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                              >  2:192456/1‑149 (MQ=255)
                           gTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAtttt                                                                                                                           >  1:949934/1‑149 (MQ=255)
                           gTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAAAAACAGGGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGCACAATATCAATCCGGCGGAGGTGTAAAACACCACCGGATTGAtttt                                                                                                                           <  2:1291649/149‑1 (MQ=25)
                           gTTCAATATAAAATGTGTCTTAATTGTTACGAAGTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAAACAATCCGGTGGTGTTTTACACCTCCGACGGACTCAtttt                                                                                                                           >  2:856770/1‑149 (MQ=255)
                             tCAATATAAAATGTGTCTTAATTTTTACGAATTTGATTTTAAATAGTAACAATAGCAGTGGGGATACTGGATGTCATACCTTTCACACAGCAGGGAAAAGGAGAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGTACTGATTttgt                                                                                                                         <  1:883810/149‑1 (MQ=255)
                             ccacaataAAATGTGTCTAAATTCTGACTAAGATGATTTTAAATAGTACCACTAACAGTGGTGATACAGGATGACACACCTTTCACACAGCAGGGAAAAGGACAAAATCACTCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgt                                                                                                                         <  1:1571614/144‑1 (MQ=255)
                               aaTATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGTATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgttg                                                                                                                       >  1:742309/1‑149 (MQ=255)
                                  ataAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTAACACAGAAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGCACTGATTTTGTTGttt                                                                                                                    >  1:717965/1‑149 (MQ=255)
                                   aaacaTGTGGCTTAAATGTTACGCATTAGATATTAAATAGTAACAATAACAGTGGGGATACTGGCTGACAGACCATTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGGGTTTTACACCTCCGCCGGACTGATTTTGTTGtttt                                                                                                                   <  1:1452339/145‑1 (MQ=255)
                                    aaaaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTc                                                                                                                  >  2:217967/1‑149 (MQ=255)
                                      aaTTTGTCTTAATTGTTTCGAATTTTATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTccc                                                                                                                <  2:483101/149‑1 (MQ=17)
                                      aaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTccc                                                                                                                <  1:574767/149‑1 (MQ=255)
                                       aTGTCTCTTAATTGTTACTAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCt                                                                                                               <  1:753072/149‑1 (MQ=255)
                                        tgtgTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTCCCTg                                                                                                              <  1:1254509/149‑1 (MQ=18)
                                          tgtCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGAc                                                                                                            <  2:1476480/149‑1 (MQ=255)
                                           gtCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGCCCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGTACGGATTTTGTTGTTTTCCCTGACa                                                                                                           >  1:310775/1‑149 (MQ=255)
                                            tCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACaa                                                                                                          <  1:1032599/149‑1 (MQ=255)
                                            tCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACaa                                                                                                          <  2:638133/149‑1 (MQ=255)
                                            tCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTCCCTgctgt                                                                                                          <  1:1248881/149‑5 (MQ=17)
                                              ttAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACg                                                                                                        >  1:1348306/1‑149 (MQ=255)
                                                       aCGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATATTGTTTTTTTCCCTGACACCGATCCTGttt                                                                                               >  2:1161375/1‑149 (MQ=255)
                                                       aCGAATTTGATTTTAAATAGTAACAATAACAGTCGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGtta                                                                                               >  2:688726/1‑148 (MQ=255)
                                                        cGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTc                                                                                              >  2:355984/1‑149 (MQ=255)
                                                           aTTTGATTTTAAATAGTAACAATAACAGTTGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATGAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGAGAACGATCCTGTTTcgcg                                                                                           <  1:731204/149‑1 (MQ=255)
                                                            tttGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGa                                                                                          <  2:1256250/149‑1 (MQ=255)
                                                                   ttAAATAGTAACAATAACAGTGGGTATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCg                                                                                   <  1:975159/149‑1 (MQ=255)
                                                                     aaaTAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGCCTTCTCGGc                                                                                 >  1:530519/1‑149 (MQ=255)
                                                                      aaTAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGTAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGcc                                                                                <  1:266243/149‑1 (MQ=255)
                                                                           tAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCTGGTAAAAGGACAAAATCAATCCGGTGGTGGTTTACACCTCCGCCGTACTGATTTTGTTTTTTTCCCTGACAACGTTCCTGTTTCGCGACTTCTCGGCCCtgtg                                                                           <  1:42128/149‑1 (MQ=255)
                                                                             aCAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCGGTGGa                                                                         >  1:1191168/1‑149 (MQ=255)
                                                                               aaTAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGAtt                                                                       <  1:290370/149‑1 (MQ=255)
                                                                                 tAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTgg                                                                     <  1:949489/149‑1 (MQ=255)
                                                                                  aaCAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGc                                                                    <  2:1470012/149‑1 (MQ=255)
                                                                                  aaCAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGc                                                                    <  1:3797/149‑1 (MQ=255)
                                                                                   aCAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGGTTGGcc                                                                   <  2:1090531/149‑1 (MQ=255)
                                                                                   aCAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGcc                                                                   <  2:949934/149‑1 (MQ=255)
                                                                                      gTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCa                                                                <  1:1161375/149‑1 (MQ=255)
                                                                                           gATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTg                                                           <  2:1348306/149‑1 (MQ=255)
                                                                                            aTACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTgg                                                          >  2:408198/1‑149 (MQ=255)
                                                                                            aTACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGACTGGCCGCACGCTgg                                                          >  1:1271206/1‑149 (MQ=255)
                                                                                                 ggATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTgg                                                     >  1:1396579/1‑149 (MQ=255)
                                                                                                 ggATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAATTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTgg                                                     >  1:90806/1‑149 (MQ=255)
                                                                                                     gACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGttt                                                 <  2:671282/149‑1 (MQ=255)
                                                                                                     gACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGGGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGttt                                                 <  1:1470648/149‑1 (MQ=255)
                                                                                                        aGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTa                                              >  2:1209278/1‑149 (MQ=255)
                                                                                                          aCCTTTAACTCAGCAGGTAAAAGGACAAAATCAATCCGGTGGTGTTTTACACATCCGCCTGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTaaa                                            <  2:1388443/149‑1 (MQ=255)
                                                                                                           ccTTTCACACAGCAGTGAAAAGGAAAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGCCAAAGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGACGCAAGCTGGACTGGTTTTCTaaaa                                           >  1:289417/1‑149 (MQ=255)
                                                                                                           ccTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTaaaa                                           <  1:222940/149‑1 (MQ=255)
                                                                                                            cTTTCACACAGCAGGGAAAAGGAGAAAATCAATCCGGGGGTGTTTTACACCTCCGCCTGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAAc                                          <  1:856770/149‑1 (MQ=255)
                                                                                                            cTTTCACACAGCAGGGAAAAGGACAAAATCAATCCTGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAAc                                          <  2:1036184/149‑1 (MQ=255)
                                                                                                             tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAAcc                                         <  2:1265026/149‑1 (MQ=255)
                                                                                                             tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAAcc                                         >  1:1253147/1‑149 (MQ=255)
                                                                                                             tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAAcc                                         >  1:1343737/1‑149 (MQ=255)
                                                                                                                 acacaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCGCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCg                                     <  2:1415126/149‑1 (MQ=255)
                                                                                                                   acaGCCGGGAAAAGGCCAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGt                                   <  2:1533330/149‑1 (MQ=255)
                                                                                                                   acaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGCTTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGCCTTCTCGGCCCTGTGGATTGGCCGCACGCTGGCCTGGGTTTCTAAAACCTTCGGt                                   >  1:915084/1‑149 (MQ=255)
                                                                                                                   acaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGt                                   >  2:958988/1‑149 (MQ=255)
                                                                                                                   acaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGt                                   <  1:217967/149‑1 (MQ=255)
                                                                                                                     agcagGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTg                                 >  1:1105960/1‑149 (MQ=255)
                                                                                                                        agGGAAAAGGACAATATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTa                              <  1:1266583/149‑1 (MQ=255)
                                                                                                                        agGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTa                              <  1:62777/149‑1 (MQ=255)
                                                                                                                         gggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTAc                             <  1:224900/149‑1 (MQ=255)
                                                                                                                          ggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACt                            <  2:887859/149‑1 (MQ=255)
                                                                                                                          ggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACt                            >  1:630984/1‑149 (MQ=255)
                                                                                                                            aaaaGGACAAAATCAATCCGGTGGTGGTTTACACCTCCGCCGGACTGATTTTGGTGTTTTCCCTGTCAACGATCCTGTTTCGCGTCTTCTCGGCCCTGTGGATTGGGCGCACGCTGGACTGGGTTTCTAAAACCTTCTGTTGGTACTAt                          <  1:856762/149‑1 (MQ=255)
                                                                                                                               aGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctg                       <  2:305646/149‑1 (MQ=255)
                                                                                                                                 gACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgct                     >  1:788565/1‑149 (MQ=255)
                                                                                                                                 gACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgct                     >  1:1508700/1‑149 (MQ=255)
                                                                                                                                 gACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgct                     >  2:1289100/1‑149 (MQ=255)
                                                                                                                                 gACAAAATCAAACCGGTCGTGTTTTACACCTCCGCCTGACTGATTTTGTTGTTATCCCTGACAACGAACCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgct                     >  2:339388/1‑149 (MQ=255)
                                                                                                                                   cAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGCTCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTgg                   >  2:816784/1‑149 (MQ=255)
                                                                                                                                   cAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCACGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTGCTAAAACATTCGGTTGTTACTATCTGCTgg                   >  2:1546748/1‑149 (MQ=255)
                                                                                                                                     aaaTCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggcg                 >  2:141362/1‑149 (MQ=255)
                                                                                                                                      aatcaatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggcgg                >  2:635460/1‑149 (MQ=255)
                                                                                                                                            tcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGc          >  1:801032/1‑149 (MQ=255)
                                                                                                                                             ccGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGct         <  1:422496/149‑1 (MQ=255)
                                                                                                                                             ccGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGct         <  1:1461077/149‑1 (MQ=255)
                                                                                                                                               ggtggtGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGctct       <  1:1494677/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGATACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATT  >  NC_000913/329364‑329660

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