Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,077,911 A→C 71.8% 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,9110AC71.8% 727.6 / 254.2 339K71N (AAA→AACefeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (41/56);  new base C (106/136);  total (147/192)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

TTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGG  >  CP009273/1077762‑1078059
                                                                                                                                                     |                                                                                                                                                    
   aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGc                                                                                                                                                   <  1:1481212/150‑1 (MQ=255)
   aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGc                                                                                                                                                   >  1:2717165/1‑150 (MQ=255)
    tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCg                                                                                                                                                  >  1:765316/1‑150 (MQ=255)
           cgccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggt                                                                                                                                           <  2:2777362/150‑1 (MQ=255)
             ccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtgg                                                                                                                                         <  1:2923327/150‑1 (MQ=255)
                cTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggtgg                                                                                                                                      >  1:410011/1‑150 (MQ=255)
                cTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggtgg                                                                                                                                      >  2:630811/1‑150 (MQ=255)
                  gATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGaa                                                                                                                                    <  1:1460478/150‑1 (MQ=255)
                  gATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGaa                                                                                                                                    >  1:4201695/1‑150 (MQ=255)
                   aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAag                                                                                                                                   <  1:1262792/150‑1 (MQ=255)
                   aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAAGGTGGTGGAag                                                                                                                                   >  2:4372498/1‑150 (MQ=255)
                   aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAag                                                                                                                                   >  2:4820864/1‑150 (MQ=255)
                    tgtgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAaga                                                                                                                                  <  1:3492507/150‑1 (MQ=255)
                     gtgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAagag                                                                                                                                 <  2:3594508/150‑1 (MQ=255)
                         cGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCggg                                                                                                                             <  1:1057004/150‑1 (MQ=255)
                          gCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGa                                                                                                                            <  1:3607100/150‑1 (MQ=255)
                           cAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                           <  2:2957414/150‑1 (MQ=255)
                           cAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                           >  2:1463059/1‑150 (MQ=255)
                           cAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                           >  2:3883528/1‑150 (MQ=255)
                            aGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaaa                                                                                                                          >  2:2495850/1‑150 (MQ=255)
                                cAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAatat                                                                                                                      <  1:2273520/150‑1 (MQ=255)
                                cAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAatat                                                                                                                      <  1:847114/150‑1 (MQ=255)
                                  aaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCg                                                                                                                    >  2:1004057/1‑150 (MQ=255)
                                   aGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGc                                                                                                                   >  2:3454598/1‑150 (MQ=255)
                                     tGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGccc                                                                                                                 <  2:3800461/150‑1 (MQ=255)
                                      gACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                                                                <  2:4503665/150‑1 (MQ=255)
                                       aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                               <  2:2536913/150‑1 (MQ=255)
                                          gTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGc                                                                                                            >  1:3750763/1‑150 (MQ=255)
                                           tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCt                                                                                                           >  1:1379154/1‑150 (MQ=255)
                                            gACGGATAAGCAGTGCGAACCGATGACCATTACGGTTACCGCCGGGAAAACACAGTACATTATTCAGAACCACCGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCtt                                                                                                          <  1:2343500/150‑1 (MQ=255)
                                             aCGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                                                         <  1:2579708/150‑1 (MQ=255)
                                             aCGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                                                         <  1:4408756/150‑1 (MQ=255)
                                                gATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                      >  1:4430659/1‑150 (MQ=255)
                                                   aaGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAg                                                                                                   >  2:739386/1‑150 (MQ=255)
                                                      cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                                >  1:181986/1‑150 (MQ=255)
                                                       aGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaaa                                                                                               >  2:841088/1‑150 (MQ=255)
                                                         tGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATg                                                                                             <  2:2010664/150‑1 (MQ=255)
                                                            gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                          >  1:1023198/1‑150 (MQ=255)
                                                             aaCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACgg                                                                                         >  2:4615536/1‑150 (MQ=255)
                                                              aCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGc                                                                                        <  2:313287/150‑1 (MQ=255)
                                                              aCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGc                                                                                        <  1:4758156/150‑1 (MQ=255)
                                                                 gATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGaa                                                                                     >  1:1729647/1‑150 (MQ=255)
                                                                  aTGACCTTTACGGTTAACGCCGGGAAAACACTGTTCATTATTCAGAACCACAGCCTGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAt                                                                                    <  1:4598023/150‑1 (MQ=255)
                                                                       cATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACa                                                                               <  1:724301/150‑1 (MQ=255)
                                                                         ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                                                                             >  2:1949999/1‑150 (MQ=255)
                                                                           aCGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCt                                                                           <  1:4290141/150‑1 (MQ=255)
                                                                              gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                                                                        <  2:1259512/150‑1 (MQ=255)
                                                                              gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                                                                        <  2:1910566/150‑1 (MQ=255)
                                                                                tAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGa                                                                      >  1:4503985/1‑150 (MQ=255)
                                                                                   cGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                  <  1:1298170/150‑1 (MQ=255)
                                                                                    gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATAc                                                                  <  2:3506963/150‑1 (MQ=255)
                                                                                     ccGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACg                                                                 >  1:1619900/1‑150 (MQ=255)
                                                                                          aaaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATg                                                            <  2:3208512/150‑1 (MQ=255)
                                                                                          aaaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAGGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATg                                                            >  2:3685207/1‑150 (MQ=255)
                                                                                           aaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGa                                                           >  2:1438816/1‑150 (MQ=255)
                                                                                            aaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGAc                                                          <  2:1337748/150‑1 (MQ=255)
                                                                                              cacaGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCt                                                        <  2:2222905/150‑1 (MQ=255)
                                                                                                 aGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCg                                                     >  2:4793844/1‑150 (MQ=255)
                                                                                                  gTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCgg                                                    >  1:2431158/1‑150 (MQ=255)
                                                                                                     cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                                 >  1:2810162/1‑150 (MQ=255)
                                                                                                         attCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgctg                                             <  2:61867/150‑1 (MQ=255)
                                                                                                           tCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAc                                           >  1:1957080/1‑150 (MQ=255)
                                                                                                            cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                          <  2:3234172/150‑1 (MQ=255)
                                                                                                            cAGAACCACAGCCAGAAGGCGCTGAGTGGGAGATCCTCAACGGCGTGATGGGGGTGGAAGATCGGGAAAATATCGCCCCTGGCTTTAGCCCGAAAATGACGGCGAATTTACAGCCTGGCGCATACGATATGACCTGCGGTCTGCTGACTa                                         >  2:2309459/1‑150 (MQ=255)
                                                                                                             aGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGACTTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTa                                         >  1:3011137/1‑150 (MQ=255)
                                                                                                               aaCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAc                                       <  2:2709757/150‑1 (MQ=255)
                                                                                                                aCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAcc                                      <  2:3844698/150‑1 (MQ=255)
                                                                                                                  cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                    <  1:1913218/150‑1 (MQ=255)
                                                                                                                  cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                    <  2:2848416/150‑1 (MQ=255)
                                                                                                                     aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 >  2:3403330/1‑150 (MQ=255)
                                                                                                                     aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 >  2:2073091/1‑150 (MQ=255)
                                                                                                                     aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                                 >  2:3665043/1‑150 (MQ=255)
                                                                                                                         aGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGa                             <  1:2973658/150‑1 (MQ=255)
                                                                                                                         aGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGa                             >  2:3960515/1‑150 (MQ=255)
                                                                                                                           aaGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAg                           <  2:1719335/150‑1 (MQ=255)
                                                                                                                           aaGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAg                           <  1:4785240/150‑1 (MQ=255)
                                                                                                                              gcgcTGGTGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTg                        <  1:202308/150‑1 (MQ=255)
                                                                                                                              gcgcTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTg                        <  2:4231105/150‑1 (MQ=255)
                                                                                                                                gcTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGAt                      >  2:393365/1‑150 (MQ=255)
                                                                                                                                gcTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGAt                      <  1:1577585/150‑1 (MQ=255)
                                                                                                                                gcTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGAt                      <  1:724913/150‑1 (MQ=255)
                                                                                                                                gcTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGAt                      <  1:3231462/150‑1 (MQ=255)
                                                                                                                                 cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTGCAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                     <  1:1701143/150‑1 (MQ=255)
                                                                                                                                 cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATc                     >  1:1185040/1‑150 (MQ=255)
                                                                                                                                  tGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                    <  2:1217952/150‑1 (MQ=255)
                                                                                                                                   ggAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGt                   >  1:3818914/1‑150 (MQ=255)
                                                                                                                                    gAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTc                  <  1:1079223/150‑1 (MQ=255)
                                                                                                                                       tGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               <  1:3950813/150‑1 (MQ=255)
                                                                                                                                       tGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               >  2:1183566/1‑150 (MQ=255)
                                                                                                                                       tGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               >  1:3151427/1‑150 (MQ=255)
                                                                                                                                       tGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               <  1:3859365/150‑1 (MQ=255)
                                                                                                                                       tGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa               >  2:1608994/1‑150 (MQ=255)
                                                                                                                                        gggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAg              <  2:3010877/150‑1 (MQ=255)
                                                                                                                                             aTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAg         >  2:1844967/1‑150 (MQ=255)
                                                                                                                                             aTCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAg         <  2:4201695/150‑1 (MQ=255)
                                                                                                                                                cTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa      <  2:263269/150‑1 (MQ=255)
                                                                                                                                                cTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa      <  1:3922262/150‑1 (MQ=255)
                                                                                                                                                cTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa      <  2:181986/150‑1 (MQ=255)
                                                                                                                                                 tCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCaa     <  2:3305191/150‑1 (MQ=255)
                                                                                                                                                 tCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCaa     <  1:4820864/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAACGG  >  CP009273/1077762‑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: