Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,313,903 G→T 23.0% L12L (CTG→CTT yjdP → putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,9030GT23.0% 258.9 / 70.5 109L12L (CTG→CTTyjdPputative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (43/41);  new base T (14/11);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 109 total aligned reads displayed.

AACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAG  >  NC_000913/4313756‑4314049
                                                                                                                                                   |                                                                                                                                                  
aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGc                                                                                                                                                   >  3:593861/1‑149 (MQ=255)
   gcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCa                                                                                                                                                >  3:666809/1‑149 (MQ=255)
      tGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGc                                                                                                                                             <  3:466403/149‑1 (MQ=255)
        gAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCt                                                                                                                                            <  2:677263/148‑1 (MQ=255)
        gAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCt                                                                                                                                            >  1:677263/1‑148 (MQ=255)
           gCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTcc                                                                                                                                        >  4:640932/1‑149 (MQ=255)
           gCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTcc                                                                                                                                        <  4:342808/149‑1 (MQ=255)
                   ggCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                                                                                                                                >  1:606062/1‑149 (MQ=255)
                     cAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                                                                                                                              >  3:486221/1‑149 (MQ=255)
                       aTAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCTCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGcc                                                                                                                            >  3:74033/1‑149 (MQ=255)
                           agcagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAgg                                                                                                                        <  3:50378/149‑1 (MQ=255)
                            gcagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGc                                                                                                                       >  2:591274/1‑149 (MQ=255)
                                      aaTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                                                                                             <  1:525715/149‑1 (MQ=255)
                                       aTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg                                                                                                            >  2:478440/1‑149 (MQ=255)
                                           agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacga                                                                                                        <  1:223018/149‑1 (MQ=255)
                                           agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacga                                                                                                        <  3:167608/149‑1 (MQ=255)
                                             aTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCgg                                                                                                                              <  4:309915/125‑1 (MQ=255)
                                             aTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCgg                                                                                                                              >  3:309915/1‑125 (MQ=255)
                                                 aTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                                                                                             >  4:802542/1‑138 (MQ=255)
                                                 aTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                                                                                             <  3:802542/138‑1 (MQ=255)
                                                 aTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                                                                                  <  1:638799/149‑1 (MQ=255)
                                                  tCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                                                                 >  1:742406/1‑149 (MQ=255)
                                                  tCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                                                                 >  4:634147/1‑149 (MQ=255)
                                                   ccACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                                                                <  3:358343/149‑1 (MQ=255)
                                                     actactTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac                                                                                                          >  1:929581/1‑137 (MQ=255)
                                                     actactTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac                                                                                                          <  2:929581/137‑1 (MQ=255)
                                                     actactTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                                                                              <  3:262026/149‑1 (MQ=255)
                                                          ttGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                                                                         >  3:139809/1‑149 (MQ=255)
                                                              aTGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                                                                     <  3:902370/149‑1 (MQ=255)
                                                               tGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                                                                    >  4:281754/1‑149 (MQ=255)
                                                                 gATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                                                                  >  4:414353/1‑149 (MQ=255)
                                                                 gATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                                                                  >  3:350953/1‑149 (MQ=255)
                                                                     aGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                                                              <  4:849133/149‑1 (MQ=255)
                                                                      gTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                                                             <  2:734232/149‑1 (MQ=255)
                                                                       tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                                                                            >  3:72366/1‑149 (MQ=255)
                                                                       tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                                                                            >  4:283775/1‑149 (MQ=255)
                                                                       tAATGATTAATGTGGACAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                                                                            <  4:621208/149‑1 (MQ=255)
                                                                        aaTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTa                                                                           >  4:383993/1‑149 (MQ=255)
                                                                         aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                                                                          <  4:488966/149‑1 (MQ=255)
                                                                           gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                                                                        >  2:847477/1‑149 (MQ=255)
                                                                             ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                                                                      >  1:267612/1‑149 (MQ=255)
                                                                                aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                                                                   <  3:283775/149‑1 (MQ=255)
                                                                                aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                                                                   >  1:818863/1‑149 (MQ=255)
                                                                                 tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                                                                          <  2:516199/141‑1 (MQ=255)
                                                                                 tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                                                                          >  1:516199/1‑141 (MQ=255)
                                                                                     gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                                                              >  4:189975/1‑149 (MQ=255)
                                                                                     gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                                                              <  4:666809/149‑1 (MQ=255)
                                                                                        agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc                                                           <  2:233259/149‑1 (MQ=255)
                                                                                        agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                                                                    <  2:601000/140‑1 (MQ=255)
                                                                                         gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt                                                          <  2:546062/149‑1 (MQ=255)
                                                                                         gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt                                                          <  3:189975/149‑1 (MQ=255)
                                                                                            tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCTGCAGGCGATTAATGACGCTTATAACCCTGAc                                                       >  1:148452/1‑149 (MQ=255)
                                                                                            tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc                                                       >  3:163919/1‑149 (MQ=255)
                                                                                            tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc                                                       >  1:770749/1‑149 (MQ=255)
                                                                                            tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc                                                       <  3:634147/149‑1 (MQ=255)
                                                                                               cTTTTTGAGGTTGGCTAATGAAACGCTGTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt                                                    >  2:790030/1‑149 (MQ=255)
                                                                                                 ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg                                                  <  4:816059/149‑1 (MQ=255)
                                                                                                     gAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCga                                              <  4:30561/149‑1 (MQ=255)
                                                                                                      aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                                                              >  4:397448/1‑132 (MQ=255)
                                                                                                      aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                                                              <  3:397448/132‑1 (MQ=255)
                                                                                                       ggTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCgaga                                            <  2:267612/149‑1 (MQ=255)
                                                                                                       ggTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCgaga                                            <  4:350953/149‑1 (MQ=255)
                                                                                                         ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACt                                          <  3:640932/149‑1 (MQ=255)
                                                                                                          tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTa                                         <  4:593861/149‑1 (MQ=255)
                                                                                                             cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGa                                      >  2:653840/1‑149 (MQ=255)
                                                                                                             cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGa                                      >  1:375261/1‑149 (MQ=255)
                                                                                                              tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGaa                                     >  3:111301/1‑149 (MQ=255)
                                                                                                              tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGaa                                     >  2:123977/1‑149 (MQ=255)
                                                                                                                aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGa                                   >  2:431150/1‑149 (MQ=255)
                                                                                                                 tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGAt                                  <  1:801915/149‑1 (MQ=255)
                                                                                                                    aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTcg                               >  2:621072/1‑149 (MQ=255)
                                                                                                                     aCGCTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTcgc                              <  4:111301/149‑1 (MQ=255)
                                                                                                                     aCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTcgc                              <  2:606062/149‑1 (MQ=255)
                                                                                                                         tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGcga                          >  2:315483/1‑149 (MQ=255)
                                                                                                                         tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGcga                          >  4:385031/1‑149 (MQ=255)
                                                                                                                             cGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTcgcgc                            >  1:27550/1‑143 (MQ=255)
                                                                                                                             cGCTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTcgcgc                            <  2:27550/143‑1 (MQ=255)
                                                                                                                               cTATTTCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACgg                    <  3:383993/149‑1 (MQ=255)
                                                                                                                                 aTTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCt                  <  4:907620/149‑1 (MQ=255)
                                                                                                                                  tttCTTTTATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTg                 >  1:437061/1‑149 (MQ=255)
                                                                                                                                      ttttATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCga                                              >  3:827791/1‑116 (MQ=255)
                                                                                                                                      ttttATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCga                                              <  4:827791/116‑1 (MQ=255)
                                                                                                                                      ttttATTCACCCTGCTCACGCTGTCCACCTTTCCGGCCCATGCCGATATTCTCGACGAAACCACTGGCAATATTCAGCAGGCGATTAATGCCGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAg             >  1:359837/1‑149 (MQ=255)
                                                                                                                                      ttttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAg             >  1:221707/1‑149 (MQ=255)
                                                                                                                                      ttttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAg             >  1:860448/1‑149 (MQ=255)
                                                                                                                                       tttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGc            <  2:810141/149‑1 (MQ=255)
                                                                                                                                        ttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCg           >  4:734084/1‑149 (MQ=255)
                                                                                                                                         tATTCCCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGt          <  4:140661/149‑1 (MQ=255)
                                                                                                                                         tATTCACCCTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGt          <  3:75904/149‑1 (MQ=255)
                                                                                                                                          aTTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTg         >  2:860397/1‑149 (MQ=255)
                                                                                                                                           ttCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGa                                      <  2:827774/119‑1 (MQ=255)
                                                                                                                                           ttCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGa                                      >  1:827774/1‑119 (MQ=255)
                                                                                                                                            tCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGaa       <  2:375261/149‑1 (MQ=255)
                                                                                                                                            tCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGaa       <  1:431150/149‑1 (MQ=255)
                                                                                                                                              aCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGt     >  3:714888/1‑149 (MQ=255)
                                                                                                                                                ccTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                                                                   <  2:392365/85‑1 (MQ=255)
                                                                                                                                                ccTTCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                                                                   >  1:392365/1‑85 (MQ=255)
                                                                                                                                                ccTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGa   >  2:418677/1‑149 (MQ=255)
                                                                                                                                                ccTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGa   <  2:860448/149‑1 (MQ=255)
                                                                                                                                                 cTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAg  <  3:917869/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAG  >  NC_000913/4313756‑4314049

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: