Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,654,383 A→T 9.4% C170S (TGC→AGC)  iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,654,3830AT 9.4% 745.3 / 51.1 229C170S (TGC→AGC) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (105/102);  new base T (13/9);  total (118/111)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

ACGTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTT  >  CP009273/2654234‑2654532
                                                                                                                                                     |                                                                                                                                                     
acgTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGca                                                                                                                                                       >  2:2763323/1‑150 (MQ=255)
  gTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGcaca                                                                                                                                                     >  2:729199/1‑150 (MQ=255)
    acatacaGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACAtt                                                                                                                                                   <  1:756171/150‑1 (MQ=255)
     catacaGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACAttt                                                                                                                                                  <  2:2674990/150‑1 (MQ=255)
     catacaGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACAttt                                                                                                                                                  <  2:2450257/150‑1 (MQ=255)
       tacaGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCGGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCg                                                                                                                                                >  1:1261467/1‑150 (MQ=255)
        acaGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATCATGCCACGAGCACGGCACATTTCGc                                                                                                                                               >  2:1982193/1‑150 (MQ=255)
               aCCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGc                                                                                                                                        <  2:677050/150‑1 (MQ=255)
                ccGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGcc                                                                                                                                       <  1:2580313/150‑1 (MQ=255)
                 cGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCcg                                                                                                                                      >  2:1409035/1‑150 (MQ=255)
                   aTACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCcgcg                                                                                                                                    <  2:1776839/150‑1 (MQ=255)
                   aTACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCcgcg                                                                                                                                    >  2:3764424/1‑150 (MQ=255)
                    tACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGa                                                                                                                                   >  1:1140848/1‑150 (MQ=255)
                    tACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGa                                                                                                                                   >  2:660274/1‑150 (MQ=255)
                    tACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGa                                                                                                                                   <  2:1212496/150‑1 (MQ=255)
                    tACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTATACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGACACGAGCACGGCACATTTCGCCGAAAGCCGCGa                                                                                                                                   >  2:906238/1‑150 (MQ=255)
                       cTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGatat                                                                                                                                >  1:2910666/1‑150 (MQ=255)
                          tCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCt                                                                                                                             <  1:3228972/150‑1 (MQ=255)
                                cATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGcacca                                                                                                                       >  2:2071751/1‑150 (MQ=255)
                                cATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGcacca                                                                                                                       >  2:3440854/1‑150 (MQ=255)
                                    gATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGcc                                                                                                                   >  2:2019912/1‑150 (MQ=255)
                                     aatatctgACCGGACAAAGACATCATGTCAACATTCAACTGGCTCACCAATAAAGTCATTAAACCCACTCTAAGGGTACAATAAAACTGATATATAAACACACGAGCACGGCACATTTCGCCGATAGCCCCGATATCCATCACAACGcct                                                                                                                  >  2:1092600/7‑149 (MQ=255)
                                     aTTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCg                                                                                                                  <  2:758452/150‑1 (MQ=255)
                                     aTTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCg                                                                                                                  >  2:1092645/1‑150 (MQ=255)
                                     aTTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCg                                                                                                                  <  2:1319703/150‑1 (MQ=255)
                                      atttGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGa                                                                                                                 >  2:478039/2‑150 (MQ=255)
                                        ttGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGAtt                                                                                                               <  1:373109/150‑1 (MQ=255)
                                         tgtgACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGAttt                                                                                                              <  2:2486365/150‑1 (MQ=255)
                                          gtgACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTc                                                                                                             <  1:3513841/150‑1 (MQ=255)
                                            gACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGt                                                                                                           >  1:1829978/1‑150 (MQ=255)
                                               cGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTAt                                                                                                        >  2:3854203/1‑150 (MQ=255)
                                                 gagaAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTAttt                                                                                                      >  1:1851521/1‑150 (MQ=255)
                                                  agaAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTa                                                                                                     >  1:1292851/1‑150 (MQ=255)
                                                    aaaGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACg                                                                                                   <  2:2227430/150‑1 (MQ=255)
                                                     aaGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGt                                                                                                  >  1:4026998/1‑150 (MQ=255)
                                                     aaGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGt                                                                                                  <  2:3092689/150‑1 (MQ=255)
                                                           tCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCatga                                                                                            >  2:3432382/1‑150 (MQ=255)
                                                            cAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCatgat                                                                                           >  1:1505808/1‑150 (MQ=255)
                                                               gTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGa                                                                                        <  2:472942/150‑1 (MQ=255)
                                                                tCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGac                                                                                       >  2:3094674/1‑150 (MQ=255)
                                                                tCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGac                                                                                       >  2:2231207/1‑150 (MQ=255)
                                                                 cAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGaca                                                                                      <  1:4076699/150‑1 (MQ=255)
                                                                  aaCTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGacac                                                                                     >  1:3246443/1‑150 (MQ=255)
                                                                      ttCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAgg                                                                                 >  2:508915/1‑150 (MQ=255)
                                                                      ttCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAgg                                                                                 >  2:878060/1‑150 (MQ=255)
                                                                       tCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGa                                                                                <  1:2293636/150‑1 (MQ=255)
                                                                         aaCTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATg                                                                              <  1:2796198/150‑1 (MQ=255)
                                                                         aaCTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATg                                                                              <  1:3350360/150‑1 (MQ=255)
                                                                          aCTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATgg                                                                             >  1:797072/1‑150 (MQ=255)
                                                                             ggCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGgtgt                                                                          <  1:2198628/150‑1 (MQ=255)
                                                                             ggCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGgtgt                                                                          >  1:2640410/1‑150 (MQ=255)
                                                                              gCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtc                                                                         >  1:1054728/1‑150 (MQ=255)
                                                                              gCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtc                                                                         >  1:2784190/1‑150 (MQ=255)
                                                                              gCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtc                                                                         >  1:390553/1‑150 (MQ=255)
                                                                               cTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtcg                                                                        >  2:1227678/1‑150 (MQ=255)
                                                                                tCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtcgt                                                                       >  1:2475793/1‑150 (MQ=255)
                                                                                 cAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtcgtc                                                                      >  1:1656953/1‑150 (MQ=255)
                                                                                 cAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTgtcgtc                                                                      <  2:3700964/150‑1 (MQ=255)
                                                                                    gTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACg                                                                   <  1:2113737/150‑1 (MQ=255)
                                                                                     tCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGc                                                                  >  1:1479964/1‑150 (MQ=255)
                                                                                      cGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCa                                                                 <  1:765511/150‑1 (MQ=255)
                                                                                       gATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCAt                                                                >  1:3588152/1‑150 (MQ=255)
                                                                                       gATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCAt                                                                >  1:901303/1‑150 (MQ=255)
                                                                                        aTAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATc                                                               <  1:4011252/150‑1 (MQ=255)
                                                                                        aTAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATc                                                               <  1:629321/150‑1 (MQ=255)
                                                                                         tAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCg                                                              >  1:229953/1‑150 (MQ=255)
                                                                                           ggCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCgct                                                            <  2:1004340/150‑1 (MQ=255)
                                                                                           ggCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCgct                                                            <  2:3737713/150‑1 (MQ=255)
                                                                                            gCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCgctg                                                           <  2:788348/150‑1 (MQ=255)
                                                                                            gCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCgctg                                                           <  2:1208035/150‑1 (MQ=255)
                                                                                            gCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCgctg                                                           <  1:547747/150‑1 (MQ=255)
                                                                                                  tACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCaa                                                     <  2:2825903/150‑1 (MQ=255)
                                                                                                   aCCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAg                                                    >  2:139327/1‑150 (MQ=255)
                                                                                                     ccACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGtt                                                  >  2:3711269/1‑150 (MQ=255)
                                                                                                       aCGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCt                                                <  1:682439/150‑1 (MQ=255)
                                                                                                        cGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCtt                                               >  1:3058093/1‑150 (MQ=255)
                                                                                                          ctctGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTc                                             >  1:945929/1‑150 (MQ=255)
                                                                                                             tGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAgg                                          >  1:3351412/1‑150 (MQ=255)
                                                                                                              gggTTGCATCAACGTGATAGATAATGCCACGAGCACGGCTCATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGt                                         >  2:2054697/1‑150 (MQ=255)
                                                                                                               ggTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTc                                        <  1:702512/150‑1 (MQ=255)
                                                                                                                    cATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATaa                                   <  1:2906617/150‑1 (MQ=255)
                                                                                                                        aaCGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGcc                               <  2:3503479/150‑1 (MQ=255)
                                                                                                                        aaCGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGcc                               <  2:3544663/150‑1 (MQ=255)
                                                                                                                        aaCGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGcc                               <  2:296223/150‑1 (MQ=255)
                                                                                                                             gatagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTAc                          <  2:4164694/150‑1 (MQ=255)
                                                                                                                              atagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACg                         >  1:283709/1‑150 (MQ=255)
                                                                                                                              atagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACg                         <  2:3815090/150‑1 (MQ=255)
                                                                                                                                agataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCt                       <  2:187222/150‑1 (MQ=255)
                                                                                                                                 gataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTg                      <  2:403959/150‑1 (MQ=255)
                                                                                                                                 gataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTg                      <  2:403971/150‑1 (MQ=255)
                                                                                                                                    aaTGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCgg                   >  2:3879569/1‑150 (MQ=255)
                                                                                                                                      tGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTg                 <  2:1612234/150‑1 (MQ=255)
                                                                                                                                         cACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCa              >  2:2320674/1‑150 (MQ=255)
                                                                                                                                          aCGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCag             <  1:3036588/150‑1 (MQ=255)
                                                                                                                                           cGAGCACGGCACATTTCGCAGATAGCCGCGATATCCTGAACCACGCCGATTTCGTTATTTACGTGCTTGATGGACACGATGATGGTGTCGTCACGCATCGATGCTTCAAGTTCTTTCAGGTCGATAATTCCTTTACGCTGCGGTGCCagg            <  2:1658605/150‑1 (MQ=255)
                                                                                                                                            gAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggt           <  1:999379/150‑1 (MQ=255)
                                                                                                                                             aGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggta          <  1:4132109/150‑1 (MQ=255)
                                                                                                                                               cACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtagg        >  2:3296996/1‑150 (MQ=255)
                                                                                                                                                 cGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTg      <  1:2049446/150‑1 (MQ=255)
                                                                                                                                                     acaTTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTAATTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTGTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACtt  >  1:3427733/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGACCGGAGAAAGACATCAGGTCAACTTTCAACTGGCTCAGGTCGATAGGCAGTTTACCCACGCTCTGGGTTGCATCAACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTT  >  CP009273/2654234‑2654532

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: