Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 705,997 C→T 8.6% A36T (GCG→ACG)  fur ← ferric iron uptake regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273705,9970CT 8.6% 1017.9 / 62.9 291A36T (GCG→ACG) furferric iron uptake regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (126/140);  new base T (15/10);  total (141/150)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

TGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAA  >  CP009273/705849‑706146
                                                                                                                                                    |                                                                                                                                                     
tGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATCTACCGTAGCCAGACCAATTTCTTCACCCAGATCGATCAGACGGTTGTATAAATCTTCCGCa                                                                                                                                                      >  2:3522461/1‑150 (MQ=255)
tGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCa                                                                                                                                                      <  2:1328760/150‑1 (MQ=255)
 gCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCAc                                                                                                                                                     >  1:4366730/1‑150 (MQ=255)
 gCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCAc                                                                                                                                                     <  1:4660013/150‑1 (MQ=255)
  cTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACt                                                                                                                                                    >  2:226633/1‑150 (MQ=255)
  cTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACt                                                                                                                                                    <  1:314935/150‑1 (MQ=255)
  cTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACt                                                                                                                                                    >  2:2902834/1‑150 (MQ=255)
    gtgtCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGa                                                                                                                                                  <  1:5035302/150‑1 (MQ=255)
     tgtCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGAc                                                                                                                                                 <  1:3015272/150‑1 (MQ=255)
     tgtCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGAc                                                                                                                                                 <  1:4872814/150‑1 (MQ=255)
            tCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATgg                                                                                                                                          >  2:3755434/1‑150 (MQ=255)
            tCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATgg                                                                                                                                          >  2:4220822/1‑150 (MQ=255)
             cAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGt                                                                                                                                         >  1:215689/1‑150 (MQ=255)
              aaaTACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGtt                                                                                                                                        >  2:1505067/1‑150 (MQ=255)
              aaaTACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGtt                                                                                                                                        <  1:932362/150‑1 (MQ=255)
               aaTACGGATTTACCGCCTTCATAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTg                                                                                                                                       >  2:4027493/1‑150 (MQ=255)
                aTACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGt                                                                                                                                      <  1:4755130/150‑1 (MQ=255)
                    ggATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCgg                                                                                                                                  >  2:3052450/1‑150 (MQ=255)
                    ggATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGGGATGGTTGTCCgg                                                                                                                                  >  2:2213408/1‑150 (MQ=255)
                        ttACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTcc                                                                                                                              <  2:4410634/150‑1 (MQ=255)
                         tACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCt                                                                                                                             <  1:63341/150‑1 (MQ=255)
                         tACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCt                                                                                                                             <  2:3168119/150‑1 (MQ=255)
                             gCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                                                                                                         <  1:3545053/150‑1 (MQ=255)
                             gCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                                                                                                         <  1:1050024/150‑1 (MQ=255)
                             gCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                                                                                                         <  1:277246/150‑1 (MQ=255)
                              ccTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaaga                                                                                                                        >  2:4872646/1‑150 (MQ=255)
                               cTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                                                                                                       >  1:1576386/1‑150 (MQ=255)
                               cTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                                                                                                       >  2:556139/1‑150 (MQ=255)
                                ttCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAAc                                                                                                                      <  1:1305569/150‑1 (MQ=255)
                                  cAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGTTGTCCGGCTCCTGAAGAACtt                                                                                                                    >  1:1263757/1‑150 (MQ=255)
                                     aaTTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCa                                                                                                                 <  1:2936771/150‑1 (MQ=255)
                                        tgtgGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGa                                                                                                              >  2:48993/1‑150 (MQ=255)
                                           ggCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAttt                                                                                                           <  1:2229602/150‑1 (MQ=255)
                                              gggTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTa                                                                                                        <  1:3360002/150‑1 (MQ=255)
                                                gTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTaaa                                                                                                      >  1:342038/1‑150 (MQ=255)
                                                  gACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACg                                                                                                    <  1:1317196/150‑1 (MQ=255)
                                                  gACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACg                                                                                                    <  2:3631023/150‑1 (MQ=255)
                                                   aCGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGa                                                                                                   >  2:4637646/1‑150 (MQ=255)
                                                      aTACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCCTATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGa                                                                                                <  1:1973941/150‑1 (MQ=255)
                                                       tACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGaa                                                                                               <  1:388679/150‑1 (MQ=255)
                                                       tACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGaa                                                                                               >  2:4939490/1‑150 (MQ=255)
                                                       tACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGaa                                                                                               >  1:1158528/1‑150 (MQ=255)
                                                        aCCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAg                                                                                              <  2:1533985/150‑1 (MQ=255)
                                                        aCCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAg                                                                                              <  2:2495619/150‑1 (MQ=255)
                                                           aGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGt                                                                                           <  2:4397614/150‑1 (MQ=255)
                                                             cgtcgtCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTa                                                                                         <  1:482874/150‑1 (MQ=255)
                                                              gtcgtcCACCTGGTTCAGGACGCGATATCCCGTAGCCAGACCAATTTCTTCACCCATCTCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAACCGAGGAAGCGTTAc                                                                                        <  2:1145215/150‑1 (MQ=255)
                                                              gtcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTAc                                                                                        >  1:3734686/1‑150 (MQ=255)
                                                                 gtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACttt                                                                                     >  1:181861/1‑150 (MQ=255)
                                                                  tcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTc                                                                                    <  2:639923/150‑1 (MQ=255)
                                                                   cAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCa                                                                                   >  2:5027792/1‑150 (MQ=255)
                                                                   cAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCa                                                                                   <  2:4031997/150‑1 (MQ=255)
                                                                    aaaCTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAg                                                                                  <  2:4974830/150‑1 (MQ=255)
                                                                     aaCTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAgg                                                                                 >  2:2679573/1‑150 (MQ=255)
                                                                           ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGct                                                                           >  1:4949349/1‑150 (MQ=255)
                                                                              aGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGACGCGTTACTTTCAGGCCTGctttc                                                                        >  2:3976006/1‑150 (MQ=255)
                                                                              aGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGctttc                                                                        >  1:256137/1‑150 (MQ=255)
                                                                                tACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGctttctt                                                                      >  2:2436067/1‑150 (MQ=255)
                                                                                 aCGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGctttcttt                                                                     <  2:1961784/150‑1 (MQ=255)
                                                                                  cgcgATATACCGTAGCCAGACCAATTTCTTCACCCCTATCGCTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTa                                                                    <  1:1669769/150‑1 (MQ=255)
                                                                                        ataCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCgg                                                              >  1:1020763/1‑150 (MQ=255)
                                                                                         taCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGt                                                             >  2:1394042/1‑150 (MQ=255)
                                                                                             gTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTg                                                         >  1:3368040/1‑150 (MQ=255)
                                                                                                gCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTa                                                      >  2:762303/1‑150 (MQ=255)
                                                                                                 ccAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTAt                                                     >  2:1073542/1‑150 (MQ=255)
                                                                                                  cAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATc                                                    <  2:1101574/150‑1 (MQ=255)
                                                                                                   aGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATcc                                                   >  2:2346568/1‑149 (MQ=255)
                                                                                                   aGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCa                                                   <  2:2061603/150‑1 (MQ=255)
                                                                                                    gACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATccg                                                  <  1:4438947/150‑3 (MQ=255)
                                                                                                        aaTTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCAt                                              >  1:1452070/1‑150 (MQ=255)
                                                                                                          tttCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGc                                            <  2:1022002/150‑1 (MQ=255)
                                                                                                             cttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGa                                         <  1:5027792/150‑1 (MQ=255)
                                                                                                                cACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATc                                      >  2:4676521/1‑150 (MQ=255)
                                                                                                                     atatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTcc                                 >  2:4505091/1‑150 (MQ=255)
                                                                                                                       atCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTg                               >  1:2290099/1‑150 (MQ=255)
                                                                                                                       atCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTg                               >  1:799332/1‑150 (MQ=255)
                                                                                                                          gATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTa                            >  1:2425415/1‑150 (MQ=255)
                                                                                                                          gATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTa                            <  2:800739/150‑1 (MQ=255)
                                                                                                                           aTCAGACGTTTGTATAAATCTTCCGTACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTAc                           >  1:177367/1‑150 (MQ=255)
                                                                                                                           aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTAc                           >  2:1692145/1‑150 (MQ=255)
                                                                                                                             cAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                         >  2:1846298/1‑150 (MQ=255)
                                                                                                                              aGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTaa                        >  1:1689894/1‑150 (MQ=255)
                                                                                                                               gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAg                       >  2:4453320/1‑150 (MQ=255)
                                                                                                                                 cGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCg                     >  1:5102804/1‑150 (MQ=255)
                                                                                                                                  gTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCgg                    <  2:3383253/150‑1 (MQ=255)
                                                                                                                                     tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                 <  1:1711631/150‑1 (MQ=255)
                                                                                                                                          aaaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTc            <  2:1158528/150‑1 (MQ=255)
                                                                                                                                           aaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCa           <  2:5141164/150‑1 (MQ=255)
                                                                                                                                             tCTTCCGCACTGACGTGATGGTTGTCCGTCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt         <  2:2755681/150‑1 (MQ=255)
                                                                                                                                             tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt         <  1:2822172/150‑1 (MQ=255)
                                                                                                                                             tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt         >  1:762285/1‑150 (MQ=255)
                                                                                                                                              cTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTa        <  2:2294131/150‑1 (MQ=255)
                                                                                                                                               ttCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTag       >  1:1824975/1‑150 (MQ=255)
                                                                                                                                                tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTaga      <  2:3457113/150‑1 (MQ=255)
                                                                                                                                                tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTaga      >  2:3843909/1‑150 (MQ=255)
                                                                                                                                                 ccGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaa     <  1:52322/150‑1 (MQ=255)
                                                                                                                                                   gCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaaga   >  1:4027817/1‑150 (MQ=255)
                                                                                                                                                    cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa  <  2:4999770/150‑1 (MQ=255)
                                                                                                                                                    cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa  >  1:2630815/1‑150 (MQ=255)
                                                                                                                                                    cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa  >  2:694840/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAA  >  CP009273/705849‑706146

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: