Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,025,296 T→G S316A (TCA→GCA)  xanQ → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,025,2960TG98.6% 256.4 / ‑5.7 71S316A (TCA→GCA) xanQxanthine permease
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (29/41);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CAATTTGCCGTTAATCACCATCCCGCATCCGTTCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTTCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAA  >  NC_000913/3025061‑3025519
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
cattttcctgTTAATCACCACCCCGCCTCCTTTCAAATCCGGTTTTAGTTTTCGCTCCCATCAGTTCCTGGTTGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCa                                                                                                                                                                                                                               <  2:243764/229‑1 (MQ=255)
   gtggctgTTAATCACCATCCCCCACCCGTCAAAATCCGTTTTTATTTTTAGTTTCCATCATTTCCTGGTGGTTGGCACGATTTACCTTATTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTAc                                                                                                                                                                                                                           <  2:892421/233‑1 (MQ=255)
             aTCACCATCCCGCATCCGTTCAAATACGGATTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCACCATTCAGGGGGAAGAGCATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCTGCTGTCGGTGCATTACCATTaacaac                                                                                                                                                                                                                 >  2:748767/1‑239 (MQ=255)
                    cccagCAGCCTTTCAAATACGGTTTTAGTTTTGGCTGCCATCAGTTCCTGGTGGTGGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTgcg                                                                                                                                                                                                          <  2:930561/235‑1 (MQ=255)
                       cGCCTCCGTTAAAAGTCGTTTTTTTTTTTAGCTTCCATCAGTCCCTGGGGTTGGCCACGATTTATCTGCTGAGGGTGCTGGAAGCGAGCGGCGATAGCGCCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCaa                                                                                                                                                                                                       <  2:613216/239‑1 (MQ=255)
                       cGCATCCGTTCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGAATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGGCTGGTTTCTGTTATCGCCTCCGCTGCCGGGTCATTACCATTAACCACGTTTgtcaaa                                                                                                                                                                                                       >  1:603071/1‑234 (MQ=255)
                         cATCCGTTCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCaaaa                                                                                                                                                                                                     >  2:373712/1‑239 (MQ=255)
                            ccGTTCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCGGCTTAGCGTGCTGGAACCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAaataa                                                                                                                                                                                                  <  1:974439/239‑1 (MQ=255)
                               ttCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATgg                                                                                                                                                                                               >  1:1029857/1‑239 (MQ=255)
                                          cTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGa                                                                                                                                                                                    >  1:1144743/1‑239 (MQ=255)
                                             tAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATg                                                                                                                                                                                  >  2:132050/1‑238 (MQ=255)
                                             tAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATg                                                                                                                                                                                  >  1:861554/1‑238 (MQ=255)
                                                ttttAGCTTCCATCAGTCCCGGGTGGTTGGCACGATTTATCTGCTTAGCGTGCGGGAAGCAGTCGGCGATACCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCa                                                                                                                                                                                                        <  2:965185/213‑1 (MQ=255)
                                                   tAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATg                                                                                                                                                                                                >  1:874016/1‑218 (MQ=255)
                                                        gccttgtgTTCGTGGTGTTTGACCCGATTTATCTGCTTAGCGTTCTGGACGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCt                                                                                                                                                                      <  2:1077618/232‑1 (MQ=255)
                                                         ccATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCg                                                                                                                                                                           <  1:284264/233‑1 (MQ=255)
                                                         ccATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCg                                                                                                                                                                           >  2:284264/1‑233 (MQ=255)
                                                          cATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATggg                                                                                                                                                                                              >  1:965185/1‑213 (MQ=255)
                                                            gcgtttCCGGGTGGTTGCCTCGTTTTATTTACTTGGCGTGCTGGAAGCAGCGGGCGATTCCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCAc                                                                                                                                                                  <  2:1144743/235‑1 (MQ=255)
                                                                     ggtggtTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCCCGTTATGTCg                                                                                                                                                         >  1:285770/1‑239 (MQ=255)
                                                                                  gATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATcg                                                                                                                                             <  1:373712/239‑1 (MQ=255)
                                                                                     ttATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCa                                                                                                                                                                                                                               >  2:917231/1‑154 (MQ=255)
                                                                                     ttATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCa                                                                                                                                                                                                                               <  1:917231/154‑1 (MQ=255)
                                                                                         tgtcTTAGCGGGCTGGAGGCAGTGGGCGATACCCCCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAAt                                                                                                                                      <  2:253922/236‑1 (MQ=255)
                                                                                                      tGGAAGCATTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATgggg                                                                                                                                                                                             <  2:382315/171‑1 (MQ=255)
                                                                                                      tGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATgggg                                                                                                                                                                                             >  1:382315/1‑171 (MQ=255)
                                                                                                              gTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGAt                                                                                                                                                                                   >  1:1002910/1‑173 (MQ=255)
                                                                                                                gggcattacccccGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCt                                                                                                                     <  2:603071/222‑1 (MQ=255)
                                                                                                                   cGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCg                                                                                                             <  1:529380/238‑1 (MQ=255)
                                                                                                                               cACGGCAATGGTTTCCGGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTggcggc                                                                                                 <  1:1116100/238‑1 (MQ=255)
                                                                                                                               cACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGct                                                                                                <  2:656185/239‑1 (MQ=255)
                                                                                                                                    cAATGGTTTCACGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGcttctt                                                                                            <  1:912172/238‑1 (MQ=255)
                                                                                                                                     aaTGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCa                                                                                          <  2:1037131/239‑1 (MQ=255)
                                                                                                                                                   gccCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGGTATTCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCg                                                                            >  1:878335/1‑239 (MQ=255)
                                                                                                                                                     cccATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCgg                                                                           <  1:911428/238‑1 (MQ=255)
                                                                                                                                                                 gAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAg                                                               <  2:1029857/238‑1 (MQ=255)
                                                                                                                                                                   agagTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAAc                                                                                                                                                                                                                    >  1:883494/1‑87 (MQ=255)
                                                                                                                                                                   agagTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAAc                                                                                                                                                                                                                    <  2:883494/87‑1 (MQ=255)
                                                                                                                                                                      gTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTggcggc                                                                                                 <  1:950815/199‑1 (MQ=255)
                                                                                                                                                                      gTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTggcggc                                                                                                 >  2:950815/1‑199 (MQ=255)
                                                                                                                                                                        aTCCGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCaa                                                         <  1:456432/237‑1 (MQ=255)
                                                                                                                                                                        aTCAGTCCCGGCTGTAAGGGGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGTCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCGCGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCaa                                                         <  2:511995/237‑1 (MQ=255)
                                                                                                                                                                          cAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAaataa                                                                                                                                                                                                  >  2:1053440/1‑98 (MQ=255)
                                                                                                                                                                          cAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAaataa                                                                                                                                                                                                  <  1:1053440/98‑1 (MQ=255)
                                                                                                                                                                              cccGGCTGAAAGGCGGCGTGCGGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTg                                                 <  2:315100/239‑1 (MQ=255)
                                                                                                                                                                               ccGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCtt                                                                                                                    <  2:523779/171‑1 (MQ=255)
                                                                                                                                                                               ccGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCtt                                                                                                                    >  1:523779/1‑171 (MQ=255)
                                                                                                                                                                               ccGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTg                                                 <  1:132050/238‑1 (MQ=255)
                                                                                                                                                                               ccGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGa                                                <  1:681304/239‑1 (MQ=255)
                                                                                                                                                                                  gCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCa                                                                                                                                                                                                        <  2:365333/84‑1 (MQ=255)
                                                                                                                                                                                  gCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCa                                                                                                                                                                                                        >  1:365333/1‑84 (MQ=255)
                                                                                                                                                                                             gcgTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATTATTGGCGGCTTCTTCACGACCATTCCCTCGGCCGTTCTGGGAGGCGCCATGACGTTGGTGTTTTCCATGAtt                                  >  1:939809/1‑239 (MQ=255)
                                                                                                                                                                                                tGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAAGCGCAATGACGGTGATGTTTTTCATGATTTcc                               >  1:855294/1‑239 (MQ=255)
                                                                                                                                                                                                     gcagtgtgtttGGTTTCTTTTATCGCCTCCGCTGTTGGTGCATTACCATTACCCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTGCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGc                          <  1:556235/229‑1 (MQ=255)
                                                                                                                                                                                                     gCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAAt                                                        <  1:710562/209‑1 (MQ=255)
                                                                                                                                                                                                          tGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCAATTCCCCGGCAGGTCCGGGGGGGGCAATGACCTTGATGTTTTCCATGATTGCCATCGCCgggg                     >  2:478504/1‑238 (MQ=255)
                                                                                                                                                                                                                 gTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATg                                                                                                                                                                                  <  2:832048/75‑1 (MQ=255)
                                                                                                                                                                                                                 gTTTCTGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATg                                                                                                                                                                                  >  1:832048/1‑75 (MQ=255)
                                                                                                                                                                                                                      tGTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGACGTTTCACGTTATGTCGGGCGGACCATTGCGGTAATGCTGGGTATCCTCGGGTTATTTTCGCAGATTGGCGGCTTCTTTCCGACAATTCCCTCGGCAGATTTGGGAGGCGGCATGGCGTTGGTGGTTTTCCCGATTGCCCTCGCAGGGAGTCGcactatat         >  2:554083/1‑233 (MQ=255)
                                                                                                                                                                                                                       gTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATgggg                                                                                                                                                                                             <  2:1032122/58‑1 (MQ=255)
                                                                                                                                                                                                                       gTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATgggg                                                                                                                                                                                             >  1:1032122/1‑58 (MQ=255)
                                                                                                                                                                                                                       gTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTTCCTCGGCAGTTCTGGGAg                                                               >  2:951499/1‑184 (MQ=255)
                                                                                                                                                                                                                       gTTATCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAg                                                               <  1:951499/184‑1 (MQ=255)
                                                                                                                                                                                                                           tCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCt                                                                                                                          <  1:110344/121‑1 (MQ=255)
                                                                                                                                                                                                                           tCGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCt                                                                                                                          >  2:110344/1‑121 (MQ=255)
                                                                                                                                                                                                                            cGCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCa   <  2:824676/239‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCaa  <  2:346768/239‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCaa  <  2:1165309/239‑1 (MQ=255)
                                                                                                                                                                                                                                tCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTggcggc                                                                                                 >  1:201600/1‑141 (MQ=255)
                                                                                                                                                                                                                                tCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTggcggc                                                                                                 <  2:201600/141‑1 (MQ=255)
                                                                                                                                                                                                                                tCCGCTGTCGGTGCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCCCGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCa   >  2:359389/1‑235 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
CAATTTGCCGTTAATCACCATCCCGCATCCGTTCAAATACGGCTTTAGTTTTAGCTTCCATCAGTTCCTGGTGGTTGGCACGATTTATCTGCTTAGCGTGCTGGAAGCAGTCGGCGATATCACCGCCACGGCAATGGTTTCCCGCCGCCCCATTCAGGGGGAAGAGTATCAGTCCCGGCTGAAAGGCGGCGTGCTGGCAGATGGTCTGGTTTCTGTTATCGCCTCCGCTGTCGGTTCATTACCATTAACCACGTTTGCGCAAAATAATGGGGTTATTCAGATGACTGGCGTCGCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAA  >  NC_000913/3025061‑3025519

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