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,3720CA96.1% 702.2 / 15.8 231pseudogene (565/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base A (106/116);  total (112/119)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 232 total aligned reads displayed.

GATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646225‑1646517
                                                                                                                                                   |                                                                                                                                                 
   aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGcc                                                                                                                                                >  2:1720154/1‑148 (MQ=255)
   aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGcc                                                                                                                                                <  1:3420865/148‑1 (MQ=255)
    aTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGccc                                                                                                                                               >  2:2671541/1‑148 (MQ=255)
    aTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGccc                                                                                                                                               >  1:1659836/1‑148 (MQ=255)
     tCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCt                                                                                                                                              >  2:218900/1‑148 (MQ=255)
     tCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCt                                                                                                                                              >  2:2285207/1‑148 (MQ=255)
       cGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTgg                                                                                                                                            <  1:806900/148‑1 (MQ=255)
       cGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTgg                                                                                                                                            <  2:1967933/148‑1 (MQ=255)
       cGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTgg                                                                                                                                            >  1:1719937/1‑148 (MQ=255)
           tAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAg                                                                                                                                        >  2:3012330/1‑148 (MQ=255)
            aaGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGc                                                                                                                                       <  1:2914196/148‑1 (MQ=255)
             aGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCt                                                                                                                                      <  1:2492261/148‑1 (MQ=255)
             aGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCt                                                                                                                                      >  2:1551126/1‑148 (MQ=255)
              ggggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtg                                                                                                                                     >  1:165157/1‑148 (MQ=255)
               gggATCTCGCCATTAAAGGAATCAAGAACCATTACGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAAATTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgt                                                                                                                                    >  2:3035844/1‑148 (MQ=255)
                ggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtg                                                                                                                                   <  2:2356485/148‑1 (MQ=255)
                 gATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgt                                                                                                                                  >  2:417494/1‑148 (MQ=255)
                  aTCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 <  1:2171268/148‑1 (MQ=255)
                   tctcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGc                                                                                                                                <  1:2045142/148‑1 (MQ=255)
                    ctcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCa                                                                                                                               <  2:3067188/148‑1 (MQ=255)
                     tcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCAc                                                                                                                              <  2:3491921/148‑1 (MQ=255)
                      cGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACt                                                                                                                             >  2:3216581/1‑148 (MQ=255)
                      cGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACt                                                                                                                             <  1:2383151/148‑1 (MQ=255)
                          aTTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGa                                                                                                                         >  2:2112414/1‑148 (MQ=255)
                            tAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCTCTTGGATg                                                                                                                       <  2:873439/148‑1 (MQ=255)
                             aaaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATgg                                                                                                                      <  1:1788230/148‑1 (MQ=255)
                                ggAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggag                                                                                                                   <  1:1669626/148‑1 (MQ=255)
                                ggAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggag                                                                                                                   <  1:1244875/148‑1 (MQ=255)
                                   aTCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGAc                                                                                                                <  1:2212335/148‑1 (MQ=255)
                                        gAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGa                                                                                                           >  1:142193/1‑148 (MQ=255)
                                          aCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATc                                                                                                         <  2:3329406/148‑1 (MQ=255)
                                             aTTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGa                                                                                                      <  1:1109997/148‑1 (MQ=255)
                                              ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                     >  1:2097078/1‑148 (MQ=255)
                                              ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                     >  1:1900057/1‑148 (MQ=255)
                                              ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                     <  1:2292762/148‑1 (MQ=255)
                                                  ggATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAAt                                                                                                 <  2:3477949/148‑1 (MQ=255)
                                                    aTCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATgg                                                                                               >  1:935446/1‑148 (MQ=255)
                                                     tCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGa                                                                                              <  2:2915252/148‑1 (MQ=255)
                                                        tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           >  1:622535/1‑148 (MQ=255)
                                                        tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           >  2:556200/1‑148 (MQ=255)
                                                        tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           <  1:376929/148‑1 (MQ=255)
                                                         ctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACaa                                                                                          <  2:920549/148‑1 (MQ=255)
                                                          tttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAat                                                                                         >  2:424033/1‑148 (MQ=255)
                                                           ttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatg                                                                                        <  2:2330361/148‑1 (MQ=255)
                                                            tCGCTGGAAGAATTTATCCGTCTTATAGATGCGTGCAGAAATCAGAAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGAAATAGATCTGAAAAATGGAACAatga                                                                                       <  1:479039/148‑1 (MQ=255)
                                                            tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  2:859900/148‑1 (MQ=255)
                                                             cgcgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatgat                                                                                      >  2:451512/1‑148 (MQ=255)
                                                              gcgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatgatg                                                                                     <  2:1041342/148‑1 (MQ=255)
                                                               cgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAAtgatga                                                                                    <  2:680956/148‑1 (MQ=255)
                                                                gggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATgatgat                                                                                   <  1:87048/148‑1 (MQ=255)
                                                                gggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATgatgat                                                                                   <  1:2804226/148‑1 (MQ=255)
                                                                    agaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAgg                                                                               >  2:2047330/1‑148 (MQ=255)
                                                                     gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                              >  2:3156758/1‑148 (MQ=255)
                                                                     gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                              >  1:2143883/1‑148 (MQ=255)
                                                                     gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                              >  2:163283/1‑148 (MQ=255)
                                                                       aTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGgaga                                                                            >  1:2987122/1‑148 (MQ=255)
                                                                        tttATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaa                                                                           <  1:1358168/148‑1 (MQ=255)
                                                                        tttATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaa                                                                           <  2:370763/148‑1 (MQ=255)
                                                                                  ttATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaaa                                                                 <  2:723269/148‑1 (MQ=255)
                                                                                    aTCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAg                                                               <  2:929097/148‑1 (MQ=255)
                                                                                    aTCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAg                                                               <  2:2143883/148‑1 (MQ=255)
                                                                                      cGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGAc                                                             <  2:2518409/148‑1 (MQ=255)
                                                                                                aGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACgg                                                   <  1:2017062/148‑1 (MQ=255)
                                                                                                aGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACgg                                                   >  2:2419183/1‑148 (MQ=255)
                                                                                                 gAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGt                                                  >  1:1014471/1‑148 (MQ=255)
                                                                                                  aaaTCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTa                                                 >  1:2729968/1‑148 (MQ=255)
                                                                                                     tCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCa                                              <  1:1327319/148‑1 (MQ=255)
                                                                                                      cAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGACATTTAGCAAAAGACCGTTTCACGGTACCaa                                             >  2:2097193/1‑148 (MQ=255)
                                                                                                      cAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaa                                             >  2:3451799/1‑148 (MQ=255)
                                                                                                        gcaagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaa                                           <  2:1770223/148‑1 (MQ=255)
                                                                                                        gcaagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaa                                           >  1:2521035/1‑148 (MQ=255)
                                                                                                         caagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaaa                                          >  1:2868688/1‑148 (MQ=255)
                                                                                                            gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       <  2:1092463/148‑1 (MQ=255)
                                                                                                            gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       <  1:3351276/148‑1 (MQ=255)
                                                                                                            gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       <  1:1750634/148‑1 (MQ=255)
                                                                                                            gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       >  2:2233621/1‑148 (MQ=255)
                                                                                                               aaaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAgg                                    >  2:1701271/1‑148 (MQ=255)
                                                                                                                  aaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCgg                                 >  2:2623851/1‑148 (MQ=255)
                                                                                                                     ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                              >  2:966336/1‑148 (MQ=255)
                                                                                                                     ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                              >  2:2246895/1‑148 (MQ=255)
                                                                                                                       aTGGTGTTTTTCTGTTTATACTGGAGTTAGCCCTTGTGAGCTGTGTGCACTTTGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTGACGGTACCAAAAACACAGGCGGGAAcc                            <  2:3137161/148‑1 (MQ=255)
                                                                                                                         ggTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCaa                          >  2:1300278/1‑148 (MQ=255)
                                                                                                                            gtgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCg                       >  1:1514748/1‑148 (MQ=255)
                                                                                                                              gtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCggg                     >  1:655442/1‑148 (MQ=255)
                                                                                                                               tttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGt                    >  1:2085784/1‑148 (MQ=255)
                                                                                                                               tttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGt                    <  2:1504123/148‑1 (MQ=255)
                                                                                                                               tttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGt                    >  1:2614056/1‑148 (MQ=255)
                                                                                                                                ttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtc                   >  2:3449263/1‑148 (MQ=255)
                                                                                                                                   tGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcatt                <  1:3300667/148‑1 (MQ=255)
                                                                                                                                     tttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattca              >  2:1603205/1‑148 (MQ=255)
                                                                                                                                     tttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattca              >  2:1080530/1‑148 (MQ=255)
                                                                                                                                      ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat             <  1:3483658/148‑1 (MQ=255)
                                                                                                                                       tataCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc            <  2:3378723/148‑1 (MQ=255)
                                                                                                                                       tataCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc            <  2:2017427/148‑1 (MQ=255)
                                                                                                                                         taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt          <  1:1156070/148‑1 (MQ=255)
                                                                                                                                          aCTGGAGTTCGCCCTGGTGAGCTGGGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         >  2:754638/1‑148 (MQ=255)
                                                                                                                                          aCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         >  2:3526389/1‑148 (MQ=255)
                                                                                                                                          aCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta         <  2:3126109/148‑1 (MQ=255)
                                                                                                                                            tGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCttatt       >  2:1672414/1‑148 (MQ=255)
                                                                                                                                                 ttCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGcc  <  1:3216552/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646225‑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: