Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,523,798 T→G 8.3% D63A (GAC→GCC)  xapB ← xanthosine:H(+) symporter XapB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,523,7980TG 8.3% 227.9 / 8.3 81D63A (GAC→GCC) xapBxanthosine:H(+) symporter XapB
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base G (5/2);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTG  >  NC_000913/2523650‑2523946
                                                                                                                                                    |                                                                                                                                                    
tAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTc                                                                                                                                                     <  1:161043/150‑1 (MQ=255)
tAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTc                                                                                                                                                     <  2:630121/150‑1 (MQ=255)
tAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTc                                                                                                                                                     <  2:1224666/150‑1 (MQ=255)
 aaCGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCt                                                                                                                                                    >  2:105714/1‑150 (MQ=255)
  aCGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTATg                                                                                                                                                   >  1:669359/1‑150 (MQ=255)
   cGCAATAGTCGGCGTAAACGCCATCGCATTGACTAACATCACCCAAAACGTCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGc                                                                                                                                                  >  2:964638/1‑150 (MQ=255)
           tCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcc                                                                                                                                          <  2:1262470/150‑1 (MQ=255)
           tCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcc                                                                                                                                          <  2:43395/150‑1 (MQ=255)
           tCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcc                                                                                                                                          <  1:581975/150‑1 (MQ=255)
             ggCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATcccc                                                                                                                                        <  2:717431/150‑1 (MQ=255)
                 tAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTATGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATta                                                                                                                                    >  1:451751/1‑150 (MQ=255)
                   aaCGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATtata                                                                                                                                  >  1:244440/1‑150 (MQ=255)
                                   cTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcg                                                                                                                  >  1:336751/1‑150 (MQ=255)
                                   cTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCACGCACACACCAGGTGACACAGCATGTATGCACGTTCTGAGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATcg                                                                                                                  >  2:1473934/1‑150 (MQ=255)
                                      aCATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCgcgg                                                                                                               <  1:738431/150‑1 (MQ=255)
                                       cATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACAAGGCATAATAATCgcggc                                                                                                              >  2:667333/1‑150 (MQ=255)
                                         tCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGa                                                                                                            >  2:1163640/1‑150 (MQ=255)
                                            cccAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATcc                                                                                                         >  1:810268/1‑150 (MQ=255)
                                            cccAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTAGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATcc                                                                                                         >  1:281750/1‑150 (MQ=255)
                                             ccAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATccc                                                                                                        >  2:1396371/1‑150 (MQ=255)
                                                aaaCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCttt                                                                                                     >  1:1136824/1‑150 (MQ=255)
                                                  aCATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTgg                                                                                                   <  1:379157/150‑1 (MQ=255)
                                                   catcatATCCGGATCAGTTACGGATGCGGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGa                                                                                                  >  2:125130/1‑150 (MQ=255)
                                                       atatCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTg                                                                                              <  2:658453/150‑1 (MQ=255)
                                                            cGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAAc                                                                                         <  1:350801/150‑1 (MQ=255)
                                                                  aGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGcc                                                                                   <  1:1538028/150‑1 (MQ=255)
                                                                   gTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACGCAGCATGTATGCACGTTATGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCa                                                                                  >  2:945081/1‑150 (MQ=255)
                                                                   gTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCa                                                                                  <  2:810268/150‑1 (MQ=255)
                                                                        ggATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACAtt                                                                             >  2:1655843/1‑150 (MQ=255)
                                                                          aTGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATTCCAACATTAg                                                                           <  2:62802/150‑1 (MQ=255)
                                                                            gCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAgcg                                                                         >  1:1173803/1‑150 (MQ=255)
                                                                            gCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAgcg                                                                         >  1:1325405/1‑150 (MQ=255)
                                                                             ccGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAgcgc                                                                        >  1:1048593/1‑150 (MQ=255)
                                                                              cGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGcc                                                                       >  1:866474/1‑150 (MQ=255)
                                                                              cGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGcc                                                                       >  2:361448/1‑150 (MQ=255)
                                                                                  tAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTg                                                                   <  1:1655843/150‑1 (MQ=255)
                                                                                  tAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTg                                                                   <  2:1508984/150‑1 (MQ=255)
                                                                                    aaaaaGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGaa                                                                 >  2:40539/1‑150 (MQ=255)
                                                                                    aaaaaGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGaa                                                                 >  2:20339/1‑150 (MQ=255)
                                                                                      aaaGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAAt                                                               <  2:336751/150‑1 (MQ=255)
                                                                                          tACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAag                                                           <  1:509618/150‑1 (MQ=255)
                                                                                            cGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAagag                                                         >  2:328064/1‑150 (MQ=255)
                                                                                              cccGCACACACCAGGTGCCCCAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGCAATGAAGAGtt                                                       >  2:604751/1‑149 (MQ=255)
                                                                                                cGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTAtt                                                     <  2:451751/150‑1 (MQ=255)
                                                                                                 gCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTa                                                    <  2:244440/150‑1 (MQ=255)
                                                                                                    cacacCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATc                                                 <  1:420722/150‑1 (MQ=255)
                                                                                                    cacacCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATc                                                 <  2:1173803/150‑1 (MQ=255)
                                                                                                      cacCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCAt                                               >  1:1644319/1‑150 (MQ=255)
                                                                                                        ccAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGt                                             >  1:1672225/1‑150 (MQ=255)
                                                                                                        ccAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGt                                             >  1:343740/1‑150 (MQ=255)
                                                                                                        acAGGTGACATATCATGTATTCACGTTTTGCGCGGAGCTTTTTGTATTCGATGATCCCCATTATACCAGGCTTAATAATCGCGGGGATCCCTTTGGAAGTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTATTCATGt                                             <  2:866474/149‑1 (MQ=255)
                                                                                                          aGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTaa                                           >  1:1641094/1‑150 (MQ=255)
                                                                                                                  caGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAg                                   <  1:681345/150‑1 (MQ=255)
                                                                                                                  caGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAg                                   >  1:1236874/1‑150 (MQ=255)
                                                                                                                    gCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAggg                                 <  2:1081256/150‑1 (MQ=255)
                                                                                                                    gCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAggg                                 <  2:1644319/150‑1 (MQ=255)
                                                                                                                     taTGTATGCACGTTCTGCGCGCAGCCATTTGTTTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGt                                <  2:669359/149‑1 (MQ=255)
                                                                                                                     cATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGt                                >  2:1266509/1‑150 (MQ=255)
                                                                                                                     cATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGt                                <  2:1136824/150‑1 (MQ=255)
                                                                                                                      aTGTATGCACGTTCTGCGCGCAGCCATTTGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTa                               <  2:520866/150‑1 (MQ=255)
                                                                                                                           tGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAccag                          >  1:1382868/1‑150 (MQ=255)
                                                                                                                             cACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAccagcc                        <  1:964638/150‑1 (MQ=255)
                                                                                                                                   cTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAccagcc                        >  2:44518/1‑144 (MQ=255)
                                                                                                                                   cTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAAccagcc                        <  1:44518/144‑1 (MQ=255)
                                                                                                                                   cTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTcc                  >  1:1403470/1‑150 (MQ=255)
                                                                                                                                      atcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccccc               <  2:313126/148‑1 (MQ=255)
                                                                                                                                      atcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccccc               >  1:801961/3‑150 (MQ=255)
                                                                                                                                      atcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTccccc               >  1:313126/3‑150 (MQ=255)
                                                                                                                                       tcgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGACGAGGGTAACCAGCCAGCTCCCCCa              >  2:827511/2‑150 (MQ=255)
                                                                                                                                       gcgcAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACAAGCCAGCTCCCCCa              <  1:1659124/150‑1 (MQ=255)
                                                                                                                                        cgcAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAg             <  2:1325405/150‑1 (MQ=255)
                                                                                                                                        cgcAGACAAATGGCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAg             >  2:1356184/1‑150 (MQ=255)
                                                                                                                                         gcAGCCATTTGTCTGCGTTGATCCCCATTATACCAGGCATATTAATCGCGGCGATCCCTTTGGAACTTTAAACCATGCCAACATTAGCGCCGGTGAAATGAATAGTATTAATCATGTAATAGCCGAGGGTAACCAGCCAGCTCCCCCAGa            <  2:469622/150‑1 (MQ=255)
                                                                                                                                          cAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTa                               <  1:1045935/130‑1 (MQ=255)
                                                                                                                                          cAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTa                               >  2:1045935/1‑130 (MQ=255)
                                                                                                                                              cATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGTCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATaaaa       >  2:477592/1‑150 (MQ=255)
                                                                                                                                              cATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATaaaa       >  1:708461/1‑150 (MQ=255)
                                                                                                                                               aTTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATg                                              <  1:649232/110‑1 (MQ=255)
                                                                                                                                               aTTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATg                                              >  2:649232/1‑110 (MQ=255)
                                                                                                                                                  tGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATAtt   >  2:685435/1‑150 (MQ=255)
                                                                                                                                                   gTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTg  <  2:1236874/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAACGCAATAGTCGGCATAAACGCCATCGCATTGACTAACATCACCCAAAACATCATATCCGGATCAGTTACGGATGCCGCATAAAAAAGTACGCCCGCACACACCAGGTGACACAGCATGTATGCACGTTCTGCGCGCAGCCATTTGTCTGCGATGATCCCCATTATACCAGGCATAATAATCGCGGCGATCCCTTTGGAACTGTAAACCATGCCAACATTAGCGCCGGTGAAATGAAGAGTATTAATCATGTAAGAGCCGAGGGTAACCAGCCAGCTCCCCCAGATAAAATATTG  >  NC_000913/2523650‑2523946

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