Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 262,114 C→A 46.4% D51Y (GAT→TAT)  yafY ← lipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273262,1140CA46.4% 192.4 / inf 297D51Y (GAT→TAT) yafYlipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (69/90);  new base A (66/72);  total (135/162)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

ATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTT  >  CP009273/261965‑262260
                                                                                                                                                     |                                                                                                                                                  
   catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatcata                                                                                                                                                 >  1:258735/1‑149 (MQ=255)
    atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATcc                                                                                                                                                >  1:1485200/1‑150 (MQ=255)
     tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCt                                                                                                                                               >  2:1812552/1‑150 (MQ=255)
     tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCt                                                                                                                                               >  2:1813311/1‑150 (MQ=255)
      caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                              <  1:5104301/150‑1 (MQ=255)
      caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                              >  1:1850087/1‑150 (MQ=255)
      caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtg                                                                                                                                              <  2:265894/150‑1 (MQ=255)
         tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggtg                                                                                                                                           <  2:230370/150‑1 (MQ=255)
          ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggtgg                                                                                                                                          >  1:4622155/1‑150 (MQ=255)
           gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGc                                                                                                                                         <  1:2684601/150‑1 (MQ=255)
             tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTg                                                                                                                                       <  1:4478848/150‑1 (MQ=255)
              ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGa                                                                                                                                      >  2:1847426/1‑150 (MQ=255)
              ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGa                                                                                                                                      >  2:130574/1‑150 (MQ=255)
               gACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGAc                                                                                                                                     >  2:71867/1‑150 (MQ=255)
                aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGAcc                                                                                                                                    <  1:929874/150‑1 (MQ=255)
                 cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACca                                                                                                                                   <  1:4700369/150‑1 (MQ=255)
                   gTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACcaca                                                                                                                                 <  1:3452667/150‑1 (MQ=255)
                      ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCa                                                                                                                              >  2:2759475/1‑150 (MQ=255)
                         tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTg                                                                                                                           >  1:1284594/1‑150 (MQ=255)
                          ccTGCATATCCAAGAGGTAGTAGTAACTGCCATCCTCGCAGCCATAGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGAACATTCTTCGCTGTCATCATCCTGGTGGCTGACCACATCACTGa                                                                                                                          <  1:1765712/150‑1 (MQ=255)
                             gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAc                                                                                                                       <  1:2099601/150‑1 (MQ=255)
                                tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGAGGCTGACCACATCACTGAAGcgcg                                                                                                                    >  1:4803763/1‑150 (MQ=255)
                                   ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTg                                                                                                                 <  1:838619/150‑1 (MQ=255)
                                   ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTg                                                                                                                 <  2:1214428/150‑1 (MQ=255)
                                   ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTg                                                                                                                 >  1:4861539/1‑150 (MQ=255)
                                     aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                                                                               >  1:2721006/1‑150 (MQ=255)
                                     aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                                                                               <  2:4835546/150‑1 (MQ=255)
                                       gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                                                                             >  1:3288996/1‑150 (MQ=255)
                                       gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                                                                             <  1:1161135/150‑1 (MQ=255)
                                           tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAgg                                                                                                         >  1:3350268/1‑150 (MQ=255)
                                           tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAgg                                                                                                         <  1:265434/150‑1 (MQ=255)
                                              tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc                                                                                                      <  1:2830433/150‑1 (MQ=255)
                                               agtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCt                                                                                                     >  1:3985942/1‑150 (MQ=255)
                                                 aaaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTa                                                                                                   >  2:3317826/2‑150 (MQ=255)
                                                  aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                                                                  >  1:4721339/1‑150 (MQ=255)
                                                   aCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaaa                                                                                                 <  2:5085419/150‑1 (MQ=255)
                                                        cATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGACCCGCGGTGGGGGGAGGTCCTTAAATTTAc                                                                                            >  2:2099038/1‑150 (MQ=255)
                                                          tCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTg                                                                                          >  1:981644/1‑150 (MQ=255)
                                                           ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                         <  1:4450591/150‑1 (MQ=255)
                                                                 aGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                                                                   <  1:142089/150‑1 (MQ=255)
                                                                    cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                                                                >  1:999469/1‑150 (MQ=255)
                                                                    cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAACTTTACTGATGGCCTTCa                                                                                >  2:3730739/1‑150 (MQ=255)
                                                                    cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                                                                <  2:5056204/150‑1 (MQ=255)
                                                                      tCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAgg                                                                              <  2:1741848/150‑1 (MQ=255)
                                                                         ggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCTCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                                                                           <  1:2352021/150‑1 (MQ=255)
                                                                         ggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                                                                           >  2:3848222/1‑150 (MQ=255)
                                                                          gCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                                                                          <  2:2541161/150‑1 (MQ=255)
                                                                           cGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcata                                                                         <  1:3002945/150‑2 (MQ=255)
                                                                            gACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatca                                                                        <  2:1348584/150‑1 (MQ=255)
                                                                              cTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt                                                                      >  1:2388484/1‑150 (MQ=255)
                                                                              cTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt                                                                      >  1:1526002/1‑150 (MQ=255)
                                                                                 gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                                                                   <  2:3191917/150‑1 (MQ=255)
                                                                                   gggTACGTAAAACCTGCAAGGTGGGACCGGTTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                                 <  1:2911707/150‑1 (MQ=255)
                                                                                   gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                                 <  2:3717592/150‑1 (MQ=255)
                                                                                   gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                                 <  1:3917296/150‑1 (MQ=255)
                                                                                   gggTAAGTATTACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                                                                 <  1:2166811/150‑1 (MQ=255)
                                                                                     gTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                                                               >  1:4131759/1‑150 (MQ=255)
                                                                                     gTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                                                               >  2:1950354/1‑150 (MQ=255)
                                                                                      tACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAt                                                              >  1:768866/1‑150 (MQ=255)
                                                                                      tACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAt                                                              >  2:614398/1‑150 (MQ=255)
                                                                                      tACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCGGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAt                                                              <  2:2241452/150‑1 (MQ=255)
                                                                                       aCGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc                                                             <  2:1554266/150‑1 (MQ=255)
                                                                                        cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                                                            >  2:187578/1‑150 (MQ=255)
                                                                                        cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                                                            <  1:5137798/150‑1 (MQ=255)
                                                                                         gTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                                                           >  2:779672/1‑150 (MQ=255)
                                                                                           aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCGTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGctt                                                         >  2:1147852/1‑148 (MQ=255)
                                                                                           aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt                                                         <  1:553256/150‑1 (MQ=255)
                                                                                             aaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGc                                                       <  1:3779119/150‑1 (MQ=255)
                                                                                                cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                    <  2:3895382/150‑1 (MQ=255)
                                                                                                cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                    <  2:2598190/150‑1 (MQ=255)
                                                                                                cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                    <  1:2289332/150‑1 (MQ=255)
                                                                                                cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                    >  2:680841/1‑150 (MQ=255)
                                                                                                      ggTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaaa                                              >  2:2990672/1‑150 (MQ=255)
                                                                                                        tGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACa                                            >  2:4459440/1‑150 (MQ=255)
                                                                                                        tGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACa                                            <  1:2425755/150‑1 (MQ=255)
                                                                                                         gggACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                           >  1:436049/1‑150 (MQ=255)
                                                                                                              cGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtgg                                      <  2:2163566/150‑1 (MQ=255)
                                                                                                              cGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtgg                                      >  2:2062296/1‑150 (MQ=255)
                                                                                                              cGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtgg                                      <  1:4977906/150‑1 (MQ=255)
                                                                                                                gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                    <  2:4924562/150‑1 (MQ=255)
                                                                                                                 aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGtggtgg                                   >  2:4714242/1‑150 (MQ=255)
                                                                                                                      tagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAcc                              >  1:4266266/1‑150 (MQ=255)
                                                                                                                        gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAg                            <  1:1922580/150‑1 (MQ=255)
                                                                                                                        gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAg                            >  2:3100341/1‑150 (MQ=255)
                                                                                                                         tCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAgc                           <  1:4989163/150‑1 (MQ=255)
                                                                                                                            aCCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGcc                        <  2:4364546/150‑1 (MQ=255)
                                                                                                                             ccTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCca                       <  1:2116130/150‑1 (MQ=255)
                                                                                                                              cTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcag                      <  2:4301775/150‑1 (MQ=255)
                                                                                                                                     cATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGc               >  2:110092/1‑150 (MQ=255)
                                                                                                                                       ttcttcGCTGACATCCTCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAAAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg             <  2:2124357/150‑1 (MQ=255)
                                                                                                                                       ttcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg             <  1:3702567/150‑1 (MQ=255)
                                                                                                                                       ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg             <  1:1017747/150‑1 (MQ=255)
                                                                                                                                          tncGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCCTCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtgt          >  1:1538257/1‑150 (MQ=255)
                                                                                                                                           tcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGtt         <  2:4355278/150‑1 (MQ=255)
                                                                                                                                           tcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGtt         <  2:4039916/150‑1 (MQ=255)
                                                                                                                                               tGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc     >  2:947149/1‑150 (MQ=255)
                                                                                                                                               tGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc     <  2:350400/150‑1 (MQ=255)
                                                                                                                                               tGACATCATCCTGGTGGCGGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc     <  1:3214561/150‑1 (MQ=255)
                                                                                                                                                gACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACg    <  2:2446630/150‑1 (MQ=255)
                                                                                                                                                 aCATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGt   <  2:517154/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTT  >  CP009273/261965‑262260

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: