Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,428,092 G→T 92.7% G190V (GGA→GTA)  ytfI → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,428,0920GT92.7% 812.5 / 55.6 285G190V (GGA→GTA) ytfIuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (8/13);  new base T (145/119);  total (153/132)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 289 total aligned reads displayed.

ATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAG  >  CP009273/4427945‑4428240
                                                                                                                                                   |                                                                                                                                                    
atatTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCagtag                                                                                                                                                    >  1:1867822/1‑150 (MQ=255)
   ttCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTg                                                                                                                                                 >  2:618819/1‑150 (MQ=255)
       aCAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAAt                                                                                                                                             >  1:4584753/1‑150 (MQ=255)
       aCAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAAt                                                                                                                                             >  1:2534009/1‑150 (MQ=255)
         aGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTa                                                                                                                                           >  2:4316720/1‑150 (MQ=255)
            aaaaaCAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAg                                                                                                                                        <  2:4253583/150‑1 (MQ=255)
            aaaaaCAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAg                                                                                                                                        <  2:2764273/150‑1 (MQ=255)
            aaaaaCAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAg                                                                                                                                        >  1:2686911/1‑150 (MQ=255)
             aaaaCAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGa                                                                                                                                       <  1:3857337/150‑1 (MQ=255)
                 cAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGc                                                                                                                                   <  1:1422139/150‑1 (MQ=255)
                    tCTGACATACCGAAAATGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgat                                                                                                                                <  2:3510668/150‑1 (MQ=255)
                    tCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgat                                                                                                                                <  2:1192205/150‑1 (MQ=255)
                     cTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATTCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatgg                                                                                                                              >  1:167528/1‑150 (MQ=255)
                     cTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatg                                                                                                                               <  2:3751148/150‑1 (MQ=255)
                     cTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatg                                                                                                                               >  1:4152049/1‑150 (MQ=255)
                      tGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatgg                                                                                                                              >  2:4559136/1‑150 (MQ=255)
                        aCATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatggat                                                                                                                            <  2:1888990/150‑1 (MQ=255)
                         cATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatggatg                                                                                                                           >  2:2378699/1‑150 (MQ=255)
                         cATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCgatggatg                                                                                                                           >  1:3665547/1‑150 (MQ=255)
                          aTACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGatggatgg                                                                                                                          >  2:1034543/1‑150 (MQ=255)
                          aTACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGatggatgg                                                                                                                          >  1:4133376/1‑150 (MQ=255)
                          aTACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGatggatgg                                                                                                                          >  1:3174912/1‑150 (MQ=255)
                           tACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGt                                                                                                                         <  2:4530854/150‑1 (MQ=255)
                             ccGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGtat                                                                                                                       >  1:737240/1‑150 (MQ=255)
                             ccGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGtat                                                                                                                       <  1:3600406/150‑1 (MQ=255)
                             ccGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGtat                                                                                                                       <  2:2307769/150‑1 (MQ=255)
                             ccGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGtat                                                                                                                       <  1:1935298/150‑1 (MQ=255)
                              cGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTACGTGGAATTAAAGAATGCGATGGATGGtata                                                                                                                      >  1:3035700/1‑150 (MQ=255)
                                 aaaGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATAttt                                                                                                                   >  1:657917/1‑150 (MQ=255)
                                 aaaGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATAttt                                                                                                                   >  2:3840735/1‑150 (MQ=255)
                                 aaaGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATAttt                                                                                                                   >  1:1226991/1‑150 (MQ=255)
                                   aGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTAc                                                                                                                 >  1:4524810/1‑150 (MQ=255)
                                         aaaGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATTTTTACGaaaaa                                                                                                           <  1:3908452/150‑1 (MQ=255)
                                         aaaGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                                                                                                           <  1:2900022/150‑1 (MQ=255)
                                         aaaGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGaaaaa                                                                                                           >  1:4373974/1‑150 (MQ=255)
                                           aGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaat                                                                                                         <  1:1096028/150‑1 (MQ=255)
                                              ttATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAACACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaataat                                                                                                      <  1:4217147/150‑1 (MQ=255)
                                              ttATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAaataat                                                                                                      >  2:3480730/1‑150 (MQ=255)
                                               tataTTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAAtt                                                                                                     <  1:3026019/150‑1 (MQ=255)
                                               tataTGCACCCGCCGAGTGGTGACGTGTCTTTAGGGTATGCCTTATTGGAGGATATAGTAAGCAAATTCAAACGGGGTAAATTGGAACTCCAAATACCAGTAGGTGGAATTACAGAATGCGATGGATGGTATATTTACGAAAAAATAAtt                                                                                                     <  2:2419669/150‑1 (MQ=255)
                                                atatTCACCCGGCGATTTGTTACTTTTCTTTATGTTATGCATTATTTGAGTATATATTAATCAAATTAAAACGTGGTAAATTTTATATCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGATAAAATTATtg                                                                                                    <  2:846824/150‑1 (MQ=255)
                                                atatTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtg                                                                                                    <  2:176593/150‑1 (MQ=255)
                                                 tatTCACCAGGCGATTGGTTACGTGTCTTTAGGTTATGCCTTATTGGACGATATAGTAATCAAATTAAAACGGGGTAAATTTGAACTCAAAATACAAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtga                                                                                                   <  2:1006234/150‑1 (MQ=255)
                                                  atTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtgat                                                                                                  <  1:3433212/150‑1 (MQ=255)
                                                   ttCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtgatg                                                                                                 <  1:2771938/150‑1 (MQ=255)
                                                   ttCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACACCAAATACCAGTAGGTGGAATTAAATAATGCGATGGATGGTATATTTACTAAAAAATAATtgatg                                                                                                 >  2:1586800/1‑150 (MQ=255)
                                                    tCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATtgatga                                                                                                <  2:3604732/150‑1 (MQ=255)
                                                     cACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGTATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTgatgat                                                                                               >  1:1656002/1‑150 (MQ=255)
                                                     cACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTgatgat                                                                                               >  2:661470/1‑150 (MQ=255)
                                                           gCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAAtttt                                                                                         <  2:111050/150‑1 (MQ=255)
                                                               ttGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCa                                                                                     <  1:3793688/150‑1 (MQ=255)
                                                                  gTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAg                                                                                  <  2:2650873/150‑1 (MQ=255)
                                                                    gACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTc                                                                                >  1:4034681/1‑150 (MQ=255)
                                                                     aCGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCa                                                                               <  2:1843697/150‑1 (MQ=255)
                                                                       gtgtCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATa                                                                             <  2:4076228/150‑1 (MQ=255)
                                                                       gtgtCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATa                                                                             <  1:2457070/150‑1 (MQ=255)
                                                                           cTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAAt                                                                         >  2:3191964/1‑150 (MQ=255)
                                                                               aGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCttt                                                                     >  2:3281106/1‑150 (MQ=255)
                                                                                 gTTATGCCTTATTGTAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTgg                                                                   <  2:3921736/150‑1 (MQ=255)
                                                                                 gTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAACTTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTgg                                                                   <  2:1559613/150‑1 (MQ=255)
                                                                                 gTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTgg                                                                   >  1:1351578/1‑150 (MQ=255)
                                                                                  ttATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTggg                                                                  <  2:1783648/150‑1 (MQ=255)
                                                                                             tGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTa                                                       >  2:3392534/1‑150 (MQ=255)
                                                                                               gaggaTATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTaaa                                                     <  2:3170070/150‑1 (MQ=255)
                                                                                               gaggaTATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTaaa                                                     <  2:497641/150‑1 (MQ=255)
                                                                                                  gaTATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGCTGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATc                                                  <  2:1506508/150‑1 (MQ=255)
                                                                                                   atatAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCt                                                 <  1:1563410/150‑1 (MQ=255)
                                                                                                   atatAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCt                                                 >  2:584799/1‑150 (MQ=255)
                                                                                                      taGTAAGCAAATTAACACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATg                                              <  1:4471420/150‑1 (MQ=255)
                                                                                                        gTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGa                                            >  1:2298177/1‑150 (MQ=255)
                                                                                                        gTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGa                                            >  2:2968648/1‑150 (MQ=255)
                                                                                                            gCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAgggg                                        >  2:1642914/1‑150 (MQ=255)
                                                                                                            gCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAgggg                                        >  1:1319423/1‑150 (MQ=255)
                                                                                                              aaaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGAt                                      <  1:2034253/150‑1 (MQ=255)
                                                                                                               aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATg                                     >  1:2211025/1‑150 (MQ=255)
                                                                                                               aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATg                                     <  1:4454729/150‑1 (MQ=255)
                                                                                                               aaTTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATg                                     >  1:3745488/1‑150 (MQ=255)
                                                                                                                aTTAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGa                                    >  1:883585/1‑150 (MQ=255)
                                                                                                                  tAAAACGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAAc                                  >  1:3292351/1‑150 (MQ=255)
                                                                                                                   aaaaCGGGGTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACg                                 >  1:4019810/1‑150 (MQ=255)
                                                                                                                       cGGGGTAAATTTTAACTAAAAATACCAGTAGGTGGAATTAAAGAATGCTATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGtttt                             <  2:3038408/150‑1 (MQ=255)
                                                                                                                        ggggTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTc                            <  1:1891073/150‑1 (MQ=255)
                                                                                                                        ggggTAAATTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTc                            <  2:401332/150‑1 (MQ=255)
                                                                                                                            tAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTaa                        <  1:2501049/150‑1 (MQ=255)
                                                                                                                             aaaTTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAAt                       <  2:3214813/150‑1 (MQ=255)
                                                                                                                               aTTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATgg                     >  2:4503978/1‑150 (MQ=255)
                                                                                                                               aTTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATgg                     <  1:4271925/150‑1 (MQ=255)
                                                                                                                               aTTTGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATgg                     >  2:3063053/1‑150 (MQ=255)
                                                                                                                                  tGAACTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTc                  <  2:1747107/150‑1 (MQ=255)
                                                                                                                                    aaCTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCct                >  2:4174213/1‑150 (MQ=255)
                                                                                                                                      cTCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCctct              >  2:1869022/1‑150 (MQ=255)
                                                                                                                                       tCCAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCtg             >  1:617281/1‑150 (MQ=255)
                                                                                                                                        ccAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCtgt            <  2:4299897/150‑1 (MQ=255)
                                                                                                                                         cAAATACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCtgtt           <  1:3463903/150‑1 (MQ=255)
                                                                                                                                          aaaTGCCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCtgttg          >  2:3638753/1‑150 (MQ=255)
                                                                                                                                          aaaTACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCtgttg          <  2:1336934/150‑1 (MQ=255)
                                                                                                                                            aTACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTgttgtt        <  2:1402926/150‑1 (MQ=255)
                                                                                                                                             tACCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGTGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGttt       <  1:967374/150‑1 (MQ=255)
                                                                                                                                              aCCAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGtttt      >  2:2090332/1‑150 (MQ=255)
                                                                                                                                                cAGTAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTgg    >  1:1775005/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAG  >  CP009273/4427945‑4428240

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: