Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,197 G→T 21.0% H174Q (CAC→CAA sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,150,1970GT21.0% 876.3 / 158.1 315H174Q (CAC→CAAsthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (123/125);  new base T (32/35);  total (155/160)
Fisher's exact test for biased strand distribution p-value = 8.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

GAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAGCGCCAA  >  CP009273/4150048‑4150346
                                                                                                                                                     |                                                                                                                                                     
gagagagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCg                                                                                                                                                       >  2:3060539/1‑150 (MQ=255)
  gagagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtg                                                                                                                                                     <  2:278751/150‑1 (MQ=255)
    gagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCttggt                                                                                                                                                   >  1:2044331/1‑150 (MQ=255)
      gaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGc                                                                                                                                                 <  2:2687834/150‑1 (MQ=255)
         tCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATg                                                                                                                                              <  2:3012689/150‑1 (MQ=255)
         tCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATg                                                                                                                                              <  2:1125581/150‑1 (MQ=255)
             aCATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTga                                                                                                                                          >  2:4577688/1‑150 (MQ=255)
                  tctTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAAtt                                                                                                                                     >  1:1777287/1‑150 (MQ=255)
                  tctTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACACCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAtt                                                                                                                                     >  2:1774320/1‑150 (MQ=255)
                   ctTGATCGAGAAATGCCAGCAGGCGATCGCTGGTGTTGATCAGATCCACTTTTACATCCATCCCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCAGAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTg                                                                                                                                    >  2:1315096/1‑150 (MQ=255)
                    ttGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                                                                                                                                   <  1:1241688/150‑1 (MQ=255)
                    ttGATCGAAAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGa                                                                                                                                   >  1:4892703/1‑150 (MQ=255)
                      gATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGt                                                                                                                                 <  2:1651527/150‑1 (MQ=255)
                       aTCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTc                                                                                                                                <  1:2797333/150‑1 (MQ=255)
                           agaAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTg                                                                                                                            >  2:2833632/1‑150 (MQ=255)
                            aaaaTGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGt                                                                                                                           >  2:1254474/2‑150 (MQ=255)
                             aaaTGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTc                                                                                                                          >  1:559076/1‑150 (MQ=255)
                              aaTGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACCTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCg                                                                                                                         <  2:3618495/150‑1 (MQ=255)
                               aTGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGt                                                                                                                        <  1:4996649/150‑1 (MQ=255)
                                   cagcagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAAt                                                                                                                    >  1:2928060/1‑150 (MQ=255)
                                    atcatgCGATCGCGGGTGTTGAAAAGATCCACTTTAAAATCCATAACGCGTAAGATCGACGCATATTCACATCCGATCACTCCAGCACCATAGATAATAACATGGCGCCGTACGTGGAGCATGCTGAGAATTGAGTCGCTGTCTTAAATg                                                                                                                   >  2:2213163/6‑150 (MQ=255)
                                      cagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCg                                                                                                                 <  2:817217/150‑1 (MQ=255)
                                         gCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTgg                                                                                                              >  2:3501469/1‑150 (MQ=255)
                                            aTCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATg                                                                                                           >  1:2957562/1‑150 (MQ=255)
                                               gcgGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGt                                                                                                        <  2:1692469/150‑1 (MQ=255)
                                                cgGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTg                                                                                                       >  2:3820469/1‑150 (MQ=255)
                                                 gggTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGa                                                                                                      >  1:352606/1‑150 (MQ=255)
                                                  ggTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGaa                                                                                                     >  2:1985491/1‑150 (MQ=255)
                                                   gtgtTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGaaa                                                                                                    >  2:668196/1‑150 (MQ=255)
                                                    tgtTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAt                                                                                                   <  1:2728118/150‑1 (MQ=255)
                                                    tgtTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAt                                                                                                   >  1:2604749/1‑150 (MQ=255)
                                                     gtTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATc                                                                                                  <  2:5133829/150‑1 (MQ=255)
                                                     gtTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATc                                                                                                  >  1:556707/1‑150 (MQ=255)
                                                        gATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAAc                                                                                               <  1:3716699/150‑1 (MQ=255)
                                                           cAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATc                                                                                            >  1:3761240/1‑150 (MQ=255)
                                                                ccACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTg                                                                                       <  2:179985/150‑1 (MQ=255)
                                                                ccACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAACAACATCTGTTg                                                                                       >  1:3082511/1‑150 (MQ=255)
                                                                   cTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGTCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTACATGCGTGGATGGGTGAAATCAACATCTGTTTgat                                                                                    >  2:3933775/1‑150 (MQ=255)
                                                                   cTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgat                                                                                    <  1:1198691/150‑1 (MQ=255)
                                                                   cTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgat                                                                                    >  1:617345/1‑150 (MQ=255)
                                                                     tttACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgatga                                                                                  <  2:3704073/150‑1 (MQ=255)
                                                                     tttACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgatga                                                                                  <  2:3704101/150‑1 (MQ=255)
                                                                     tttACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgatga                                                                                  >  2:4228832/1‑150 (MQ=255)
                                                                        aCATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGatat                                                                               >  2:1209150/1‑150 (MQ=255)
                                                                        aCATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGatat                                                                               <  2:5041958/150‑1 (MQ=255)
                                                                         cATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATg                                                                              >  1:3017182/1‑150 (MQ=255)
                                                                         cATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATg                                                                              <  1:2366280/150‑1 (MQ=255)
                                                                          aTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATgg                                                                             <  1:4231816/150‑1 (MQ=255)
                                                                          aTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATgg                                                                             >  1:3995102/1‑150 (MQ=255)
                                                                          aTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATgg                                                                             >  2:4935907/1‑150 (MQ=255)
                                                                          aTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATgg                                                                             <  2:3726624/150‑1 (MQ=255)
                                                                            ccATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGAc                                                                           <  1:1774320/150‑1 (MQ=255)
                                                                             cATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACg                                                                          <  1:4967646/150‑1 (MQ=255)
                                                                                aCCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACgaga                                                                       <  2:5076973/150‑1 (MQ=255)
                                                                                  cgcgGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGc                                                                     >  1:4363212/1‑150 (MQ=255)
                                                                                    cgGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGAGATGGACGAGAGCCg                                                                   >  2:3102968/1‑150 (MQ=255)
                                                                                     ggAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGc                                                                  <  2:545520/150‑1 (MQ=255)
                                                                                     ggAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGc                                                                  >  1:3195444/1‑150 (MQ=255)
                                                                                     ggAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGc                                                                  <  1:955530/150‑1 (MQ=255)
                                                                                       acgATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACCTGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAg                                                                <  2:2338363/148‑1 (MQ=255)
                                                                                           tCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCaa                                                            >  1:2496708/1‑150 (MQ=255)
                                                                                             gaagCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATa                                                          >  2:315256/4‑150 (MQ=255)
                                                                                             gACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATa                                                          >  1:336904/1‑150 (MQ=255)
                                                                                             gACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATa                                                          >  1:1150406/1‑150 (MQ=255)
                                                                                             gACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATa                                                          >  1:3201918/1‑150 (MQ=255)
                                                                                              aCGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATaa                                                         <  1:4764630/150‑1 (MQ=255)
                                                                                               cGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAAc                                                        <  2:3384150/150‑1 (MQ=255)
                                                                                                gCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATTATATGGACGAGAGCCGCAGGCAATAACa                                                       >  1:3854112/1‑150 (MQ=255)
                                                                                                gCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACa                                                       <  1:2476605/150‑1 (MQ=255)
                                                                                                  atatTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACaaa                                                     <  1:5147254/150‑1 (MQ=255)
                                                                                                  atatTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACaaa                                                     >  1:167948/1‑150 (MQ=255)
                                                                                                  atatTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACaaa                                                     <  2:1891246/150‑1 (MQ=255)
                                                                                                  atatTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACaaa                                                     >  1:927215/1‑150 (MQ=255)
                                                                                                               aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGtt                                        >  1:78154/1‑150 (MQ=255)
                                                                                                               aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGtt                                        >  2:2927751/1‑150 (MQ=255)
                                                                                                                tCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTa                                       <  2:207972/150‑1 (MQ=255)
                                                                                                                 cACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAg                                      >  1:2140898/1‑150 (MQ=255)
                                                                                                                 cACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAg                                      <  1:955197/150‑1 (MQ=255)
                                                                                                                  aCTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAgt                                     <  2:4361972/150‑1 (MQ=255)
                                                                                                                   cTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAgtg                                    <  2:1428747/150‑1 (MQ=255)
                                                                                                                       aGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTc                                <  2:2052854/150‑1 (MQ=255)
                                                                                                                       aGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTc                                <  2:4122561/150‑1 (MQ=255)
                                                                                                                       aGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTc                                <  2:1560816/150‑1 (MQ=255)
                                                                                                                        gCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCa                               <  1:1026161/150‑1 (MQ=255)
                                                                                                                         cACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCaa                              >  1:1282907/1‑150 (MQ=255)
                                                                                                                          aCCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAAc                             <  1:4456182/150‑1 (MQ=255)
                                                                                                                          aCCATAGATAAGTACATGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAAc                             <  2:3312644/150‑1 (MQ=255)
                                                                                                                             aTAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCt                          <  1:2198697/150‑1 (MQ=255)
                                                                                                                                      tACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCgg                 <  1:1385438/150‑1 (MQ=255)
                                                                                                                                      tACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCgg                 <  1:4096521/150‑1 (MQ=255)
                                                                                                                                       aCATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCggg                >  1:3625095/1‑150 (MQ=255)
                                                                                                                                         aTGGCGCGGTTCTTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCa              <  1:1939197/150‑1 (MQ=255)
                                                                                                                                         aTGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCa              <  1:504826/150‑1 (MQ=255)
                                                                                                                                            gcgcGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATc           <  2:932025/150‑1 (MQ=255)
                                                                                                                                             cgcgGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATcc          >  2:4890964/1‑150 (MQ=255)
                                                                                                                                               cgGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCACTCCAg        >  1:4605699/1‑150 (MQ=255)
                                                                                                                                               cgGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAg        >  2:4712778/1‑150 (MQ=255)
                                                                                                                                               cgGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAg        >  2:4712720/1‑150 (MQ=255)
                                                                                                                                               cgGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAg        <  2:3247000/150‑1 (MQ=255)
                                                                                                                                                     gtggtgCATGCTGAGAATTGAGTCGCGGTCGTAAATGCGTGGATGGGTGAAATCCACCTCTGTTGGATGATATGGACGAGAGCCGCAGGCAAGAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAGCGCCaa  <  2:2504975/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCGGGCAATCCAGCGCCAA  >  CP009273/4150048‑4150346

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: