Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 619,410 T→C 7.0% T264A (ACG→GCG)  fepC ← ferric enterobactin ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,4100TC 7.0% 254.0 / 7.6 86T264A (ACG→GCG) fepCferric enterobactin ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (43/37);  new base C (3/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

GGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCA  >  NC_000913/619262‑619559
                                                                                                                                                    |                                                                                                                                                     
ggTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGtt                                                                                                                                                      <  1:1085159/150‑1 (MQ=255)
  tAAGGCGATTATTGTGATCCTATCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCg                                                                                                                                                    <  2:516057/150‑1 (MQ=255)
    aGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTc                                                                                                                                                  >  2:1482599/1‑150 (MQ=255)
       cGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCGGa                                                                                                                                               >  1:1117875/1‑150 (MQ=255)
       cGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGa                                                                                                                                               >  2:531650/1‑150 (MQ=255)
          ttattGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCACGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCg                                                                                                                                            <  2:297849/150‑1 (MQ=255)
              tgtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcac                                                                                                                                        <  1:1283001/150‑1 (MQ=255)
              tgtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcac                                                                                                                                        <  1:529591/150‑1 (MQ=255)
              tgtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcac                                                                                                                                        <  2:1352373/150‑1 (MQ=255)
               gtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacc                                                                                                                                       >  1:899588/1‑150 (MQ=255)
                 gATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcaccac                                                                                                                                     >  1:1256851/1‑150 (MQ=255)
                  aTCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACcaca                                                                                                                                    >  1:130064/1‑150 (MQ=255)
                     cTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgc                                                                                                                                 <  1:1293266/150‑1 (MQ=255)
                     cTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcaccttcaac                                                                                                                                 <  1:713181/150‑7 (MQ=255)
                        tCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcggc                                                                                                                              >  1:1149184/1‑150 (MQ=255)
                         ccGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcggcg                                                                                                                             >  2:927703/1‑150 (MQ=255)
                             gTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCGTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTcc                                                                                                                         <  1:1224534/150‑1 (MQ=255)
                                gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGc                                                                                                                      >  2:561325/1‑150 (MQ=255)
                                gATCGGCGGGATGGGGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGc                                                                                                                      <  2:839126/150‑1 (MQ=255)
                                 aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTATTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCa                                                                                                                     <  2:859664/150‑1 (MQ=255)
                                 aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCa                                                                                                                     >  1:636959/1‑150 (MQ=255)
                                 aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCa                                                                                                                     >  1:599007/1‑150 (MQ=255)
                                     gcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAATTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTg                                                                                                                 <  1:969004/150‑1 (MQ=255)
                                      cggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTgg                                                                                                                <  1:1035483/150‑1 (MQ=255)
                                         gATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATc                                                                                                             <  1:40797/150‑1 (MQ=255)
                                          attttggCTTCTTGTAGCGTACTATTGATATATCCTATGGCATACAGAAAACAGGATGACATTATAACAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCg                                                                                                            >  2:408371/5‑150 (MQ=255)
                                           tggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGt                                                                                                           >  2:703054/1‑150 (MQ=255)
                                            ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTc                                                                                                          <  2:295748/150‑1 (MQ=255)
                                            ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGAGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTc                                                                                                          >  1:569946/1‑150 (MQ=255)
                                             gtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCa                                                                                                         <  2:1149184/150‑1 (MQ=255)
                                             gtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCa                                                                                                         <  2:358274/150‑1 (MQ=255)
                                                gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATg                                                                                                      >  1:677834/1‑150 (MQ=255)
                                                gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATg                                                                                                      <  1:871435/150‑1 (MQ=255)
                                                 cTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGa                                                                                                     >  2:423209/1‑150 (MQ=255)
                                                       gtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGc                                                                                               >  2:977866/1‑150 (MQ=255)
                                                       gtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGc                                                                                               <  1:1210338/150‑1 (MQ=255)
                                                        tggCGTGCTGTTTCGCTATGCGAGGGCATCCAGAAAACAGGATGCGATGATGGCAGACCACTTAGTTTAAATATTTTCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCa                                                                                              <  2:239390/150‑1 (MQ=255)
                                                                 gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGAcc                                                                                     <  2:425763/150‑1 (MQ=255)
                                                                 gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGAcc                                                                                     <  1:531650/150‑1 (MQ=255)
                                                                   tGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCAt                                                                                   <  2:329794/150‑1 (MQ=255)
                                                                    gcgcTATGCGATGGCATCCAGAAAACATGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATa                                                                                  >  1:945341/1‑150 (MQ=255)
                                                                           gCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcgc                                                                           >  1:97200/1‑150 (MQ=255)
                                                                            tcaTGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCt                                                                          >  2:1411781/3‑150 (MQ=255)
                                                                            cGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCt                                                                          >  1:1056871/1‑150 (MQ=255)
                                                                             gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTc                                                                         <  2:338558/150‑1 (MQ=255)
                                                                             gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTc                                                                         <  2:636959/150‑1 (MQ=255)
                                                                              aTGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCa                                                                        >  2:1083361/1‑150 (MQ=255)
                                                                                  cATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCa                                                                    >  1:614589/1‑150 (MQ=255)
                                                                                   aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAg                                                                   >  2:165586/1‑150 (MQ=255)
                                                                                   aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAg                                                                   <  2:368420/150‑1 (MQ=255)
                                                                                    tCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGt                                                                  >  1:445653/1‑150 (MQ=255)
                                                                                       aGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTca                                                               <  1:309408/150‑1 (MQ=255)
                                                                                        gAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcag                                                              <  1:1482599/150‑1 (MQ=255)
                                                                                             cAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTg                                                         >  2:294828/1‑150 (MQ=255)
                                                                                               ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGAc                                                       <  2:167290/150‑1 (MQ=255)
                                                                                                 aTGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACaa                                                     >  2:1472496/1‑150 (MQ=255)
                                                                                                  tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                    >  1:137992/1‑150 (MQ=255)
                                                                                                  tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                    >  1:298504/1‑150 (MQ=255)
                                                                                                  tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                    >  2:528136/1‑150 (MQ=255)
                                                                                                   gCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAtc                                                   >  1:1216402/1‑150 (MQ=255)
                                                                                                    ccATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAtct                                                  >  1:1184904/1‑150 (MQ=255)
                                                                                                       tgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCt                                               <  1:1411781/150‑1 (MQ=255)
                                                                                                        gatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCtt                                              >  1:230136/1‑150 (MQ=255)
                                                                                                           ggCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCgg                                           <  1:610200/150‑1 (MQ=255)
                                                                                                            gCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgc                                          <  2:445653/150‑1 (MQ=255)
                                                                                                                   aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATTATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCt                                   <  2:1025428/150‑1 (MQ=255)
                                                                                                                   aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgc                                          <  1:422195/143‑1 (MQ=255)
                                                                                                                   aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgc                                          >  2:422195/1‑143 (MQ=255)
                                                                                                                      tAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAg                                <  1:408371/150‑1 (MQ=255)
                                                                                                                       aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGc                               <  2:1256851/150‑1 (MQ=255)
                                                                                                                       aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCAGAGc                               >  2:1483004/1‑150 (MQ=255)
                                                                                                                           tAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaaca                           >  2:1402625/1‑150 (MQ=255)
                                                                                                                            aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaacaa                          >  2:823238/1‑150 (MQ=255)
                                                                                                                              cTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtt                        >  2:430985/1‑150 (MQ=255)
                                                                                                                                aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtttt                      >  2:555556/1‑150 (MQ=255)
                                                                                                                                aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtttt                      >  2:861757/1‑150 (MQ=255)
                                                                                                                                  ttGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                    <  2:899588/150‑1 (MQ=255)
                                                                                                                                      ggTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAACCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttc                >  2:678452/1‑150 (MQ=255)
                                                                                                                                           aaGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgca           >  1:1259442/1‑150 (MQ=255)
                                                                                                                                            tggTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaa          >  1:600526/2‑150 (MQ=255)
                                                                                                                                             ggTGCCGTTCGTCCGAGCGGCCCCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCCATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaat         >  2:477631/1‑150 (MQ=255)
                                                                                                                                             ggTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaat         <  1:823238/150‑1 (MQ=255)
                                                                                                                                              gTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatg        <  1:977866/150‑1 (MQ=255)
                                                                                                                                              gTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatg        >  1:933291/1‑150 (MQ=255)
                                                                                                                                               tGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgc       >  2:1261715/1‑150 (MQ=255)
                                                                                                                                                gCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgca      <  1:1083361/150‑1 (MQ=255)
                                                                                                                                                 ccGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgcaa     <  2:1184904/150‑1 (MQ=255)
                                                                                                                                                    ttCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCa  <  1:555556/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCA  >  NC_000913/619262‑619559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: