Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,077,911 A→C 93.7% K71N (AAA→AAC efeO → inactive ferrous ion transporter EfeUOB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,077,9110AC93.7% 1212.2 / 32.1 307K71N (AAA→AACefeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base C (126/161);  total (136/171)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGG  >  CP009273/1077764‑1078059
                                                                                                                                                   |                                                                                                                                                    
 aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGc                                                                                                                                                   >  2:3740325/1‑150 (MQ=255)
   ggCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGt                                                                                                                                                 <  2:4020659/150‑1 (MQ=255)
    gCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTg                                                                                                                                                <  1:1340758/150‑1 (MQ=255)
      tAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAt                                                                                                                                              >  1:2128613/1‑150 (MQ=255)
        aCGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtgg                                                                                                                                            <  2:368112/150‑1 (MQ=255)
         cgccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggt                                                                                                                                           <  2:3999032/150‑1 (MQ=255)
         cgccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggt                                                                                                                                           <  1:4022569/150‑1 (MQ=255)
          gccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtg                                                                                                                                          >  2:1594269/1‑150 (MQ=255)
            cgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggt                                                                                                                                        >  2:1510034/1‑150 (MQ=255)
            cgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggt                                                                                                                                        >  1:5108926/1‑150 (MQ=255)
             gcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggtg                                                                                                                                       >  1:919963/1‑150 (MQ=255)
              cTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggtgg                                                                                                                                      >  2:4683869/1‑150 (MQ=255)
                 aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAag                                                                                                                                   <  2:1142798/150‑1 (MQ=255)
                   gtgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAagag                                                                                                                                 >  1:784843/1‑150 (MQ=255)
                      ccGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCgg                                                                                                                              <  1:4401305/150‑1 (MQ=255)
                          aGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaaa                                                                                                                          >  1:393874/1‑150 (MQ=255)
                           ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaaaa                                                                                                                         >  1:3502392/1‑150 (MQ=255)
                             tctaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACTCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGTGCGGGAAAata                                                                                                                       <  2:2524933/147‑1 (MQ=255)
                               aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATc                                                                                                                     <  1:3581382/150‑1 (MQ=255)
                                   tGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGccc                                                                                                                 >  2:662616/1‑150 (MQ=255)
                                     aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGACCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                               <  2:1369386/150‑1 (MQ=255)
                                     aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                               <  2:319918/150‑1 (MQ=255)
                                      ccGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTg                                                                                                              <  2:4490772/150‑1 (MQ=255)
                                        gTGACGGATCAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGc                                                                                                            <  2:2779904/150‑1 (MQ=255)
                                            cGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTa                                                                                                        <  2:35616/150‑1 (MQ=255)
                                             ggATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAg                                                                                                       <  2:1855594/150‑1 (MQ=255)
                                               aTAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGcc                                                                                                     <  2:3204058/150‑1 (MQ=255)
                                                tAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCa                                                                                                    >  1:4667242/1‑150 (MQ=255)
                                                  aGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGa                                                                                                  <  1:2877208/150‑1 (MQ=255)
                                                   gCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaa                                                                                                 >  1:3051644/1‑150 (MQ=255)
                                                    cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                                <  1:135561/150‑1 (MQ=255)
                                                    cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                                >  2:4630630/1‑150 (MQ=255)
                                                      gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                                                                              >  1:1629360/1‑150 (MQ=255)
                                                       tGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATg                                                                                             >  1:2024004/1‑150 (MQ=255)
                                                       tGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATg                                                                                             <  2:4749011/150‑1 (MQ=255)
                                                        gCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGa                                                                                            >  1:1285044/1‑150 (MQ=255)
                                                         cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                           <  2:5087127/150‑1 (MQ=255)
                                                         cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                           >  1:1035984/1‑150 (MQ=255)
                                                          gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                          <  2:3851013/150‑1 (MQ=255)
                                                          gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                          >  1:1659047/1‑150 (MQ=255)
                                                          gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                          <  1:2355257/150‑1 (MQ=255)
                                                          gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                          <  1:1381093/150‑1 (MQ=255)
                                                             ccGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCg                                                                                       >  2:4323595/1‑150 (MQ=255)
                                                             ccGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCg                                                                                       >  2:4323555/1‑150 (MQ=255)
                                                              cGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGa                                                                                      <  2:2611452/150‑1 (MQ=255)
                                                              cGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGa                                                                                      >  1:5059023/1‑150 (MQ=255)
                                                                  gACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAttt                                                                                  <  2:4311432/150‑1 (MQ=255)
                                                                    ccATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTAc                                                                                <  2:2112421/150‑1 (MQ=255)
                                                                     cATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACa                                                                               >  2:1488579/1‑150 (MQ=255)
                                                                      aTTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                              >  2:3407380/1‑150 (MQ=255)
                                                                      aTTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                              <  2:4333854/150‑1 (MQ=255)
                                                                       ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                                                                             <  2:1077251/150‑1 (MQ=255)
                                                                       ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                                                                             >  2:1452711/1‑150 (MQ=255)
                                                                         aCGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCt                                                                           <  2:2070535/150‑1 (MQ=255)
                                                                            gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                                                                        >  1:1530056/1‑150 (MQ=255)
                                                                             ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCg                                                                       <  2:1730933/150‑1 (MQ=255)
                                                                               aaCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGaa                                                                     >  2:4861411/1‑150 (MQ=255)
                                                                                  gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                  >  1:849950/1‑150 (MQ=255)
                                                                                  gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                  <  2:2992399/150‑1 (MQ=255)
                                                                                    cGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                <  2:3768548/150‑1 (MQ=255)
                                                                                    cGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                                <  1:660103/150‑1 (MQ=255)
                                                                                     gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat                                                               >  1:382444/1‑150 (MQ=255)
                                                                                       gAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGatat                                                             <  1:3135599/150‑1 (MQ=255)
                                                                                         aaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGa                                                           <  1:3382897/150‑1 (MQ=255)
                                                                                            cacaGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCt                                                        >  2:1777353/1‑150 (MQ=255)
                                                                                             acaGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTg                                                       <  2:423934/150‑1 (MQ=255)
                                                                                                   cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                 <  2:1939069/150‑1 (MQ=255)
                                                                                                   cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                 <  1:2508877/150‑1 (MQ=255)
                                                                                                   cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                 >  1:36887/1‑150 (MQ=255)
                                                                                                   cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                 <  2:4918091/150‑1 (MQ=255)
                                                                                                    attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctg                                                <  1:2924201/150‑1 (MQ=255)
                                                                                                      tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                              >  1:3786029/1‑150 (MQ=255)
                                                                                                      tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                              <  2:501394/150‑1 (MQ=255)
                                                                                                       attCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgctg                                             <  2:2437959/150‑1 (MQ=255)
                                                                                                          cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                          <  1:4238694/150‑1 (MQ=255)
                                                                                                           aGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTa                                         <  2:3089995/150‑1 (MQ=255)
                                                                                                            gAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTaa                                        >  1:2482378/1‑150 (MQ=255)
                                                                                                            gAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTaa                                        >  1:1305361/1‑150 (MQ=255)
                                                                                                             aaCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGACTTTACAGCCTGGCGACTACGATATGACCTGCGGTCTGCTGACTCAc                                       >  1:2734646/1‑150 (MQ=255)
                                                                                                                cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                    <  2:849950/150‑1 (MQ=255)
                                                                                                                cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                    <  2:3327261/150‑1 (MQ=255)
                                                                                                                cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                    <  1:3773731/150‑1 (MQ=255)
                                                                                                                   aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 <  2:2701700/150‑1 (MQ=255)
                                                                                                                   aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 <  1:2394634/150‑1 (MQ=255)
                                                                                                                   aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 <  1:2912643/150‑1 (MQ=255)
                                                                                                                         aaGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAg                           >  1:4964008/1‑150 (MQ=255)
                                                                                                                          aGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGt                          <  1:3740325/150‑1 (MQ=255)
                                                                                                                               cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                     >  2:3642560/1‑150 (MQ=255)
                                                                                                                               cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                     <  1:4904659/150‑1 (MQ=255)
                                                                                                                                tGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                    >  2:1569222/1‑150 (MQ=255)
                                                                                                                                     tGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               >  1:4806930/1‑150 (MQ=255)
                                                                                                                                        gagaTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGt            >  2:582390/1‑150 (MQ=255)
                                                                                                                                          gaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa          <  2:1933587/150‑1 (MQ=255)
                                                                                                                                          gaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa          <  2:1117028/150‑1 (MQ=255)
                                                                                                                                           aTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAg         <  1:2756889/150‑1 (MQ=255)
                                                                                                                                              cTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa      >  1:1105664/1‑150 (MQ=255)
                                                                                                                                               tCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCaa     >  1:398372/1‑150 (MQ=255)
                                                                                                                                                cAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAAc    >  1:2414801/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAAcgg  >  2:3309565/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTGGATCGTCAAAGGTGAGGCAAcgg  >  2:732467/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGG  >  CP009273/1077764‑1078059

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: