Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,523,798 T→G 7.4% D63A (GAC→GCC)  xapB ← xanthosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,523,7980TG 7.4% 255.4 / 7.4 91D63A (GAC→GCC) xapBxanthosine transporter
Reads supporting (aligned to +/- strand):  ref base T (45/39);  new base G (5/2);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTGC  >  NC_000913/2523649‑2523947
                                                                                                                                                     |                                                                                                                                                     
aTAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGt                                                                                                                                                       <  2:856174/150‑1 (MQ=255)
   aCGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTg                                                                                                                                                    <  2:1021248/150‑1 (MQ=255)
   aCGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTg                                                                                                                                                    >  2:801972/1‑150 (MQ=255)
    cGCAATAGTCGGCATAAATGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGc                                                                                                                                                   >  1:29168/1‑150 (MQ=255)
     gCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCg                                                                                                                                                  <  1:1071189/150‑1 (MQ=255)
     gCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCACACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCg                                                                                                                                                  <  2:1140437/150‑1 (MQ=255)
      cAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCga                                                                                                                                                 <  1:801972/150‑1 (MQ=255)
        aTAGTCGGCATAAACGCCATCGCAATGACTACCATCACCCAAAACATCATATCCGGATCAGTTACGGGTGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCCCGTACTGCGCGCAGCCATTTGTCTGCgatg                                                                                                                                               >  1:450703/1‑150 (MQ=255)
         tAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGACAAATGGCTGCgcgca                                                                                                                                              <  2:834010/150‑5 (MQ=255)
           gTCGGCATAAACGCCATTGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACTGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATc                                                                                                                                            <  1:167131/150‑1 (MQ=255)
           gTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATc                                                                                                                                            >  1:608828/1‑150 (MQ=255)
           gTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATc                                                                                                                                            >  1:1402875/1‑150 (MQ=255)
              ggCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcccc                                                                                                                                         <  2:41586/150‑1 (MQ=255)
              ggCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcccc                                                                                                                                         <  2:419670/150‑1 (MQ=255)
               gCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTCCGCCCGCACACCCCAGGTGCCACCGCATGTATGCAAGTTCTGCGCGCAGCCATTTTTCTGCGCTGATcaccc                                                                                                                                        >  1:321122/1‑149 (MQ=255)
                       gCCATCGCATTGACTAACCTCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCa                                                                                                                                <  2:921354/150‑1 (MQ=255)
                        ccATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTCCGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAg                                                                                                                               >  2:637630/1‑150 (MQ=255)
                        ccATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAg                                                                                                                               <  1:1207797/150‑1 (MQ=255)
                         cATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAgg                                                                                                                              >  2:1697016/1‑150 (MQ=255)
                             gCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCata                                                                                                                          >  2:1711458/1‑150 (MQ=255)
                              cATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCataa                                                                                                                         >  2:1396869/1‑150 (MQ=255)
                                ttGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCAGGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCataata                                                                                                                       <  2:775923/150‑1 (MQ=255)
                                  gACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATaataat                                                                                                                     >  2:960268/1‑150 (MQ=255)
                                     tAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcgc                                                                                                                  <  1:1575405/150‑1 (MQ=255)
                                     tAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcgc                                                                                                                  <  2:1402875/150‑1 (MQ=255)
                                     tAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcgc                                                                                                                  <  1:1218760/150‑1 (MQ=255)
                                      aaCATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcgcg                                                                                                                 >  2:342910/1‑150 (MQ=255)
                                      aaCATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcgcg                                                                                                                 <  2:1321342/150‑1 (MQ=255)
                                       aCATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCgcgg                                                                                                                <  2:629485/150‑1 (MQ=255)
                                            aCCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACGAGGCATAATAATCGCGGCGATc                                                                                                           >  2:1594466/1‑150 (MQ=255)
                                              ccAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATccc                                                                                                         >  2:1363559/1‑150 (MQ=255)
                                                aaaaCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCtt                                                                                                       <  2:1096382/150‑1 (MQ=255)
                                                  aaCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTg                                                                                                     <  2:1255167/150‑1 (MQ=255)
                                                   aCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTgg                                                                                                    >  1:1586375/1‑150 (MQ=255)
                                                   aCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGAGCCCTTTgg                                                                                                    >  2:314098/1‑150 (MQ=255)
                                                   aCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTGCTGCGCGCAGCCATTTGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTgg                                                                                                    <  2:216939/150‑1 (MQ=255)
                                                    catcatATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGCTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGa                                                                                                   >  1:1321670/1‑150 (MQ=255)
                                                     atcatATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGaa                                                                                                  <  1:191183/150‑1 (MQ=255)
                                                        atatCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTg                                                                                               <  2:970792/150‑1 (MQ=255)
                                                        atatCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTg                                                                                               <  2:320866/150‑1 (MQ=255)
                                                         tatCCGTAGCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACCGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGt                                                                                              <  1:87192/150‑1 (MQ=255)
                                                          atCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGTTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTa                                                                                             <  2:906663/150‑1 (MQ=255)
                                                             gggATCAGTTACGGATGACGCATAAAAAAGTACGCCCGCACACACCATGTGACACAGCATGTATGCACGTTCTGGGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAAc                                                                                          <  2:1686499/149‑1 (MQ=255)
                                                               gATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCa                                                                                        >  2:1688478/1‑150 (MQ=255)
                                                                     ttACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCCGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCaa                                                                                  >  2:1298111/1‑150 (MQ=255)
                                                                     ttACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCaa                                                                                  >  2:1353708/1‑150 (MQ=255)
                                                                      tACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTAAGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGAACCCTTTGGAACTGTAAACCATGCCAAc                                                                                 >  2:326590/1‑150 (MQ=255)
                                                                        cGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACAt                                                                               <  1:466668/150‑1 (MQ=255)
                                                                         ggctGCCGCATAAAAAAGTACGCCCGCCCACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCCTTTGTCTGCGATGATCCCCATTCTACCAGGCCTAATAATCGCGGCGATCCCTTTGCAACTGTCACCCCTGCCATCAtt                                                                              >  2:1506873/4‑150 (MQ=255)
                                                                         ggATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACAtt                                                                              >  2:797588/1‑150 (MQ=255)
                                                                          gATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTa                                                                             >  2:643249/1‑150 (MQ=255)
                                                                           aTGCAGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAg                                                                            <  1:1178951/150‑1 (MQ=255)
                                                                              ccGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAgcgc                                                                         >  1:98272/1‑150 (MQ=255)
                                                                              ccGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAgcgc                                                                         >  1:444125/1‑150 (MQ=255)
                                                                                ggaTAAAAAAGTACGCCGGCACACACCAGGTGACACAGCATGTATGCACTTTCTGAGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCGTTTGGAACTGTATACCATGCCAACATTAGCGCCg                                                                       <  1:912127/148‑1 (MQ=255)
                                                                                gCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCg                                                                       <  2:98272/150‑1 (MQ=255)
                                                                                  aTAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGt                                                                     >  2:1698199/1‑150 (MQ=255)
                                                                                     aaaaaGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGaa                                                                  <  1:1396869/150‑1 (MQ=255)
                                                                                       aacgTACGCCCGCACACACCAGGTGACACAGCTTGTAGGCACGTTCTGCGCGCAGCCATTTGTCTGCTATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAAt                                                                <  2:1511259/147‑1 (MQ=255)
                                                                                       aaaGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAAt                                                                <  2:325611/150‑1 (MQ=255)
                                                                                          gTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGaa                                                             >  2:1379902/1‑150 (MQ=255)
                                                                                            aCGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCAATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAaga                                                           >  1:346766/1‑150 (MQ=255)
                                                                                              gCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGt                                                         >  1:844205/1‑150 (MQ=255)
                                                                                              acccGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGt                                                         >  1:1712278/2‑150 (MQ=255)
                                                                                                ccGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCGGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTAt                                                       <  2:1685041/150‑1 (MQ=255)
                                                                                                  gCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTa                                                     >  1:14207/1‑150 (MQ=255)
                                                                                                         ccAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGTCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGt                                              <  1:342910/150‑1 (MQ=255)
                                                                                                         ccAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGAGATCCCTTTGGAACTGTAAACCATGCCAAAATTAGCGCCGGTGAAATGAAGAGTATTAATCATGt                                              >  1:1550385/1‑150 (MQ=255)
                                                                                                            ggTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATAAGCGCCGGTGAAATGAAGAGTATTAATCATGTAag                                           >  2:342164/1‑150 (MQ=255)
                                                                                                                 cacaGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGATGAGTATTAATCATGTAAGAGCCg                                      >  1:862925/1‑150 (MQ=255)
                                                                                                                 cacaGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCg                                      >  2:523552/1‑150 (MQ=255)
                                                                                                                  acaGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGa                                     >  1:413754/1‑150 (MQ=255)
                                                                                                                   ctgtctttATTCAAGTTCTGAGCGCATCAATTTGTCTGCGATTATCCAAATTATACAAGGCATAATAATAGCGGAGATCCCTTTTGAACTGTAAACCATGCAAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAg                                    <  2:1060513/143‑1 (MQ=255)
                                                                                                                   caGCATGTATGCACGTACTGCGCGCACCCATTTGTCTGCGATGATGCCCATTATAGCAGGGTAATTAATCGCGGCGATCCCTTTGGATCTGTTTTCCATGCCAACATTAGCGCCGGTGAAATGAAGCGTTTTATTCATGTAAGAGCCGAg                                    <  2:1066038/150‑1 (MQ=255)
                                                                                                                    aGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAgg                                   <  2:1182061/150‑1 (MQ=255)
                                                                                                                      cATGTATGCACGTTCTGCGCGCAACCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGt                                 >  2:1503399/1‑150 (MQ=255)
                                                                                                                      cATGTAAGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGt                                 <  1:326590/150‑1 (MQ=255)
                                                                                                                          tATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAcc                             >  1:1027230/1‑150 (MQ=255)
                                                                                                                           aTGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCGGGTGAAATGAAGAGTATTTATCATGTAAGAGCCGAGGGTATcca                            >  2:1143409/1‑150 (MQ=255)
                                                                                                                              cACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAccagcc                         <  2:29168/150‑1 (MQ=255)
                                                                                                                                   gatcatcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTc                    >  2:1544677/7‑150 (MQ=255)
                                                                                                                                   gatcatcgcAGACAAATGGCTGCGATGATCCCCATTAAACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCTAGGGTAACCACCCAGCTc                    >  2:305643/7‑150 (MQ=255)
                                                                                                                                    cTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTcc                   >  2:904914/1‑150 (MQ=255)
                                                                                                                                     tGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccc                  >  1:589512/1‑150 (MQ=255)
                                                                                                                                      catcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTcccc                 >  1:605042/4‑150 (MQ=255)
                                                                                                                                       cgcgcAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccccc                >  1:935336/1‑150 (MQ=255)
                                                                                                                                       atcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccccc                >  2:799035/3‑150 (MQ=255)
                                                                                                                                       atcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCAGGTGAAATGAAGATTATTATTAATGTAAGATCCGAGGGTAACCAGCCAGCTcccac                >  1:1487278/3‑148 (MQ=255)
                                                                                                                                         cgcAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAg              >  2:431109/1‑150 (MQ=255)
                                                                                                                                         cgcAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAg              <  2:1027230/150‑1 (MQ=255)
                                                                                                                                             gCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATATTCGCGGCGAGCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATaa          >  2:1108148/1‑150 (MQ=255)
                                                                                                                                                    gTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTg   <  1:637630/150‑1 (MQ=255)
                                                                                                                                                     tCTGCGATGATCCCCATTATACCATGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTGc  <  1:1711458/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTGC  >  NC_000913/2523649‑2523947

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