Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 331,406 C→A 7.5% T648T (ACC→ACA betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913331,4060CA 7.5% 333.3 / 9.0 100T648T (ACC→ACAbetTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (44/47);  new base A (2/6);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTC  >  NC_000913/331269‑331542
                                                                                                                                         |                                                                                                                                        
ccGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  2:1660789/140‑1 (MQ=255)
ccGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  2:1413420/140‑1 (MQ=255)
ccGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  1:1192748/140‑1 (MQ=255)
ccGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  1:1408645/140‑1 (MQ=255)
     gcaATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTc                                                                                                                                   <  1:115831/140‑1 (MQ=255)
          atTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAc                                                                                                                              <  1:652596/140‑1 (MQ=255)
          atTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAc                                                                                                                              <  2:1256376/140‑1 (MQ=255)
          atTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAc                                                                                                                              <  1:1517846/140‑1 (MQ=255)
           ttCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAcc                                                                                                                             <  1:1823662/140‑1 (MQ=255)
              ggTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                                                                          >  2:367121/1‑140 (MQ=255)
              ggTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt                                                                                                                          <  1:367121/140‑1 (MQ=255)
                tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACGGCAGAGAGCAGGTGATCACCGATAGTCGTGACCAGTAc                                                                                                                        >  1:447477/1‑140 (MQ=255)
                tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTAc                                                                                                                        >  1:1340541/1‑140 (MQ=255)
                tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTAc                                                                                                                        >  1:489684/1‑140 (MQ=255)
                tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTAc                                                                                                                        >  1:1110211/1‑140 (MQ=255)
                tGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGAAc                                                                                                                        >  1:1019772/1‑140 (MQ=255)
                     ggCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCg                                                                                                                   >  1:1151731/1‑140 (MQ=255)
                          tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGTAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCAcc                                                                                                              >  2:847019/1‑140 (MQ=255)
                          tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCAcc                                                                                                              >  2:700223/1‑140 (MQ=255)
                          tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCAcc                                                                                                              >  1:1716967/1‑140 (MQ=255)
                          tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCAcc                                                                                                              >  2:846576/1‑140 (MQ=255)
                          tacctaccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCATGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCAcc                                                                                                              >  1:794821/1‑140 (MQ=255)
                              taccGAGCAAGCAGCGGTTAATCGACCTTCTACCGGCTGGAAACCTTCCTGTTAGAATGCTGCCAGGGCATCGACCTGATGGACTACAGCTATGTGCAGGTGATCTCCGATATTCTTGACCAGTACGAGGGGCTCCTTaa                                                                                                          <  2:522926/140‑1 (MQ=255)
                               accGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCGGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAAc                                                                                                         >  1:352521/1‑140 (MQ=255)
                                ccGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACt                                                                                                        >  2:980753/1‑140 (MQ=255)
                                 cgcgCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACtt                                                                                                       >  1:1523009/1‑140 (MQ=255)
                                   cgcacgcaGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTa                                                                                                     <  2:700142/140‑1 (MQ=255)
                                             ggTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCa                                                                                           <  2:1649835/140‑1 (MQ=255)
                                                      aCCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGTCCAGTACGAGCGGCACCTTAACTTTATTCATCTCCAACGTGAAgc                                                                                  >  2:817706/1‑140 (MQ=255)
                                                      aCCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAgc                                                                                  >  2:1336172/1‑140 (MQ=255)
                                                      aCCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAgc                                                                                  >  1:1470897/1‑140 (MQ=255)
                                                        ctactaCCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAgcgc                                                                                <  2:1523009/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGCAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCTCCGATATTCTTGACCTGTACGAGCGGCACCTTATCTTTTTTCATCTCCATCGTGAAGCGcc                                                                               <  2:114363/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGcc                                                                               <  1:1762122/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGcc                                                                               <  1:1761684/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGcc                                                                               <  1:1520979/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGcc                                                                               <  1:1520790/140‑1 (MQ=255)
                                                         tactacCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGcc                                                                               <  1:1520537/140‑1 (MQ=255)
                                                              ccGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGcc                                                                          <  2:1244349/140‑1 (MQ=255)
                                                              ccGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGcc                                                                          <  2:1200502/140‑1 (MQ=255)
                                                               cGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCa                                                                         <  1:732179/140‑1 (MQ=255)
                                                                 gCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATa                                                                       >  2:842097/1‑140 (MQ=255)
                                                                 gCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATa                                                                       >  2:841673/1‑140 (MQ=255)
                                                                     gAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGt                                                                   >  2:382588/1‑140 (MQ=255)
                                                                     gAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGt                                                                   >  2:302516/1‑140 (MQ=255)
                                                                            tCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        >  2:221290/1‑64 (MQ=255)
                                                                            tCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  1:221290/64‑1 (MQ=255)
                                                                             ccTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTcc                                                           >  2:1647714/1‑140 (MQ=255)
                                                                                gTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCgg                                                        >  2:1680431/1‑140 (MQ=255)
                                                                                    gAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGAcgc                                                    >  1:456343/1‑140 (MQ=255)
                                                                                         cAGCCAGTGCAACGTCCTGATGTACTACTGAAAAGTGCAGGTGCTCTCCGATAATCTTGACCAGTTCGAGCGGCACCTTAACTTTATTCTTCTCCATCGTGAAGCGCAGGGCCATAGCGTGATGTTCCCGGACGCGTGAt                                               <  2:1590110/140‑1 (MQ=255)
                                                                                           gCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        >  1:340170/1‑49 (MQ=255)
                                                                                           gCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                                                                                                                                        <  2:340170/49‑1 (MQ=255)
                                                                                           gCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTg                                             >  1:50829/1‑140 (MQ=255)
                                                                                           gCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTg                                             >  1:51081/1‑140 (MQ=255)
                                                                                            ccAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGt                                            >  1:32771/1‑140 (MQ=255)
                                                                                              aGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  2:1796117/140‑1 (MQ=255)
                                                                                              aGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  1:1559666/140‑1 (MQ=255)
                                                                                              aGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  2:1855758/140‑1 (MQ=255)
                                                                                              aGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  2:1039484/140‑1 (MQ=255)
                                                                                              aGGGCAACGACCTGATGGACAACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  2:1854774/140‑1 (MQ=255)
                                                                                                ggCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTAtt                                        <  1:1784370/140‑1 (MQ=255)
                                                                                                ggCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCGTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTAtt                                        <  1:1784554/140‑1 (MQ=255)
                                                                                                     cGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATa                                   >  2:717912/1‑140 (MQ=255)
                                                                                                         cTGATGGACGACAGCAAAGAGCAGGTGATCACGGATATTCTTGACCAGTCCGAGCGGCACCTTAACATTATTCATCTCCATCGAGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATACAAc                               <  2:1500504/140‑1 (MQ=255)
                                                                                                          tGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAAcc                              >  1:364258/1‑140 (MQ=255)
                                                                                                          tGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAAcc                              >  2:315934/1‑140 (MQ=255)
                                                                                                              ggtcTACAGCTTAGTTCATGTGATCATAGATATTCTTGACCAGTACGAGCGGCTCCTTAACTTTATTCATCTACATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCTGACGCGTGATTGTTATTGCATAAAAgcgcgc                          <  2:1044719/137‑1 (MQ=255)
                                                                                                              ggACTACAGCAAAGAGCAGGTGTTCACCGATATTCTTGTCCAGTAGGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGGCGGGCCTTAGCGTGATGTTCCCGGACGGGTGATTGTTATTGCATAAAACCgcgc                          <  2:1425577/140‑1 (MQ=255)
                                                                                                              ggACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTAATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCgcgc                          <  1:846576/140‑1 (MQ=255)
                                                                                                              ggACTACAGCAAAGAGCAGGTGATCACAGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTAATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCgcgc                          <  1:847019/140‑1 (MQ=255)
                                                                                                               gACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGcc                         <  1:588703/140‑1 (MQ=255)
                                                                                                               gACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGcc                         <  1:111587/140‑1 (MQ=255)
                                                                                                               gACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCAGGACGCGTGATTGTTATTGCATAAAACCGCGcc                         <  2:676355/140‑1 (MQ=255)
                                                                                                                aatACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCa                        <  2:277755/138‑1 (MQ=255)
                                                                                                                aCTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGCCGCGTGATTGTTATTGCATAAAACCGCGCCa                        <  2:230544/140‑1 (MQ=255)
                                                                                                                aCTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCa                        <  2:1113562/140‑1 (MQ=255)
                                                                                                                aCTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCa                        <  1:633734/140‑1 (MQ=255)
                                                                                                                aCTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCa                        <  2:698989/140‑1 (MQ=255)
                                                                                                                  tACAGCAAATAGAAGTTGATATCCGTTATTCTTTTCCAGTTCGAGCTGCTCCTTAACTTTATTCATCTAAATCGTGAAGCGCAGGGCCATAGGGTGATGTTACCGGTCGCGTGTTTTTTATTGCTTAAAACCGCGCCATg                      <  2:130121/140‑1 (MQ=255)
                                                                                                                  tACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATg                      <  1:1533178/140‑1 (MQ=255)
                                                                                                                  tACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATg                      <  1:1042179/140‑1 (MQ=255)
                                                                                                                  tACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCATGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATg                      >  2:584530/1‑140 (MQ=255)
                                                                                                                  tACAGCAAAGAGCAGGTGATCACCGATAGTCTTGACCAGTACGAGCGGCACCTTAACGTTAGTCATCTCCATCGTGAAGCGCCGGGCCAAATCGTGTGGTGCCCGGACGCGGGAGGGTTATTGCTGAAAGCGGCGCCGTg                      >  2:958575/1‑140 (MQ=255)
                                                                                                                   aCAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGt                     <  1:1922485/140‑1 (MQ=255)
                                                                                                                    cAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          >  1:662547/1‑118 (MQ=255)
                                                                                                                    cAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTa                                          <  2:662547/118‑1 (MQ=255)
                                                                                                                      gCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAgcg                                                                                 >  2:760352/1‑77 (MQ=255)
                                                                                                                      gCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAgcg                                                                                 <  1:760352/77‑1 (MQ=255)
                                                                                                                      gCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTg                  >  1:520099/1‑140 (MQ=255)
                                                                                                                         aaGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCat               >  2:1378347/1‑140 (MQ=255)
                                                                                                                         aaGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGCGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCat               >  2:1379879/1‑140 (MQ=255)
                                                                                                                          agagCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTAATCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTGATTTCATAAAACCGCGCCATGTCTGCata              >  2:1913149/1‑140 (MQ=255)
                                                                                                                               aGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGgcg         <  1:1081107/140‑1 (MQ=255)
                                                                                                                                ggTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGgcgc        >  2:606047/1‑140 (MQ=255)
                                                                                                                                 gTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCt       <  2:773996/140‑1 (MQ=255)
                                                                                                                                      cACCGATATTCTTGACCAGTACGAGCGGCACCTTTACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCTTGTCTGCATATGGCGCTGtttg  >  2:1407397/1‑139 (MQ=255)
                                                                                                                                      cACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTtc  >  2:1294864/1‑140 (MQ=255)
                                                                                                                                      cACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTtc  >  2:207459/1‑140 (MQ=255)
                                                                                                                                      cACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTtc  >  2:206904/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCGGCACCTTAACTTTATTCATCTCCATCGTGAAGCGCCGGGCCATAGCGTGATGTTCCCGGACGCGTGATTGTTATTGCATAAAACCGCGCCATGTCTGCATATGGCGCTGTTTC  >  NC_000913/331269‑331542

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