Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,091 G→T 8.2% R210S (CGC→AGC)  sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,150,0910GT 8.2% 1184.8 / 52.6 315R210S (CGC→AGC) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (136/153);  new base T (5/21);  total (141/174)
Fisher's exact test for biased strand distribution p-value = 6.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

TTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGAT  >  CP009273/4149945‑4150240
                                                                                                                                                  |                                                                                                                                                     
       gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtg                                                                                                                                             <  2:3053166/150‑1 (MQ=255)
        aCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGAGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCCGCAGGCGATCGCGGgtgt                                                                                                                                            >  2:2506503/1‑150 (MQ=255)
           tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGa                                                                                                                                         <  1:4533040/150‑1 (MQ=255)
              gATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCa                                                                                                                                      <  2:2182070/150‑1 (MQ=255)
              gATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCa                                                                                                                                      <  1:2540956/150‑1 (MQ=255)
                 gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGAt                                                                                                                                   >  2:4361174/1‑150 (MQ=255)
                 gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGAt                                                                                                                                   <  2:3265425/150‑1 (MQ=255)
                   aTGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATcc                                                                                                                                 <  1:2113004/150‑1 (MQ=255)
                    tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                                                                <  1:4408451/150‑1 (MQ=255)
                    tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                                                                <  2:3130761/150‑1 (MQ=255)
                      aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACt                                                                                                                              <  2:2859513/150‑1 (MQ=255)
                      aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACt                                                                                                                              <  2:1543974/150‑1 (MQ=255)
                        cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                                                                                            >  1:2052854/1‑150 (MQ=255)
                         acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                                                           <  2:182056/150‑1 (MQ=255)
                         acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                                                           <  1:3600381/150‑1 (MQ=255)
                          cacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                                                          >  1:1651527/1‑150 (MQ=255)
                              aTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                                                                      <  1:4355605/150‑1 (MQ=255)
                              aTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                                                                      >  2:2994500/1‑150 (MQ=255)
                                cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCa                                                                                                                    >  1:2345989/1‑150 (MQ=255)
                                 gTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                                                                                   <  1:4139829/150‑1 (MQ=255)
                                  tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                                                                  <  2:2625160/150‑1 (MQ=255)
                                  tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                                                                  >  2:2060058/1‑150 (MQ=255)
                                  tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                                                                  <  2:961494/150‑1 (MQ=255)
                                   cacaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                                                                                 >  2:3016728/1‑150 (MQ=255)
                                   cacaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                                                                                 >  2:4764684/1‑150 (MQ=255)
                                        ccTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCgg                                                                                                            <  2:4650004/150‑1 (MQ=255)
                                         cTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                                                           <  2:4092489/150‑1 (MQ=255)
                                          ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGaa                                                                                                          <  2:1266561/150‑1 (MQ=255)
                                            cGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                                                                        <  1:3455288/150‑1 (MQ=255)
                                            cGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                                                                        <  1:1445907/150‑1 (MQ=255)
                                            cGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                                                                        <  1:1935589/150‑1 (MQ=255)
                                              aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                                                                                      <  1:3497994/150‑1 (MQ=255)
                                               tcttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                                                                                                     >  1:5076973/1‑150 (MQ=255)
                                                    tcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCa                                                                                                >  2:1228741/1‑150 (MQ=255)
                                                      gTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                                                                              >  2:4949669/1‑150 (MQ=255)
                                                      gTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                                                                              <  1:3672782/150‑1 (MQ=255)
                                                        aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATAtt                                                                                            >  2:737516/1‑150 (MQ=255)
                                                         ctctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTc                                                                                           <  2:1164989/150‑1 (MQ=255)
                                                             tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGATATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                                                                       <  1:1245392/150‑1 (MQ=255)
                                                             tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATCCCGCGGAAGATCGACGCATATTCACAg                                                                                       >  2:1867151/1‑150 (MQ=255)
                                                             tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                                                                       >  1:3977329/1‑150 (MQ=255)
                                                              cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                                                                      <  1:3199405/150‑1 (MQ=255)
                                                              cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                                                                      >  2:3788203/1‑150 (MQ=255)
                                                              cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                                                                      <  1:3039107/150‑1 (MQ=255)
                                                              cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                                                                      <  2:312802/150‑1 (MQ=255)
                                                               gTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                                                                     <  1:86415/150‑1 (MQ=255)
                                                               gTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                                                                     >  2:1420605/1‑150 (MQ=255)
                                                                ttGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCg                                                                                    >  2:4480714/1‑150 (MQ=255)
                                                                 tgtgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCTGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa                                                                                   >  2:955512/1‑150 (MQ=255)
                                                                   tgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                                                                 >  1:108718/1‑150 (MQ=255)
                                                                      cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCTGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACt                                                                              >  2:2728118/1‑150 (MQ=255)
                                                                        aaTCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                                                            <  1:2690637/150‑1 (MQ=255)
                                                                        aaTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                                                            <  2:1053103/150‑1 (MQ=255)
                                                                             cTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAc                                                                       <  1:4730039/150‑1 (MQ=255)
                                                                             cTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAc                                                                       <  1:243144/150‑1 (MQ=255)
                                                                               aCGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCa                                                                     >  1:2142216/1‑150 (MQ=255)
                                                                                  ccACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAg                                                                  >  2:1939151/1‑150 (MQ=255)
                                                                                   cNCTGTTCCAGAAGTGATAGGAGAGAGACTCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATCCCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAccctcga                                                                 >  1:2504975/1‑145 (MQ=255)
                                                                                    aCTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAt                                                                >  2:504826/1‑150 (MQ=255)
                                                                                      tGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                                                              >  2:1111998/1‑150 (MQ=255)
                                                                                        ttCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                                                            <  2:2475290/150‑1 (MQ=255)
                                                                                         tCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTa                                                           >  2:1642097/1‑150 (MQ=255)
                                                                                           cAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa                                                         <  1:3247196/150‑1 (MQ=255)
                                                                                            aGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACAt                                                        <  2:3227946/150‑1 (MQ=255)
                                                                                             gAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATg                                                       <  2:1732889/150‑1 (MQ=255)
                                                                                             gAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATg                                                       >  2:600583/1‑150 (MQ=255)
                                                                                              aaGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATgg                                                      <  1:1327056/150‑1 (MQ=255)
                                                                                               aGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgc                                                     <  2:2336683/150‑1 (MQ=255)
                                                                                                gTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcg                                                    >  1:4213511/1‑150 (MQ=255)
                                                                                                 tGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcgc                                                   <  1:1399337/150‑1 (MQ=255)
                                                                                                    tAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt                                                >  2:1323762/1‑150 (MQ=255)
                                                                                                        agagagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTctt                                            <  2:1990988/150‑3 (MQ=255)
                                                                                                         gagagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtg                                           <  2:278751/150‑1 (MQ=255)
                                                                                                           gagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCttggt                                         >  1:2044331/1‑150 (MQ=255)
                                                                                                           gagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtggt                                         >  2:4042564/1‑150 (MQ=255)
                                                                                                                 cTGACATCTCTTGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGc                                   >  2:360595/1‑150 (MQ=255)
                                                                                                                  tGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCt                                  >  2:3716764/1‑150 (MQ=255)
                                                                                                                       tctcTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGaa                             <  2:3919885/150‑1 (MQ=255)
                                                                                                                        ctctTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAAt                            >  1:3788512/1‑150 (MQ=255)
                                                                                                                         tctTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAAtt                           >  1:1777287/1‑150 (MQ=255)
                                                                                                                         tctTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAtt                           >  2:1777494/1‑150 (MQ=255)
                                                                                                                           ttGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                         <  1:1283673/150‑1 (MQ=255)
                                                                                                                           ttGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                         <  1:4502967/150‑1 (MQ=255)
                                                                                                                           ttGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                         >  1:4892703/1‑150 (MQ=255)
                                                                                                                           ttGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAAAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                         >  1:5122058/1‑150 (MQ=255)
                                                                                                                             gATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGt                       >  2:528909/1‑150 (MQ=255)
                                                                                                                             gATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGt                       <  1:2651381/150‑1 (MQ=255)
                                                                                                                              aTCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTc                      <  1:2797659/150‑1 (MQ=255)
                                                                                                                               tCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCg                     <  2:2939560/150‑1 (MQ=255)
                                                                                                                                cGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATCCCGCGGAAGATCGACGCATAGTCACAGCCGCTCACTCCAGCACCCTCGCTCAGTACATGGCGCGGTGCGGGGGGCCTGCTGAGACTTGAGTCGc                    >  2:4399689/1‑150 (MQ=255)
                                                                                                                                   aaaaTGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGACTTGAGTCGCTGt                 >  1:4543792/2‑150 (MQ=255)
                                                                                                                                    aaaTGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTc                >  1:559076/1‑150 (MQ=255)
                                                                                                                                      aTGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGt              <  1:4996256/150‑1 (MQ=255)
                                                                                                                                         ccAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTaaa           >  2:350050/1‑150 (MQ=255)
                                                                                                                                          cagcagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAAt          <  2:4169309/150‑1 (MQ=255)
                                                                                                                                          cagcagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAAt          >  1:2928060/1‑150 (MQ=255)
                                                                                                                                           atcatgCGATCGCGGGTGTTGAAAAGATCCACTTTAAAATCCATAACGCGTAAGATCGACGCATATTCACATCCGATCACTCCAGCACCATAGATAATAACATGGCGCCGTACGTGGAGCATGCTGAGAATTGAGTCGCTGTCTTAAATg         >  2:2213163/6‑150 (MQ=255)
                                                                                                                                           agcagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAGGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATg         >  1:3311728/1‑150 (MQ=255)
                                                                                                                                                 cGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGa   >  1:757733/1‑150 (MQ=255)
                                                                                                                                                  gATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGAt  <  2:1152047/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGAT  >  CP009273/4149945‑4150240

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: