Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,646,372 C→A 100% pseudogene (565/1158 nt) intQ → pseudogene, Qin prophage; phage integrase family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,646,3720CA95.6% 682.4 / 18.7 226pseudogene (565/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base A (119/97);  total (125/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 226 total aligned reads displayed.

ATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646226‑1646517
                                                                                                                                                  |                                                                                                                                                 
   aTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCt                                                                                                                                              >  2:1974439/1‑149 (MQ=255)
   aTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCt                                                                                                                                              <  1:470544/149‑1 (MQ=255)
       gTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTg                                                                                                                                          <  2:143021/149‑1 (MQ=255)
        tttAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGa                                                                                                                                         >  2:2373252/1‑149 (MQ=255)
          tAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGc                                                                                                                                       >  2:1778740/1‑149 (MQ=255)
          tAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGc                                                                                                                                       <  1:2244649/149‑1 (MQ=255)
           aaGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCt                                                                                                                                      >  2:1865934/1‑149 (MQ=255)
            aGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtg                                                                                                                                     >  1:2420596/1‑149 (MQ=255)
             ggggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgt                                                                                                                                    <  1:3843118/149‑1 (MQ=255)
              gggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtg                                                                                                                                   <  2:2089266/149‑1 (MQ=255)
                gATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 >  2:1700257/1‑149 (MQ=255)
                gATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 <  2:265146/149‑1 (MQ=255)
                gATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 >  2:3306655/1‑149 (MQ=255)
                gATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 >  1:2811084/1‑149 (MQ=255)
                  tctcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCa                                                                                                                               >  2:1776411/1‑149 (MQ=255)
                      gCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTg                                                                                                                           <  2:3602564/149‑1 (MQ=255)
                       ccATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTgg                                                                                                                          >  1:3812475/1‑149 (MQ=255)
                        cATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGa                                                                                                                         <  1:2166164/149‑1 (MQ=255)
                        cATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGa                                                                                                                         >  2:1357356/1‑149 (MQ=255)
                         aTTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACt                                                                                                                             <  1:319112/144‑1 (MQ=255)
                         aTTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACt                                                                                                                             >  2:319112/1‑144 (MQ=255)
                          ttAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATg                                                                                                                       >  1:1388731/1‑149 (MQ=255)
                            aaaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATggg                                                                                                                     <  2:3315187/149‑1 (MQ=255)
                            aaaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATggg                                                                                                                     <  2:2503208/149‑1 (MQ=255)
                              aGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggag                                                                                                                   >  1:3387786/1‑149 (MQ=255)
                               ggAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggagg                                                                                                                  <  2:3354198/149‑1 (MQ=255)
                                gAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggagga                                                                                                                 <  1:1010332/149‑1 (MQ=255)
                                 aaTCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGAc                                                                                                                <  1:960028/149‑1 (MQ=255)
                                   tCAAGAACCATTCCTGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATAGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGAACTTGGATGGGAGGACAt                                                                                                              <  1:2893286/149‑1 (MQ=255)
                                   tCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACAt                                                                                                              <  2:3708258/149‑1 (MQ=255)
                                    cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATa                                                                                                             >  1:1343186/1‑149 (MQ=255)
                                    cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATa                                                                                                             <  1:3492371/149‑1 (MQ=255)
                                     aaGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAg                                                                                                            <  2:1537986/149‑1 (MQ=255)
                                      aGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGa                                                                                                           <  1:2218885/149‑1 (MQ=255)
                                      aGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTGTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGa                                                                                                           >  2:3797204/1‑149 (MQ=255)
                                        aaCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATc                                                                                                         >  2:317938/1‑149 (MQ=255)
                                         aCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCt                                                                                                        <  1:2725889/149‑1 (MQ=255)
                                         aCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCt                                                                                                        >  1:1253850/1‑149 (MQ=255)
                                         aCCAGTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCt                                                                                                        <  2:2606894/149‑1 (MQ=255)
                                                 ggATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATg                                                                                                >  2:595581/1‑149 (MQ=255)
                                                 ggATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATg                                                                                                >  2:916738/1‑149 (MQ=255)
                                                 ggATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATGTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATg                                                                                                >  2:3615995/1‑149 (MQ=255)
                                                  gATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATgg                                                                                               >  1:2190238/1‑149 (MQ=255)
                                                   aTCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGa                                                                                              <  1:2332143/149‑1 (MQ=255)
                                                   aTCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGa                                                                                              >  2:2032586/1‑149 (MQ=255)
                                                   aTCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGa                                                                                              >  2:4068150/1‑149 (MQ=255)
                                                    tCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGaa                                                                                             >  2:3928097/1‑149 (MQ=255)
                                                     ccTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAAc                                                                                            >  2:3719439/1‑149 (MQ=255)
                                                     ccTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAAc                                                                                            >  2:1156175/1‑149 (MQ=255)
                                                      ctctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           >  1:3197830/1‑149 (MQ=255)
                                                       tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACaa                                                                                          <  1:3725927/149‑1 (MQ=255)
                                                         tttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatg                                                                                        >  1:1189253/1‑149 (MQ=255)
                                                          ttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCCATCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTTTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  1:1823452/149‑1 (MQ=255)
                                                          ttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  2:3789669/149‑1 (MQ=255)
                                                           tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatgat                                                                                      <  1:1865547/149‑1 (MQ=255)
                                                                 gaagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAg                                                                                <  2:688791/149‑1 (MQ=255)
                                                                  aagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAgg                                                                               >  1:3426486/1‑149 (MQ=255)
                                                                   agaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                              >  1:3055825/1‑149 (MQ=255)
                                                                    gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGgag                                                                             >  2:2752626/1‑149 (MQ=255)
                                                                      aTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaa                                                                           >  2:3751294/1‑149 (MQ=255)
                                                                       tttATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaaa                                                                          >  1:421844/1‑149 (MQ=255)
                                                                         tATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAtt                                                                        <  2:2615843/149‑1 (MQ=255)
                                                                                   aTCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGa                                                              >  2:1279786/1‑149 (MQ=255)
                                                                                   aTCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGa                                                              <  1:1896125/149‑1 (MQ=255)
                                                                                    tCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGAc                                                             >  2:3221835/1‑149 (MQ=255)
                                                                                       aTGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGt                                                          <  1:661433/149‑1 (MQ=255)
                                                                                           gTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCa                                                      >  1:538218/1‑149 (MQ=255)
                                                                                            tGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCAc                                                     <  1:36321/149‑1 (MQ=255)
                                                                                            tGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCAc                                                     >  1:3698689/1‑149 (MQ=255)
                                                                                                gAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTa                                                 <  1:3472219/149‑1 (MQ=255)
                                                                                                gAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTa                                                 >  1:4025342/1‑149 (MQ=255)
                                                                                                gAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGATGTGTGCACTTTGATGTGATGAAATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTa                                                 <  1:1564030/149‑1 (MQ=255)
                                                                                                  aaTCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTaac                                               >  2:4085902/1‑147 (MQ=255)
                                                                                                      agcaagcaAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaa                                           >  1:3271676/1‑149 (MQ=255)
                                                                                                       gcaagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaaa                                          >  1:3776018/1‑149 (MQ=255)
                                                                                                         aagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCTAAAGACCGTTTCACGGTACCAAAAaca                                        >  1:1416161/1‑149 (MQ=255)
                                                                                                          agcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       >  1:1464496/1‑149 (MQ=255)
                                                                                                            caaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAg                                     <  2:2811084/149‑1 (MQ=255)
                                                                                                             aaaaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAgg                                    >  1:704097/1‑149 (MQ=255)
                                                                                                              aaaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGc                                   >  1:486954/1‑149 (MQ=255)
                                                                                                                aaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCgg                                 >  2:3072963/1‑149 (MQ=255)
                                                                                                                 aaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCggg                                >  2:3009459/1‑149 (MQ=255)
                                                                                                                   tttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                              >  1:820327/1‑149 (MQ=255)
                                                                                                                     tATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAcc                            >  1:1722748/1‑149 (MQ=255)
                                                                                                                      aTGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCa                           >  1:3879016/1‑149 (MQ=255)
                                                                                                                        ggTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAt                         <  1:2548605/149‑1 (MQ=255)
                                                                                                                           gtgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCgg                      <  1:2373252/149‑1 (MQ=255)
                                                                                                                              tttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtc                   >  1:912453/1‑149 (MQ=255)
                                                                                                                              tttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtc                   <  1:1610268/149‑1 (MQ=255)
                                                                                                                               ttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtca                  <  2:89568/149‑1 (MQ=255)
                                                                                                                                  tGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattc               >  1:2147061/1‑149 (MQ=255)
                                                                                                                                  tGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattc               >  1:642203/1‑149 (MQ=255)
                                                                                                                                    tttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat             >  2:2093295/1‑149 (MQ=255)
                                                                                                                                     ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc            <  2:3812475/149‑1 (MQ=255)
                                                                                                                                     ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc            <  1:61004/149‑1 (MQ=255)
                                                                                                                                     ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc            >  1:1057751/1‑149 (MQ=255)
                                                                                                                                        taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         >  2:3296486/1‑149 (MQ=255)
                                                                                                                                        taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         >  1:552885/1‑149 (MQ=255)
                                                                                                                                        taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         <  1:916738/149‑1 (MQ=255)
                                                                                                                                               gTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTATCACTGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAgaa  >  1:2471381/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646226‑1646517

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