Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,666 A→C 5.2% intergenic (‑156/‑219) fepB ← / → entC ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,6660AC 5.2% 294.1 / 5.0 88intergenic (‑156/‑219)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (41/42);  new base C (2/3);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGC  >  NC_000913/624518‑624814
                                                                                                                                                    |                                                                                                                                                    
cTTCCTTTTATTAATTTGTTTAGGGGCTTAACCAAAATTACCAAAAAGCCCAAAACCCGCCCCCCCGCCCTTTGGGGGTGTGTGTTCGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                                                                                     <  1:262825/150‑1 (MQ=255)
cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                                                                                     >  1:1329577/1‑150 (MQ=255)
cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                                                                                     >  1:1329581/1‑150 (MQ=255)
 ttCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAg                                                                                                                                                    >  1:707600/1‑150 (MQ=255)
  tCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                                                                                                   >  1:140404/1‑150 (MQ=255)
      gTTATTACTAAGGTTCGGGGCGAACCCCCACCTCCGACAACGCGCCCAACCCGGCCCCGCGCCCCGTGGGGGTGGGGGTTTGGGTGTGGGTCACCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                                                                                               <  1:484276/150‑1 (MQ=255)
          ttaATCAGGTTCTGGGCTTAATCACAAAGTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGCGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                                                                                                           <  1:1016455/150‑1 (MQ=25)
           taatacGGCTAAGGGCTTAATGAAAAAGTCGCCAAAGCGCACAAGCCGTCCCCTCGCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCTCCCTTCCAGGGTGAgggct                                                                                                                                          <  2:1234120/150‑3 (MQ=18)
           taataaGGTTACGGGCGTAATGACAAATGCGACAAAGCGCACAATCCGTCCCGTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCTCCCTTCCAGGGTGAgggct                                                                                                                                          <  2:1234122/150‑3 (MQ=25)
           taataaGGTTAAGGGCGTAATGACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCTCCCTTCCAGGGTGAgggct                                                                                                                                          <  2:429883/150‑3 (MQ=25)
              taaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAg                                                                                                                                       >  1:1225948/1‑150 (MQ=255)
                 ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAAc                                                                                                                                    >  1:936514/1‑150 (MQ=255)
                  gTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                                                                                                   >  2:831501/1‑150 (MQ=255)
                       gggCTTAATTACAAATTCCAAAACGCGCACAACCCGTCCCCTCGCCCCTTTGGGTTGTGGTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                              <  1:1120281/150‑1 (MQ=25)
                       gggCTTAATGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCGTTTGGGGTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                              <  1:100161/150‑1 (MQ=255)
                       gggCGTAATGACAAATCCGACAACGCGCACAATCCGTCCCCTGGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                              <  1:1384971/150‑1 (MQ=255)
                             aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                                                                                                        >  1:780415/1‑150 (MQ=255)
                             aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAGTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                                                                                                        >  2:1284285/1‑150 (MQ=255)
                               tGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATca                                                                                                                      <  2:235163/150‑1 (MQ=255)
                               tGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATca                                                                                                                      >  2:475160/1‑150 (MQ=255)
                                gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcac                                                                                                                     >  2:246290/1‑150 (MQ=255)
                                gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAACCCAGCCCTCACCCTGGCAGGGAGAGGGGGCCGACCGGCGCAGGACCTccc                                                                                                                     >  2:246236/1‑148 (MQ=25)
                                 aaaaaTTCTACAATGCGAAAAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                                                                                                    <  1:967735/148‑1 (MQ=255)
                                 aaaaaTCCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGTTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                                                                                                    <  1:967730/148‑1 (MQ=255)
                                    aaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTg                                                                                                                 >  1:686603/1‑150 (MQ=255)
                                      ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACGGGTGCGAACATGAACCCTCACCCCAGCCCTCACCCGGGACGGGAGAGGGGGCCGAACGGCGCAGGCCCTCACCTTgcg                                                                                                               >  2:1264285/1‑150 (MQ=25)
                                       tCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc                                                                                                              >  2:489640/1‑150 (MQ=255)
                                       tCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc                                                                                                              >  2:489662/1‑150 (MQ=255)
                                            aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTAGGGGGAGGGGAACCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAGAGCGCTTATg                                                                                                         <  2:718532/150‑1 (MQ=255)
                                             aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                                                                                                        <  2:51424/150‑1 (MQ=255)
                                               gcgcACAATCGGGCCCCTCGCCCTTTGGGGGTGTGGTTTCGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGa                                                                                                      <  2:140404/150‑1 (MQ=255)
                                                cgcACAAGCCGTCCCCTCGCCCCTTTGGGGATAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                                                                                                     <  1:740475/150‑1 (MQ=255)
                                                  cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACCTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                                                                                                   >  1:695543/1‑150 (MQ=255)
                                                  cacaATCCGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                                                                                                   >  2:581534/1‑150 (MQ=255)
                                                   acaATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTTGGTTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                                                                                                  <  2:69681/150‑1 (MQ=255)
                                                     aaTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                                                                                                <  2:983363/150‑1 (MQ=255)
                                                       tCCGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                                                                                              <  1:581534/150‑1 (MQ=255)
                                                         cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaat                                                                                            >  1:780451/1‑150 (MQ=255)
                                                          gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata                                                                                           >  2:383669/1‑150 (MQ=255)
                                                          gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAGTGCGCTTATGCGAATCCATCaata                                                                                           >  2:383668/1‑150 (MQ=255)
                                                           tCCCCTCGCCCGTTTGGGGTGTGGGTTTGGGGGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaataa                                                                                          <  1:604826/150‑1 (MQ=255)
                                                              ccTCGCCCTTTTGGGTTGTGGTTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc                                                                                       <  1:522763/150‑1 (MQ=255)
                                                                            ggtgtgggTTCGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                                                                         <  2:1103634/145‑1 (MQ=255)
                                                                            ggAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                                                                         <  2:926990/150‑1 (MQ=255)
                                                                            ggAGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                                                                         <  1:544496/150‑1 (MQ=255)
                                                                             gagaGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTg                                                                        >  1:633411/1‑150 (MQ=255)
                                                                               gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTa                                                                      <  2:1163206/150‑1 (MQ=255)
                                                                                  ggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAAcc                                                                   >  2:1381828/1‑150 (MQ=255)
                                                                                         gTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCa                                                            <  1:1307739/150‑1 (MQ=255)
                                                                                           gcggggAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGa                                                          <  1:1063214/148‑1 (MQ=255)
                                                                                           gAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGa                                                          <  1:867211/150‑1 (MQ=255)
                                                                                            aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt                                                         <  1:189134/150‑1 (MQ=255)
                                                                                             ggggAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATg                                                        <  1:741622/150‑1 (MQ=255)
                                                                                                 aaCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAc                                                    <  2:512858/150‑1 (MQ=255)
                                                                                                 aaCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAc                                                    <  2:512846/150‑1 (MQ=255)
                                                                                                    agccagcACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAcccc                                                 >  1:470349/1‑150 (MQ=255)
                                                                                                               gTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGAt                                      >  2:642765/1‑150 (MQ=255)
                                                                                                                tGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGAtt                                     >  1:1386768/1‑150 (MQ=255)
                                                                                                                tGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGAtg                                     >  1:290325/1‑149 (MQ=255)
                                                                                                                tGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAAGCAACCCCGAGTTGCAGAtt                                     >  1:1386770/1‑150 (MQ=255)
                                                                                                                   gAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCACCCCCGAGTTGCAGATTGCg                                  >  1:1166121/1‑150 (MQ=255)
                                                                                                                    aaCATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGt                                 >  1:307740/1‑150 (MQ=255)
                                                                                                                     aCATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAGTGCGCGTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGtt                                >  2:145473/1‑150 (MQ=255)
                                                                                                                       aTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAc                              <  2:215678/150‑1 (MQ=255)
                                                                                                                       aTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAc                              <  1:271394/150‑1 (MQ=255)
                                                                                                                       aTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAc                              <  1:271400/150‑1 (MQ=255)
                                                                                                                        ttAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAcc                             <  1:55751/150‑1 (MQ=255)
                                                                                                                          aaCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTc                           <  2:340540/150‑1 (MQ=255)
                                                                                                                          aaCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTc                           <  2:359582/150‑1 (MQ=255)
                                                                                                                           aCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCa                          >  2:745872/1‑150 (MQ=255)
                                                                                                                           aCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGGTACCTCa                          >  2:719541/1‑150 (MQ=255)
                                                                                                                             ccTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCCGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAag                        >  2:1086007/1‑150 (MQ=255)
                                                                                                                             ccTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCacg                        >  2:1086010/1‑148 (MQ=255)
                                                                                                                              cTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAaga                       >  1:545459/1‑150 (MQ=255)
                                                                                                                                   cccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGa                  <  1:1073688/150‑1 (MQ=255)
                                                                                                                                   cccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGa                  >  1:781812/1‑150 (MQ=255)
                                                                                                                                   cccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGa                  <  1:1073690/150‑1 (MQ=255)
                                                                                                                                       gCCCCCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATa              <  1:145552/150‑1 (MQ=255)
                                                                                                                                        cccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAg             >  1:1387629/1‑150 (MQ=255)
                                                                                                                                          cTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTg           >  1:755060/1‑150 (MQ=255)
                                                                                                                                           tCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTgc          >  2:592681/1‑150 (MQ=255)
                                                                                                                                           tCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAGTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCGCAAGAGTTGACATAGTgc          >  1:593753/1‑150 (MQ=255)
                                                                                                                                             aCCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTgcgc        >  2:541853/1‑150 (MQ=255)
                                                                                                                                             aCCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTgcgc        >  2:541899/1‑150 (MQ=255)
                                                                                                                                             aCCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAATGTAACCAAAACCAGATGCAACACCGAGTTGCAGATTGCGATACCTCAAGAGTTGACATAGTgcgc        >  2:541836/1‑150 (MQ=255)
                                                                                                                                              cccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCATATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGcgcg       >  2:836108/1‑150 (MQ=255)
                                                                                                                                                  ggAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCTGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTg   <  2:1312041/150‑1 (MQ=255)
                                                                                                                                                  ggAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTg   <  2:1312081/150‑1 (MQ=255)
                                                                                                                                                  ggAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTg   <  2:1312042/150‑1 (MQ=255)
                                                                                                                                                   gAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGc  <  1:154969/150‑1 (MQ=255)
                                                                                                                                                   gAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGc  <  1:1024673/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGC  >  NC_000913/624518‑624814

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