Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,754,950 T→A 89.9% 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,9500TA89.9% 817.2 / 72.0 281I367F (ATT→TTT) sufDcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base A (122/123);  total (138/143)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

Only 100 of 281 total aligned reads displayed.

GACCGATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGTTG  >  CP009273/1754803‑1755098
                                                                                                                                                   |                                                                                                                                                    
 aCCGATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACg                                                                                                                                                   <  2:2813142/150‑1 (MQ=255)
     aTTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCTGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGTTAGAATATCTGTTCATCATCAAAACGcccc                                                                                                                                               <  2:4420442/150‑1 (MQ=255)
      ttCGGGCCAGCACCTGCTTTTTAAGCCCCTCATCACGCTGTGCTTCCGTCTGTTCGGCAGCGAAGGCGTAAATGATAATCTGCTGGGCATCCTGCTTTTTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCa                                                                                                                                              <  1:551550/150‑1 (MQ=255)
      ttCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCa                                                                                                                                              <  1:119802/150‑1 (MQ=255)
         gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCg                                                                                                                                           <  2:2364914/150‑1 (MQ=255)
         gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCg                                                                                                                                           >  2:35787/1‑150 (MQ=255)
         gggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCg                                                                                                                                           >  2:4277077/1‑150 (MQ=255)
          ggCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGt                                                                                                                                          >  1:3098708/1‑150 (MQ=255)
            ccAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTcg                                                                                                                                        >  1:1567071/1‑150 (MQ=255)
             cAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTcgc                                                                                                                                       >  2:1186107/1‑150 (MQ=255)
             cAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTcgc                                                                                                                                       >  2:2409925/1‑150 (MQ=255)
                cACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGcc                                                                                                                                    >  2:3690352/1‑150 (MQ=255)
                   ctgctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTg                                                                                                                                 <  1:580454/150‑1 (MQ=255)
                   ctgctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTg                                                                                                                                 <  2:199572/150‑1 (MQ=255)
                   ctgctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTg                                                                                                                                 >  2:3900319/1‑150 (MQ=255)
                     gctgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGc                                                                                                                               >  1:451109/1‑150 (MQ=255)
                       tgTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTg                                                                                                                             <  2:4663443/150‑1 (MQ=255)
                          ttAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAt                                                                                                                          <  1:4396606/150‑1 (MQ=255)
                            aaGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCAttt                                                                                                                        <  1:1378711/150‑1 (MQ=255)
                             aGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCGGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTc                                                                                                                       <  2:590449/150‑1 (MQ=255)
                              gCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTca                                                                                                                      <  2:1590915/150‑1 (MQ=255)
                               ccccTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTcac                                                                                                                     >  2:2920016/1‑150 (MQ=255)
                                 ccTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCAcat                                                                                                                   <  2:1523818/150‑1 (MQ=255)
                                 ccTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCAcat                                                                                                                   <  2:722676/150‑1 (MQ=255)
                                    catcaCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCAcatcat                                                                                                                >  1:3586448/1‑150 (MQ=255)
                                     atcaCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACatcatc                                                                                                               >  2:2634619/1‑150 (MQ=255)
                                      tcaCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCAACCGTCGCGCCGTGGCTGCATTTCACATCATCt                                                                                                              >  1:386229/1‑150 (MQ=255)
                                       caCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTg                                                                                                             >  1:4605617/1‑150 (MQ=255)
                                        aCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGc                                                                                                            <  1:4105093/150‑1 (MQ=255)
                                        aCGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGc                                                                                                            <  1:637017/150‑1 (MQ=255)
                                            aGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTACTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAg                                                                                                        <  1:4296098/150‑1 (MQ=255)
                                             gTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGa                                                                                                       >  2:3104418/1‑150 (MQ=255)
                                              tgattCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATGTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCTTATAGATTATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGAt                                                                                                      <  1:3668661/147‑1 (MQ=255)
                                              tGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGAt                                                                                                      <  1:246112/150‑1 (MQ=255)
                                               gCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGAtt                                                                                                     <  2:3022770/150‑1 (MQ=255)
                                               gCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGAtt                                                                                                     <  2:4043392/150‑1 (MQ=255)
                                                  tCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTcc                                                                                                  >  1:204319/1‑150 (MQ=255)
                                                   ccGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCa                                                                                                 >  1:1487985/1‑150 (MQ=255)
                                                   ccGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAACACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCa                                                                                                 <  1:883939/150‑1 (MQ=255)
                                                   ccGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCa                                                                                                 <  2:4768784/150‑1 (MQ=255)
                                                      tCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCt                                                                                              <  2:2261742/150‑1 (MQ=255)
                                                       cAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCCGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTg                                                                                             >  2:1545352/1‑150 (MQ=255)
                                                        aGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGc                                                                                            >  2:54198/1‑150 (MQ=255)
                                                         gTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCg                                                                                           >  2:922432/1‑150 (MQ=255)
                                                             ggCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGttt                                                                                       >  2:1739537/1‑150 (MQ=255)
                                                               cAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCg                                                                                     <  1:695281/150‑1 (MQ=255)
                                                                  cGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTAt                                                                                  <  2:1862410/150‑1 (MQ=255)
                                                                    aaGGCGTAAATGATCATCTGCTGGTCATCCTGCTGATTGTTCCCGCTCTTGCGCAGATTGAATATCTGTTCATCATCAAACCGCCCCACCGTCGCGCCGTGGCTGCATTTCGCATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATcc                                                                                <  1:1079890/150‑1 (MQ=255)
                                                                    aaGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATcc                                                                                >  1:2534235/1‑150 (MQ=255)
                                                                       gCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACt                                                                             <  2:2940349/150‑1 (MQ=255)
                                                                       gCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACt                                                                             <  1:342849/150‑1 (MQ=255)
                                                                            aaTGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGc                                                                        >  2:1409394/1‑150 (MQ=255)
                                                                             aTGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGcc                                                                       >  1:1769646/1‑150 (MQ=255)
                                                                               gATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAg                                                                     >  1:3025108/1‑150 (MQ=255)
                                                                                  catcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGttt                                                                  <  1:4121600/150‑1 (MQ=255)
                                                                                  catcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGttt                                                                  <  2:2980089/150‑1 (MQ=255)
                                                                                  catcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGttt                                                                  <  2:214223/150‑1 (MQ=255)
                                                                                   atcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTg                                                                 <  1:2582301/150‑1 (MQ=255)
                                                                                   atcTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTg                                                                 <  2:3398403/150‑1 (MQ=255)
                                                                                        ctgGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCAt                                                            >  2:3776/1‑150 (MQ=255)
                                                                                                ccTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAt                                                    <  2:104563/150‑1 (MQ=255)
                                                                                                 ctgctgATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAtt                                                   <  2:2774316/150‑1 (MQ=255)
                                                                                                   gctgATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttgt                                                 >  1:795744/1‑150 (MQ=255)
                                                                                                      gattgatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGAttgttgt                                              <  2:1784071/150‑1 (MQ=255)
                                                                                                       attgatCCCGCGCGAGCGCATATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCAACTGCCTAGATTTCCAGCTGCGGTTTCGTATCCACTACCGACAGTTTGCCCAGCAGCAGAgtgttgtt                                             >  1:2175643/1‑150 (MQ=255)
                                                                                                          gatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGt                                          <  2:4793489/150‑1 (MQ=255)
                                                                                                          gatCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGt                                          >  1:2731880/1‑150 (MQ=255)
                                                                                                              ccGCGCGAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATc                                      >  2:3639350/1‑150 (MQ=255)
                                                                                                                 cgcgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGa                                   <  2:4393148/150‑1 (MQ=255)
                                                                                                                  gcgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGAc                                  >  2:4427473/1‑150 (MQ=255)
                                                                                                                   cgAGCGCAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGAcc                                 >  1:2117258/1‑150 (MQ=255)
                                                                                                                         cAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCTTGGCTGCATTTCAAATCATATGCATAGATTTCCAGCTGAGGTTTCGTATCAAATTAAGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGAACATCAGt                           >  1:4376773/1‑150 (MQ=255)
                                                                                                                         cAGATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGt                           >  1:2378969/1‑150 (MQ=255)
                                                                                                                           gATAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGttt                         <  1:2980394/150‑1 (MQ=255)
                                                                                                                            aTAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGtttt                        <  1:4712955/150‑1 (MQ=255)
                                                                                                                            aTAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGtttt                        <  1:445169/150‑1 (MQ=255)
                                                                                                                            aTAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGtttt                        <  2:271461/150‑1 (MQ=255)
                                                                                                                             tAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTg                       <  1:853189/150‑1 (MQ=255)
                                                                                                                             tAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTg                       >  2:3499397/1‑150 (MQ=255)
                                                                                                                             tAGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTg                       >  1:1035899/1‑150 (MQ=255)
                                                                                                                              aGAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGa                      >  1:740712/1‑150 (MQ=255)
                                                                                                                               gAATATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGAt                     >  2:532566/1‑150 (MQ=255)
                                                                                                                                aaTATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATg                    <  1:3900319/150‑1 (MQ=255)
                                                                                                                                aaTATCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATg                    >  1:3647081/1‑150 (MQ=255)
                                                                                                                                  tatCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGc                  >  1:4012967/1‑150 (MQ=255)
                                                                                                                                   atCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTACCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCg                 >  2:409316/1‑150 (MQ=255)
                                                                                                                                    tCTGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGt                <  2:3360673/150‑1 (MQ=255)
                                                                                                                                      tGTTCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGtgc              >  2:1269367/1‑150 (MQ=255)
                                                                                                                                        ttCATCATCAAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGtgctg            <  1:3308456/150‑1 (MQ=255)
                                                                                                                                         tcatcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGtgctgc           >  1:2664350/1‑150 (MQ=255)
                                                                                                                                           atcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTgcgc         >  2:2932185/1‑150 (MQ=255)
                                                                                                                                           atcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTgcgc         <  1:1943070/150‑1 (MQ=255)
                                                                                                                                            tcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGcgcg        <  1:35729/150‑1 (MQ=255)
                                                                                                                                            tcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGcgcg        >  1:4770242/1‑150 (MQ=255)
                                                                                                                                            tcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGcgcg        >  2:1670308/1‑150 (MQ=255)
                                                                                                                                            tcatcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGcgcg        >  1:1121268/1‑150 (MQ=255)
                                                                                                                                             catcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGa       >  2:1924926/1‑150 (MQ=255)
                                                                                                                                              atcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGAc      >  2:2667927/1‑150 (MQ=255)
                                                                                                                                               tcaAAACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACg     <  1:4644369/150‑1 (MQ=255)
                                                                                                                                                  aTACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGTTg  >  1:3510284/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCGATTCGGGCCAGCACCTGCTGTTTAAGCCCCTCATCACGCAGTGCTTCCGTCAGTTCGGCAGCGAAGGCGTAAATGATCATCTGCTGGGCATCCTGCTGATTGATCCCGCGCGAGCGCAGATAGAATATCTGTTCATCATCAATACGCCCCACCGTCGCGCCGTGGCTGCATTTCACATCATCTGCATAGATTTCCAGCTGCGGTTTCGTATCCACTTCCGCCAGTTTGCCCATCAGCAGATTGTTGTTGGTCATCTGACCATCCGTTTTGATGGCGTGCTGCGCGACGTTG  >  CP009273/1754803‑1755098

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: