Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 579,853 C→A 35.9% G239G (GGC→GGA)  appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273579,8530CA35.9% 398.7 / 248.8 253G239G (GGC→GGA) appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (78/84);  new base A (40/51);  total (118/135)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 257 total aligned reads displayed.

AGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACA  >  CP009273/579704‑580001
                                                                                                                                                     |                                                                                                                                                    
aGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGc                                                                                                                                                      <  2:627770/150‑1 (MQ=255)
aGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGc                                                                                                                                                      >  2:1866511/1‑150 (MQ=255)
   tATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTc                                                                                                                                                   >  2:1845420/1‑150 (MQ=255)
      gCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTcaca                                                                                                                                                <  1:426242/150‑1 (MQ=255)
       cAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAg                                                                                                                                               >  1:2496357/1‑150 (MQ=255)
        aaaaaaaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGa                                                                                                                                              >  1:4042183/1‑150 (MQ=255)
        aaaaaaaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGa                                                                                                                                              >  1:1871831/1‑150 (MQ=255)
        aaaaaaaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGa                                                                                                                                              <  2:361299/150‑1 (MQ=255)
         aaaaaaCTCATAACTTCAAACTCTTATTCTNTCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGAt                                                                                                                                             >  1:4283872/1‑150 (MQ=255)
          aaaaaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATg                                                                                                                                            >  2:1142894/1‑150 (MQ=255)
            aaaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGa                                                                                                                                          >  2:4077448/1‑150 (MQ=255)
             aaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGaa                                                                                                                                         <  1:413761/150‑1 (MQ=255)
             aaCTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGaa                                                                                                                                         <  1:4641014/150‑1 (MQ=255)
             aaCTCATAACTTCAAACTCTTATACACTCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTACATTTTTCATATGGGCATTTAAAGATTATTATGGTTTCACGCCATCTCATTAATTGGCGAAAATAATCGTCGTCACAGATGGaa                                                                                                                                         >  2:2708752/1‑150 (MQ=255)
               cTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATa                                                                                                                                       <  2:888791/150‑1 (MQ=255)
                  aTAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATaaac                                                                                                                                    <  1:2066164/150‑1 (MQ=255)
                   tAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATaaaca                                                                                                                                   >  2:4458020/1‑150 (MQ=255)
                    aaCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATaaacaa                                                                                                                                  <  2:3202098/150‑1 (MQ=255)
                    aaCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATaaacaa                                                                                                                                  >  1:526944/1‑150 (MQ=255)
                     aCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATaaacaaa                                                                                                                                 >  1:1313835/1‑150 (MQ=255)
                     aCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATaaacaaa                                                                                                                                 >  1:1313824/1‑150 (MQ=255)
                     aCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATaaacaaa                                                                                                                                 >  2:401790/1‑150 (MQ=255)
                     aCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATaaacaaa                                                                                                                                 <  1:3044667/150‑1 (MQ=255)
                       ttCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATaaacaaaac                                                                                                                               >  1:3038392/1‑150 (MQ=255)
                          aaaCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAAt                                                                                                                            <  1:384764/150‑1 (MQ=255)
                           aaCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAAtt                                                                                                                           <  1:1096922/150‑1 (MQ=255)
                           aaCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAAtt                                                                                                                           <  1:3523450/150‑1 (MQ=255)
                            aCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTg                                                                                                                          <  2:3579247/150‑1 (MQ=255)
                            aCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTg                                                                                                                          <  1:1065663/150‑1 (MQ=255)
                             ctctTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGa                                                                                                                         <  2:959124/150‑1 (MQ=255)
                             ctctTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGa                                                                                                                         >  1:2320250/1‑150 (MQ=255)
                              tctTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGAc                                                                                                                        <  1:3910404/150‑1 (MQ=255)
                              tctTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGAc                                                                                                                        <  2:3179116/150‑1 (MQ=255)
                               ctTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACt                                                                                                                       <  1:4512003/150‑1 (MQ=255)
                                ttATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTg                                                                                                                      <  1:995489/150‑1 (MQ=255)
                                ttATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTg                                                                                                                      >  1:2216693/1‑150 (MQ=255)
                                  aTTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGAt                                                                                                                    >  1:3263395/1‑150 (MQ=255)
                                  aTTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGAt                                                                                                                    >  1:2370367/1‑150 (MQ=255)
                                  aTTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGAt                                                                                                                    <  1:2344274/150‑1 (MQ=255)
                                     cTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAAt                                                                                                                 <  1:1919076/150‑1 (MQ=255)
                                     cTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAAt                                                                                                                 <  1:3188915/150‑1 (MQ=255)
                                       aTCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGt                                                                                                               >  1:4050700/1‑150 (MQ=255)
                                         cAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGttt                                                                                                             >  1:2325385/1‑150 (MQ=255)
                                             gTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTttatta                                                                                                         >  2:3780771/1‑150 (MQ=255)
                                             gTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTttatta                                                                                                         >  1:3235098/1‑150 (MQ=255)
                                             gTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGAGTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTttatta                                                                                                         >  2:228532/1‑150 (MQ=255)
                                              tCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTAc                                                                                                        >  2:3041305/1‑150 (MQ=255)
                                              tCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTAc                                                                                                        <  1:2650499/150‑1 (MQ=255)
                                                gTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACaa                                                                                                      <  1:2609481/150‑1 (MQ=255)
                                                    cccAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTAAAAGTTg                                                                                                  >  2:4352606/1‑150 (MQ=255)
                                                     ccAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGt                                                                                                 <  1:212033/150‑1 (MQ=255)
                                                       aGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACTCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCt                                                                                               <  1:3758861/150‑1 (MQ=255)
                                                       aGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCt                                                                                               <  1:1320419/150‑1 (MQ=255)
                                                           aTGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACAt                                                                                           >  1:2900442/1‑150 (MQ=255)
                                                            tgtgGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAAGAATCGGCGTCACAGATGGACTAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATg                                                                                          >  2:842077/1‑150 (MQ=255)
                                                                gCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTaa                                                                                      <  2:2393422/150‑1 (MQ=255)
                                                                  tataACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAAtt                                                                                    <  1:202389/150‑1 (MQ=255)
                                                                    taACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATtat                                                                                  <  2:1181325/150‑1 (MQ=255)
                                                                    taACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATtat                                                                                  <  2:2866245/150‑1 (MQ=255)
                                                                      aCAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTataa                                                                                >  2:2493169/1‑150 (MQ=255)
                                                                      aCAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTataa                                                                                >  2:661351/1‑150 (MQ=255)
                                                                       cAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTataat                                                                               >  2:1813879/1‑150 (MQ=255)
                                                                        aGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTataata                                                                              >  1:1423025/1‑150 (MQ=255)
                                                                         gTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAatat                                                                             >  2:1545506/1‑150 (MQ=255)
                                                                          tACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAAtatt                                                                            <  1:3698737/150‑1 (MQ=255)
                                                                           aCTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAAtatta                                                                           >  2:150950/1‑150 (MQ=255)
                                                                            cTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAAtattat                                                                          >  1:2973017/1‑150 (MQ=255)
                                                                              tcatatttcatatGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATAc                                                                        <  1:352713/150‑1 (MQ=255)
                                                                                   tttcatatGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGt                                                                   <  2:1647775/150‑1 (MQ=255)
                                                                                    ttcatatGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGtt                                                                  <  2:3423302/150‑1 (MQ=255)
                                                                                        tatGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGtttttt                                                              >  2:1244319/1‑150 (MQ=255)
                                                                                          tgtgCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGtttttttt                                                            >  2:2174423/1‑150 (MQ=255)
                                                                                           gtgCATTTAAAGATTATTCTGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTg                                                           <  2:399490/150‑1 (MQ=255)
                                                                                             gCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGAt                                                         <  2:2379893/150‑1 (MQ=255)
                                                                                               aTTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGAtgt                                                       >  2:3464788/1‑150 (MQ=255)
                                                                                                  tAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGat                                                    <  1:329011/150‑1 (MQ=255)
                                                                                                  tAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGat                                                    >  2:3337225/1‑150 (MQ=255)
                                                                                                    aaGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGatat                                                  >  2:3303873/1‑150 (MQ=255)
                                                                                                     cgATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATAtt                                                 <  1:2541086/149‑1 (MQ=255)
                                                                                                     aGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATAtt                                                 >  2:2003527/1‑150 (MQ=255)
                                                                                                     aGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATAtt                                                 <  1:2522871/150‑1 (MQ=255)
                                                                                                       attattATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCt                                               <  2:1822066/150‑1 (MQ=255)
                                                                                                       attattATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCt                                               <  1:1128284/150‑1 (MQ=255)
                                                                                                               ggTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCAt                                       >  1:833526/1‑150 (MQ=255)
                                                                                                                     aCGCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAAttt                                 >  1:1545731/1‑150 (MQ=255)
                                                                                                                       gCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGt                               >  1:4049116/1‑150 (MQ=255)
                                                                                                                       gCCATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGt                               >  1:2831176/1‑150 (MQ=255)
                                                                                                                        ccATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTa                              >  2:774596/1‑150 (MQ=255)
                                                                                                                         cATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTaa                             <  1:3265650/150‑1 (MQ=255)
                                                                                                                         cATCTCATTATTTTGAGAAAATAATCGGAGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTaa                             <  2:526944/150‑1 (MQ=255)
                                                                                                                                ttattTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGtt                      <  1:1865176/150‑1 (MQ=255)
                                                                                                                                ttattTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGtt                      <  2:693850/150‑1 (MQ=255)
                                                                                                                                  attTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAATACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAAGGATTAATTTGTAATTGGGTTGc                    >  2:2221911/1‑150 (MQ=255)
                                                                                                                                    tttGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTAACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTg                  <  2:2648836/150‑1 (MQ=255)
                                                                                                                                     ttGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGt                 <  2:4274017/150‑1 (MQ=255)
                                                                                                                                           aaaTAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTa           <  2:2856609/150‑1 (MQ=255)
                                                                                                                                           aaaTAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTa           >  1:3936906/1‑150 (MQ=255)
                                                                                                                                              taatCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTAttt        <  1:1959556/150‑1 (MQ=255)
                                                                                                                                               aatCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTAtttt       <  1:2169909/150‑1 (MQ=255)
                                                                                                                                               aatCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTAtttt       <  1:1195207/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGAAAATAATCGGCGTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACA  >  CP009273/579704‑580001

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: ‑