Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,754,950 T→A 92.2% I367F (ATT→TTT)  sufD ← component of SufBCD Fe‑S cluster assembly scaffold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,754,9500TA92.2% 804.9 / 45.2 257I367F (ATT→TTT) sufDcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base A (131/106);  total (141/116)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 257 total aligned reads displayed.

ACCGATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGT  >  CP009273/1754804‑1755096
                                                                                                                                                  |                                                                                                                                                  
   gATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGccc                                                                                                                                              >  2:908809/1‑150 (MQ=255)
    aTTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGcccc                                                                                                                                             <  2:2097715/150‑1 (MQ=255)
    aTTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGcccc                                                                                                                                             <  2:635322/150‑1 (MQ=255)
      tCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCAc                                                                                                                                           <  2:2563484/150‑1 (MQ=255)
      tCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCAc                                                                                                                                           >  2:264333/1‑150 (MQ=255)
       cGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCAcc                                                                                                                                          >  1:1964587/1‑150 (MQ=255)
        gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCg                                                                                                                                         >  1:1073522/1‑150 (MQ=255)
        gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCg                                                                                                                                         >  1:4495763/1‑150 (MQ=255)
        gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCg                                                                                                                                         >  1:780308/1‑150 (MQ=255)
         ggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGt                                                                                                                                        <  1:2551590/150‑1 (MQ=255)
             aGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTcgcg                                                                                                                                    <  2:435278/150‑1 (MQ=255)
                aactgctgTTTAAGCCCCTCATAACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCTTCTTCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCg                                                                                                                                 <  1:4540111/148‑1 (MQ=255)
                aCCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCg                                                                                                                                 <  1:1511294/150‑1 (MQ=255)
                  ctgctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTg                                                                                                                               <  1:2091755/150‑1 (MQ=255)
                   tgctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTgg                                                                                                                              >  2:2152902/1‑150 (MQ=255)
                   tgcggTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTgg                                                                                                                              >  2:2189841/1‑150 (MQ=255)
                        tttAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCa                                                                                                                         >  2:3680142/1‑150 (MQ=255)
                         ttAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAt                                                                                                                        <  1:2230190/150‑1 (MQ=255)
                         ttAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAt                                                                                                                        >  1:1813371/1‑150 (MQ=255)
                           aaGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAttt                                                                                                                      >  2:2655012/1‑150 (MQ=255)
                           aaGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAttt                                                                                                                      >  1:3012196/1‑150 (MQ=255)
                             gCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTca                                                                                                                    >  2:2901210/1‑150 (MQ=255)
                                 cTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCAcatc                                                                                                                >  2:1135373/1‑150 (MQ=255)
                                 cTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCAcatc                                                                                                                >  2:1911625/1‑150 (MQ=255)
                                    atcaCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACatcatc                                                                                                             >  2:3964062/1‑150 (MQ=255)
                                      caCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTg                                                                                                           >  2:435532/1‑150 (MQ=255)
                                       aCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGc                                                                                                          <  2:2344340/150‑1 (MQ=255)
                                       aCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGc                                                                                                          <  2:2344419/150‑1 (MQ=255)
                                        cGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCa                                                                                                         >  2:3061842/1‑150 (MQ=255)
                                         gCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCAt                                                                                                        <  1:714601/150‑1 (MQ=255)
                                         gCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCAt                                                                                                        <  1:665227/150‑1 (MQ=255)
                                                ttCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTc                                                                                                 >  2:928465/1‑150 (MQ=255)
                                                ttCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTc                                                                                                 >  2:4311300/1‑150 (MQ=255)
                                                 tCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTcc                                                                                                <  1:3989057/150‑1 (MQ=255)
                                                 tCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTcc                                                                                                <  1:3989140/150‑1 (MQ=255)
                                                  ccGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCa                                                                                               >  2:259085/1‑150 (MQ=255)
                                                    gTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGc                                                                                             <  2:307355/150‑1 (MQ=255)
                                                    gTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGc                                                                                             >  2:1777348/1‑150 (MQ=255)
                                                     tCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCt                                                                                            >  2:457444/1‑150 (MQ=255)
                                                     tCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGAGTTCCAGCt                                                                                            >  2:3765905/1‑150 (MQ=255)
                                                      cAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTCTTCATCATCAAAACGCCCCACCGACGCGCCGTGGCTGCATTTCACAAAATCTGCATAGATTTCCAGcag                                                                                           >  1:4387520/1‑148 (MQ=255)
                                                       aGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGc                                                                                          >  2:4002474/1‑150 (MQ=255)
                                                       aGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGc                                                                                          <  2:3687941/150‑1 (MQ=255)
                                                        gTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCg                                                                                         <  2:2830029/150‑1 (MQ=255)
                                                             gCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTc                                                                                    >  1:1603196/1‑150 (MQ=255)
                                                             gCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTc                                                                                    <  2:1689678/150‑1 (MQ=255)
                                                             gCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTc                                                                                    >  2:1337452/1‑150 (MQ=255)
                                                               aGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGt                                                                                  <  1:2909315/150‑1 (MQ=255)
                                                               aGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGAATGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGATGCGGTTTCGt                                                                                  >  1:124626/1‑150 (MQ=255)
                                                                 cGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGTTTTCGTAt                                                                                >  2:4413963/1‑150 (MQ=255)
                                                                 cGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTAt                                                                                <  2:2243809/150‑1 (MQ=255)
                                                                   aaGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATcc                                                                              >  2:799930/1‑150 (MQ=255)
                                                                   aaGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATcc                                                                              <  1:1536241/150‑1 (MQ=255)
                                                                     ggCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCAc                                                                            >  1:3318221/1‑150 (MQ=255)
                                                                      gCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACt                                                                           <  1:3531560/150‑1 (MQ=255)
                                                                       cGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCTGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACtt                                                                          <  2:2714527/150‑1 (MQ=255)
                                                                       cGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACtt                                                                          <  2:782178/150‑1 (MQ=255)
                                                                        gTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTc                                                                         <  1:1967359/150‑1 (MQ=255)
                                                                         tAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTcc                                                                        <  2:599793/150‑1 (MQ=255)
                                                                         tAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTcc                                                                        <  2:608410/150‑1 (MQ=255)
                                                                         tAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTcc                                                                        >  1:1000574/1‑150 (MQ=255)
                                                                          aaaTGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCg                                                                       <  2:4168003/150‑1 (MQ=255)
                                                                                 catcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGttt                                                                >  1:2032941/1‑150 (MQ=255)
                                                                                 catcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGttt                                                                >  1:2032856/1‑150 (MQ=255)
                                                                                    ctgctgGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGcc                                                             <  1:1775399/150‑1 (MQ=255)
                                                                                     tgctgGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGccc                                                            <  2:1120583/150‑1 (MQ=255)
                                                                                       ctgGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCAt                                                          <  1:384649/150‑1 (MQ=255)
                                                                                       ctgGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCAt                                                          >  2:3901604/1‑150 (MQ=255)
                                                                                           gCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATcagc                                                      >  1:2378690/1‑150 (MQ=255)
                                                                                                ctgctgATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAtt                                                 >  1:896643/1‑150 (MQ=255)
                                                                                                 tgctgATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttg                                                <  1:3598621/150‑1 (MQ=255)
                                                                                                   ctgATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttgtt                                              <  2:4453110/150‑1 (MQ=255)
                                                                                                    tgattgatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttgttg                                             >  2:810965/1‑150 (MQ=255)
                                                                                                      attgatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttgttgtt                                           >  2:2417307/1‑150 (MQ=255)
                                                                                                         gatCCCGCGCGAGCGCAGATTGAATATCTGGTCAAGATCTAATCGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGTTTTCAAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGt                                        <  2:2297231/150‑1 (MQ=255)
                                                                                                         gatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGt                                        >  2:2299291/1‑150 (MQ=255)
                                                                                                         gatCCCGCGCGAGCGCAGATAGAATATCTGATCATCATCAAAACTCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCAAAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAATGTAGTTGGt                                        >  1:603132/1‑150 (MQ=255)
                                                                                                           tCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCa                                      >  1:3226642/1‑150 (MQ=255)
                                                                                                           tCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCa                                      <  1:4179933/150‑1 (MQ=255)
                                                                                                              cgcgcgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCt                                   <  1:853432/150‑1 (MQ=255)
                                                                                                                cgcgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGa                                 <  2:1855533/150‑1 (MQ=255)
                                                                                                                  cgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGAcc                               <  1:2670101/150‑1 (MQ=255)
                                                                                                                   gAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCa                              >  1:2722645/1‑150 (MQ=255)
                                                                                                                    aGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCAt                             >  1:4259812/1‑150 (MQ=255)
                                                                                                                      cgcAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATcc                           >  2:250201/1‑150 (MQ=255)
                                                                                                                      cgcAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATcc                           >  2:250640/1‑150 (MQ=255)
                                                                                                                         aGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGtt                        >  1:1560435/1‑150 (MQ=255)
                                                                                                                          gATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGttt                       <  1:1978475/150‑1 (MQ=255)
                                                                                                                           aTAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGtttt                      <  1:2417783/150‑1 (MQ=255)
                                                                                                                               aaTATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGAGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGAAGATTTTTGTTGGTCATCTGACCATCCGTTTTGATg                  >  1:3639849/1‑150 (MQ=255)
                                                                                                                                 tatCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGc                <  1:3349749/150‑1 (MQ=255)
                                                                                                                                 tatCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGc                <  1:143881/150‑1 (MQ=255)
                                                                                                                                    cTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGGTTTGATGGCGtg             >  2:2939369/1‑150 (MQ=255)
                                                                                                                                     tGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGtgc            <  2:3568268/150‑1 (MQ=255)
                                                                                                                                      gTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGtgct           <  1:2179296/150‑1 (MQ=255)
                                                                                                                                         catcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTgcg        <  1:3306861/150‑1 (MQ=255)
                                                                                                                                         catcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTgcg        >  2:3343435/1‑150 (MQ=255)
                                                                                                                                             atcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATGGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGAc    >  2:1293645/1‑150 (MQ=255)
                                                                                                                                              tcaATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACg   <  1:3680142/150‑1 (MQ=255)
                                                                                                                                               caAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGt  >  1:58618/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCGATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGT  >  CP009273/1754804‑1755096

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: ‑