Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,577,827 A→G 28.0% intergenic (+158/‑292) yrhA → / → yrhB pseudogene, interrupted by IS1E/stable heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,577,8270AG28.0% 641.6 / inf 415intergenic (+158/‑292)yrhA/yrhBpseudogene, interrupted by IS1E/stable heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (170/128);  new base G (53/63);  total (223/192)
Fisher's exact test for biased strand distribution p-value = 4.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 419 total aligned reads displayed.

AGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACT  >  CP009273/3577678‑3577976
                                                                                                                                                     |                                                                                                                                                     
     aCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTa                                                                                                                                                  <  2:3710693/150‑1 (MQ=255)
        aaaaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGcc                                                                                                                                               <  1:2447132/150‑1 (MQ=255)
           aTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGaa                                                                                                                                            >  2:146486/1‑150 (MQ=255)
             cATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAAt                                                                                                                                          <  1:66319/150‑1 (MQ=255)
             aaTTGTAAACCATGCCATCTTTTTAGAGATTACCAGCATAATATTCCTTGTTCTAAATTTAGCTGTGCATGCTATGTAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAAt                                                                                                                                          <  1:3140599/149‑1 (MQ=255)
               ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAATCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAg                                                                                                                                        >  2:5181140/1‑150 (MQ=255)
               ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAg                                                                                                                                        >  2:2091697/1‑150 (MQ=255)
               ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAg                                                                                                                                        >  1:2577838/1‑150 (MQ=255)
               ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAg                                                                                                                                        <  2:4652038/150‑1 (MQ=255)
                tGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGc                                                                                                                                       >  1:1499780/1‑150 (MQ=255)
                  tAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTa                                                                                                                                     >  1:1616959/1‑150 (MQ=255)
                  tAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTa                                                                                                                                     >  2:4972942/1‑150 (MQ=255)
                    aaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAAt                                                                                                                                   <  1:3581942/150‑1 (MQ=255)
                     aCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATg                                                                                                                                  <  2:1914170/150‑1 (MQ=255)
                     aCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATg                                                                                                                                  >  1:3083586/1‑150 (MQ=255)
                     aCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATg                                                                                                                                  <  1:2410243/150‑1 (MQ=255)
                        aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCt                                                                                                                               >  1:1902698/1‑150 (MQ=255)
                        aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCt                                                                                                                               >  2:323052/1‑150 (MQ=255)
                          gCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTat                                                                                                                             <  1:2241003/150‑1 (MQ=255)
                            cATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTatga                                                                                                                           >  1:1200377/1‑150 (MQ=255)
                            cATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatga                                                                                                                           <  1:1102369/150‑1 (MQ=255)
                            cATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatgat                                                                                                                          <  2:949811/150‑1 (MQ=255)
                             aTCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatgat                                                                                                                          >  2:3029642/1‑150 (MQ=255)
                              tCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCTGTGAAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTatgatg                                                                                                                         <  2:70897/150‑1 (MQ=255)
                               cTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGACCTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTAtgatga                                                                                                                        >  1:2624331/1‑150 (MQ=255)
                               cTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTAtgatga                                                                                                                        <  1:866739/150‑1 (MQ=255)
                                 ttttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGAtt                                                                                                                      >  2:3598617/1‑150 (MQ=255)
                                  tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAttt                                                                                                                     <  2:4825847/150‑1 (MQ=255)
                                   ttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAtttt                                                                                                                    <  2:4147420/150‑1 (MQ=255)
                                     agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                  >  1:3848147/1‑150 (MQ=255)
                                      gagaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAAt                                                                                                                 >  1:1746405/1‑150 (MQ=255)
                                      gagaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAAt                                                                                                                 <  2:5033236/150‑1 (MQ=255)
                                       agaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATg                                                                                                                <  1:4653016/150‑1 (MQ=255)
                                          ttACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACa                                                                                                             <  1:1364278/150‑1 (MQ=255)
                                           tACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACAt                                                                                                            >  2:2507325/1‑150 (MQ=255)
                                               agcagAAAATTCCTTGTACGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTg                                                                                                        <  2:1231375/150‑1 (MQ=255)
                                                gcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGa                                                                                                       <  2:265329/150‑1 (MQ=255)
                                                  agAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                     <  1:805828/150‑1 (MQ=255)
                                                     aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                  >  1:1675920/1‑150 (MQ=255)
                                                          ccTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGa                                                                                             >  1:5309010/1‑150 (MQ=255)
                                                           cTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGaa                                                                                            >  2:1877120/1‑150 (MQ=255)
                                                            ttGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAt                                                                                           <  1:4290684/150‑1 (MQ=255)
                                                            ttGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAt                                                                                           >  1:337324/1‑150 (MQ=255)
                                                             tGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtt                                                                                          <  2:1546070/150‑1 (MQ=255)
                                                              gTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAttt                                                                                         >  1:4331494/1‑150 (MQ=255)
                                                              gTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAttt                                                                                         <  1:695068/150‑1 (MQ=255)
                                                              gTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAttt                                                                                         >  1:5157349/1‑150 (MQ=255)
                                                               ttCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtttt                                                                                        >  1:401884/1‑150 (MQ=255)
                                                               ttCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtttt                                                                                        <  1:478620/150‑1 (MQ=255)
                                                                         aGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCa                                                                              >  1:4398452/1‑150 (MQ=255)
                                                                         aGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCa                                                                              >  1:3311630/1‑150 (MQ=255)
                                                                            aGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATaa                                                                           <  2:3067972/150‑1 (MQ=255)
                                                                                cAAGCTTTGGAAAAACTTACAACGGAACTACACTACTTAAGTGAGATGGACCGTAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATTGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGAtt                                                                       <  1:3695548/150‑1 (MQ=255)
                                                                                 aaGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTc                                                                      >  2:1098016/1‑150 (MQ=255)
                                                                                  aGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACAGGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCt                                                                     >  1:4839183/1‑150 (MQ=255)
                                                                                     tATTGAAAAACTTACTACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGat                                                                  <  2:1715893/150‑1 (MQ=255)
                                                                                      aTGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGata                                                                 <  2:130107/150‑1 (MQ=255)
                                                                                       tGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGatat                                                                <  2:3791497/150‑1 (MQ=255)
                                                                                         gAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCt                                                              >  2:4458496/1‑150 (MQ=255)
                                                                                          aaaaaCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCtt                                                             <  1:1111141/150‑1 (MQ=255)
                                                                                            aaaCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTaa                                                           >  1:4715181/1‑150 (MQ=255)
                                                                                              aCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAAt                                                         <  2:2806917/150‑1 (MQ=255)
                                                                                               cTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATa                                                        >  1:1810678/1‑150 (MQ=255)
                                                                                                ttACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATaa                                                       <  2:4584686/150‑1 (MQ=255)
                                                                                                 tACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAg                                                      <  1:4086971/150‑1 (MQ=255)
                                                                                                   cAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATTATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaa                                                    <  2:2709971/150‑1 (MQ=255)
                                                                                                     aCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGAGTCTGATATCTTAAATAAGaaaa                                                  >  2:2882556/1‑150 (MQ=255)
                                                                                                      cGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAAt                                                 >  2:4749260/1‑150 (MQ=255)
                                                                                                       ggAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                                >  2:1618917/1‑150 (MQ=255)
                                                                                                       ggAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                                <  1:5040242/150‑1 (MQ=255)
                                                                                                        gAACTACACTACTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATcc                                               <  1:5019758/150‑1 (MQ=255)
                                                                                                          aCTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCg                                             <  2:4038112/150‑1 (MQ=255)
                                                                                                            tACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGtt                                           >  1:4162257/1‑150 (MQ=255)
                                                                                                             acacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTc                                          >  2:4875876/1‑150 (MQ=255)
                                                                                                              cacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCg                                         <  2:1026359/150‑1 (MQ=255)
                                                                                                              cacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCg                                         <  2:341782/150‑1 (MQ=255)
                                                                                                               acTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGa                                        <  1:2923264/150‑1 (MQ=255)
                                                                                                                   cTTAAGTGAGAAGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACACGAATCATAAGATTCTGATATATTAAATAATAAAATCCCGTTCGACGtt                                    >  1:3779251/1‑150 (MQ=255)
                                                                                                                    ttAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTc                                   >  1:3932265/1‑150 (MQ=255)
                                                                                                                      aaGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTg                                 <  2:420053/150‑1 (MQ=255)
                                                                                                                         tGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaa                              >  1:3000074/1‑150 (MQ=255)
                                                                                                                         tGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaa                              <  1:2526677/150‑1 (MQ=255)
                                                                                                                         tGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAGTTTAATGACCTGTTGATTTATATAGAATTTTACCCGACTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaa                              >  1:196848/1‑150 (MQ=255)
                                                                                                                             aTGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGa                          >  1:3957702/1‑150 (MQ=255)
                                                                                                                               ggACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaaa                        >  1:2157772/1‑150 (MQ=255)
                                                                                                                                gACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAAt                       >  2:619771/1‑150 (MQ=255)
                                                                                                                                  ccGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATTGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCAGATCGACGTTCTGAAAAAGAAATCa                     >  1:2335744/1‑150 (MQ=255)
                                                                                                                                   cGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCAt                    <  2:4331494/150‑1 (MQ=255)
                                                                                                                                   cGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCAt                    <  1:2091697/150‑1 (MQ=255)
                                                                                                                                    gAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATa                   >  2:3525334/1‑150 (MQ=255)
                                                                                                                                      aGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaaa                 >  2:4014415/1‑150 (MQ=255)
                                                                                                                                       gACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaaaa                >  1:138241/1‑150 (MQ=255)
                                                                                                                                          aTGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAatat             <  1:3253327/150‑1 (MQ=255)
                                                                                                                                          aTGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAatat             <  2:855296/150‑1 (MQ=255)
                                                                                                                                            gTAGTCTCTTTTTTAGCCGAAATAGCTAATTTCTATGATTATTTTAATGACATTTTGATTTATTTAGAATATTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAAACATAAAATATTa           <  2:5261860/150‑1 (MQ=255)
                                                                                                                                             tAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAg          >  2:778264/1‑150 (MQ=255)
                                                                                                                                               gTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCa        >  1:1016673/1‑150 (MQ=255)
                                                                                                                                                 ctctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGa      >  2:5106993/1‑150 (MQ=255)
                                                                                                                                                   ctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATa    >  2:1079650/1‑150 (MQ=255)
                                                                                                                                                     aTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACt  <  2:2311763/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACT  >  CP009273/3577678‑3577976

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