Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,757,175 G→T 9.4% I364I (ATC→ATA sufB ← 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,757,1750GT 9.4% 648.6 / 28.7 191I364I (ATC→ATAsufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base G (91/82);  new base T (8/10);  total (99/92)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 193 total aligned reads displayed.

ATTAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTGATTGCGTCC  >  CP009273/1757028‑1757320
                                                                                                                                                   |                                                                                                                                                 
aTTAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCg                                                                                                                                                   >  2:1826029/1‑148 (MQ=255)
 ttAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGa                                                                                                                                                  >  2:311307/1‑148 (MQ=255)
 ttAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGa                                                                                                                                                  >  2:3508440/1‑148 (MQ=255)
  tAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGAt                                                                                                                                                 <  2:3439029/148‑1 (MQ=255)
  tAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGAt                                                                                                                                                 <  2:244638/148‑1 (MQ=255)
     cgcCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGAtgtg                                                                                                                                              <  2:1458819/148‑1 (MQ=255)
       ccAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGa                                                                                                                                            >  1:1682166/1‑148 (MQ=255)
        cAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGat                                                                                                                                           >  2:2608252/1‑148 (MQ=255)
         aaTCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGatc                                                                                                                                          >  2:1927846/1‑148 (MQ=255)
         aaTCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGatc                                                                                                                                          <  1:1889723/148‑1 (MQ=255)
          aTCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGatca                                                                                                                                         <  2:3275882/148‑1 (MQ=255)
           tCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGatcat                                                                                                                                        <  2:476734/148‑1 (MQ=255)
           tCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGatcat                                                                                                                                        <  2:3271862/148‑1 (MQ=255)
            cAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGatcatc                                                                                                                                       >  1:3530138/1‑148 (MQ=255)
              gCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCtt                                                                                                                                     >  1:2661824/1‑148 (MQ=255)
               cATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTg                                                                                                                                    <  2:2866889/148‑1 (MQ=255)
               cATTGAGTCGCACTGAGTGAAATTGCGCGCATAGGTTGCCGACGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTg                                                                                                                                    >  1:108220/1‑148 (MQ=255)
                   gAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGc                                                                                                                                <  2:1990736/148‑1 (MQ=255)
                     gTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGATATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCg                                                                                                                              >  1:2288731/1‑148 (MQ=255)
                         cACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTAt                                                                                                                          >  1:3626136/1‑148 (MQ=255)
                           cTGAGTGAAATAGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCAGCGATAACTGTTCTGACTATTTAAGTCAGATATACCTATAGACATAATCGTAGATTTTGTGTTTTTACAGATGTGGATCATATTTGTGAAGGTATcc                                                                                                                        >  2:421752/1‑148 (MQ=255)
                            tGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCg                                                                                                                       >  1:2816379/1‑148 (MQ=255)
                             gAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGc                                                                                                                      >  1:805071/1‑148 (MQ=255)
                              aGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCt                                                                                                                     <  2:802042/148‑1 (MQ=255)
                              aGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCt                                                                                                                     <  2:2718022/148‑1 (MQ=255)
                              aGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCt                                                                                                                     >  1:2123582/1‑148 (MQ=255)
                              aGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCt                                                                                                                     >  2:1668722/1‑148 (MQ=255)
                              aGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCt                                                                                                                     >  1:2268760/1‑148 (MQ=255)
                                tGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCTTg                                                                                                                   >  2:2059134/1‑148 (MQ=255)
                                 gacaTTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGGGTTTTTACCGATGTGGATCATCGTGGTGCCGGTATCCGCTTGc                                                                                                                  >  1:3421860/4‑148 (MQ=255)
                                 gAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGc                                                                                                                  >  2:1270991/1‑148 (MQ=255)
                                 gAAAGTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGc                                                                                                                  <  2:3165824/148‑1 (MQ=255)
                                       gcgcgcATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCtgatga                                                                                                            <  2:1155440/148‑1 (MQ=255)
                                            cATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCt                                                                                                       >  2:729236/1‑148 (MQ=255)
                                            cATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCt                                                                                                       >  2:633109/1‑148 (MQ=255)
                                               tGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGt                                                                                                    <  2:2055788/148‑1 (MQ=255)
                                                  ttGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAg                                                                                                 >  1:52201/1‑148 (MQ=255)
                                                    gCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAgcg                                                                                               <  2:3449754/148‑1 (MQ=255)
                                                       gTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGTGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCa                                                                                            <  1:1620172/148‑1 (MQ=255)
                                                            cATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAg                                                                                       >  2:295111/1‑148 (MQ=255)
                                                            cATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAg                                                                                       >  2:2223100/1‑148 (MQ=255)
                                                            cATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCAGCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAg                                                                                       >  2:2283765/1‑148 (MQ=255)
                                                             aTGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGt                                                                                      >  1:2388962/1‑148 (MQ=255)
                                                              tGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTa                                                                                     >  1:447771/1‑148 (MQ=255)
                                                              tGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTa                                                                                     >  1:3356893/1‑148 (MQ=255)
                                                              tGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTa                                                                                     >  1:1157477/1‑148 (MQ=255)
                                                                    tCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTc                                                                               <  2:2812570/148‑1 (MQ=255)
                                                                          aGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAAt                                                                         >  2:3120932/1‑148 (MQ=255)
                                                                            ccGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATgg                                                                       <  1:1878505/148‑1 (MQ=255)
                                                                                 aTAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTa                                                                  <  2:1987338/148‑1 (MQ=255)
                                                                                 aTAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATAATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTa                                                                  <  2:2386138/148‑1 (MQ=255)
                                                                                   aaCTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATc                                                                >  1:3264322/1‑148 (MQ=255)
                                                                                    aCTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATcg                                                               >  1:3199884/1‑148 (MQ=255)
                                                                                    aCTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATcg                                                               >  1:2144294/1‑148 (MQ=255)
                                                                                      tGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATcgcc                                                             >  2:885195/1‑148 (MQ=255)
                                                                                       gTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATcgccg                                                            >  1:2822190/1‑148 (MQ=255)
                                                                                        ttCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATTTGGATCATCTTGGTGCCGGTATCAGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATcgccgc                                                           >  1:597745/1‑148 (MQ=255)
                                                                                            gACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaa                                                       >  2:2310306/1‑148 (MQ=255)
                                                                                             aCTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaa                                                      <  2:881221/148‑1 (MQ=255)
                                                                                             aCTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaa                                                      <  2:2652429/148‑1 (MQ=255)
                                                                                             aCTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaa                                                      <  1:1492176/148‑1 (MQ=255)
                                                                                              cTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaaa                                                     <  2:1143930/148‑1 (MQ=255)
                                                                                              cTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaaa                                                     <  1:1499007/148‑1 (MQ=255)
                                                                                              cTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCaaaa                                                     <  1:3398594/148‑1 (MQ=255)
                                                                                                 tGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGc                                                  <  1:830327/148‑1 (MQ=255)
                                                                                                 tGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGc                                                  <  1:3280126/148‑1 (MQ=255)
                                                                                                 tGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGc                                                  <  1:3280116/148‑1 (MQ=255)
                                                                                                  gTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCa                                                 >  1:378386/1‑148 (MQ=255)
                                                                                                   tCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAg                                                <  1:1634741/148‑1 (MQ=255)
                                                                                                   tCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAg                                                <  2:3477856/148‑1 (MQ=255)
                                                                                                    ccGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGc                                               <  2:2231665/148‑1 (MQ=255)
                                                                                                       gCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTgg                                            <  1:19719/148‑1 (MQ=255)
                                                                                                       gCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTgg                                            >  2:2968893/1‑148 (MQ=255)
                                                                                                        cAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTggg                                           <  1:3083789/148‑1 (MQ=255)
                                                                                                            gaTCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGatat                                       >  1:1802057/1‑148 (MQ=255)
                                                                                                                ccTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTcc                                   <  2:1658674/148‑1 (MQ=255)
                                                                                                                   ttCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACg                                >  2:2892990/1‑148 (MQ=255)
                                                                                                                     cGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTa                              <  1:1270991/148‑1 (MQ=255)
                                                                                                                      gagaTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTaa                             <  1:2608252/148‑1 (MQ=255)
                                                                                                                      gagaTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTaa                             <  1:764690/148‑1 (MQ=255)
                                                                                                                        gaTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATc                           >  1:3372769/1‑148 (MQ=255)
                                                                                                                         aTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAAGTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCg                          <  1:1753843/148‑1 (MQ=255)
                                                                                                                         aTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCg                          >  1:1064504/1‑148 (MQ=255)
                                                                                                                         aTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCg                          <  2:38763/148‑1 (MQ=255)
                                                                                                                         aTAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATATTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCAACGTAATCg                          >  2:3198378/1‑148 (MQ=255)
                                                                                                                          tAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGc                         >  2:1006919/1‑148 (MQ=255)
                                                                                                                           aaTGGTCGATTTGGTGTTTTTACCGATGTGGAACATCTTGGAGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGAAATctct                        >  2:3806013/1‑148 (MQ=255)
                                                                                                                               gTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGAcc                    <  2:660054/148‑1 (MQ=255)
                                                                                                                               gTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGAcc                    >  2:3162776/1‑148 (MQ=255)
                                                                                                                                  gATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCgg                 <  2:1700635/148‑1 (MQ=255)
                                                                                                                                   aTTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGt                <  2:3575173/148‑1 (MQ=255)
                                                                                                                                     ttGGTGTTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGttt              <  2:772228/148‑1 (MQ=255)
                                                                                                                                     ttGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGttt              <  2:2760046/148‑1 (MQ=255)
                                                                                                                                     ttGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGttt              <  2:352650/148‑1 (MQ=255)
                                                                                                                                     ttGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGttt              <  2:1565961/148‑1 (MQ=255)
                                                                                                                                        gtgtTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTg           <  2:3696418/148‑1 (MQ=255)
                                                                                                                                          gtTTTTACCTATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTGAt         <  2:330631/148‑1 (MQ=255)
                                                                                                                                          gtTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTGAt         >  1:3564815/1‑148 (MQ=255)
                                                                                                                                           tttttACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTGAtt        >  1:3010530/1‑148 (MQ=255)
                                                                                                                                                aCCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTGCACGTAATCGCTGACCCGGTTTCTGATTGCGTc   >  2:1625894/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATTAGCGCCAATCAGCATTGAGTCGCACTGAGTGAAATTGCGCGCATTGGTTGCCGTCGGCATGATTTTCACTAAGCCGCGATAACTGTTCTGACTATGTCCGGCAGAGATCCCTTTCGAGATAATGGTCGATTTGGTGTTTTTACCGATGTGGATCATCTTGGTGCCGGTATCCGCTTGCTGATGACCGCTGGTCAGCGCCACTGAGTAAAACTCACCAATGGAGTTATCGCCGCGCAAAATGCAGCTGGGATATTTCCACGTAATCGCTGACCCGGTTTCTGATTGCGTCC  >  CP009273/1757028‑1757320

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: