Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,515,328 T→G 5.5% E450D (GAA→GAC) ‡ fecA ← TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,3280TG 5.5% 393.4 / 3.7 113E450D (GAA→GAC) ‡fecATonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Reads supporting (aligned to +/- strand):  ref base T (49/57);  new base G (5/2);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGG  >  NC_000913/4515186‑4515475
                                                                                                                                              |                                                                                                                                                   
ccGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATg                                                                                                                                                <  1:77353/148‑1 (MQ=255)
   gCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTACCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGtt                                                                                                                                            >  4:623486/1‑149 (MQ=255)
     gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTg                                                                                                                                           <  1:232924/148‑1 (MQ=255)
      tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGat                                                                                                                                         <  3:292006/149‑1 (MQ=255)
      tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGat                                                                                                                                         <  4:870781/149‑1 (MQ=255)
      tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATGAAATGCGttac                                                                                                                                          <  1:859684/148‑3 (MQ=255)
       gATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGatt                                                                                                                                        <  4:431162/149‑1 (MQ=255)
       gATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGat                                                                                                                                         <  2:141359/148‑1 (MQ=255)
        aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGatt                                                                                                                                        <  1:913705/148‑1 (MQ=255)
          ggTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattca                                                                                                                                      <  1:177873/148‑1 (MQ=255)
          ggTCCAGTCGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattca                                                                                                                                      <  1:547816/148‑1 (MQ=255)
           gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcatt                                                                                                                                    <  3:430891/149‑1 (MQ=255)
           gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcat                                                                                                                                     <  2:880483/148‑1 (MQ=255)
               aGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                                                                                <  3:261383/149‑1 (MQ=255)
               aGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaa                                                                                                                                 >  1:145126/1‑148 (MQ=255)
                   gCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGc                                                                                                                            <  3:341954/149‑1 (MQ=255)
                     cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGa                                                                                                                          >  3:486180/1‑149 (MQ=255)
                            gATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgccca                                                                                                                   >  3:299764/1‑149 (MQ=255)
                             aTTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTTATGCGTTGATTCATTCAAATAGCGATAgcccac                                                                                                                  >  4:839620/1‑149 (MQ=255)
                               tttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgc                                                                                                                <  4:467176/149‑1 (MQ=255)
                                      tCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                          >  2:187531/1‑148 (MQ=255)
                                      tCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                          >  2:699041/1‑148 (MQ=255)
                                      tCCAGATACCAGGCGTACGCCTCGGTGCCGGAACGCGTCTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAAAAGCGATAGCCCACGCCCACTTc                                                                                                         >  4:462249/1‑149 (MQ=255)
                                          gATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTg                                                                                                      <  2:843291/148‑1 (MQ=255)
                                          gATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTg                                                                                                      >  2:5916/1‑148 (MQ=255)
                                           aTACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTgc                                                                                                     >  1:820127/1‑148 (MQ=255)
                                                aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               >  3:173722/1‑149 (MQ=255)
                                                aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               >  3:767648/1‑149 (MQ=255)
                                                 ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               <  2:820127/148‑1 (MQ=255)
                                                 ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               >  2:461401/1‑148 (MQ=255)
                                                       cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcat                                                                                                                                     <  1:850291/104‑1 (MQ=255)
                                                       cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcat                                                                                                                                     >  2:850291/1‑104 (MQ=255)
                                                        gcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGAt                                                                                       <  3:45461/149‑1 (MQ=255)
                                                        gcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGAt                                                                                       <  4:224925/149‑1 (MQ=255)
                                                         ccTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATc                                                                                      >  3:378377/1‑149 (MQ=255)
                                                          cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCa                                                                                     <  3:472812/149‑1 (MQ=255)
                                                           tCGGTGCCGGAACGCGTTTCGCGGTCGTAAGGGCTTGAGCCTGACGGCAGCTGCCCGCTGCTGGTGGCGCTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGACCAGGCCCACTTCGTGCGCGGATGGGCCGATCAt                                                                                    >  4:590180/1‑149 (MQ=255)
                                                            cGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCAt                                                                                    >  1:685515/1‑148 (MQ=255)
                                                            cGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATa                                                                                   <  3:374060/149‑1 (MQ=255)
                                                              gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCGGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGCTAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                                 <  3:623486/149‑1 (MQ=255)
                                                              gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                                 <  4:414186/149‑1 (MQ=255)
                                                              gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                                 <  4:5589/149‑1 (MQ=255)
                                                              gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                                 <  3:100361/149‑1 (MQ=255)
                                                              gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                                  <  2:692491/148‑1 (MQ=255)
                                                                gCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGa                                                                               >  3:103280/1‑149 (MQ=255)
                                                                   ggAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCt                                                                            <  4:301156/149‑1 (MQ=255)
                                                                    gAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                           >  3:126355/1‑149 (MQ=255)
                                                                     aaCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                              >  2:927619/1‑129 (MQ=255)
                                                                     aaCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                              <  1:927619/129‑1 (MQ=255)
                                                                     aaCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTgg                                                                          <  3:819333/149‑1 (MQ=255)
                                                                        gcgTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCt                                                                        >  2:84099/1‑148 (MQ=255)
                                                                            aTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAg                                                                    <  2:40338/148‑1 (MQ=255)
                                                                              cgcgGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg                                                                  <  1:84099/148‑1 (MQ=255)
                                                                              cgcgGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGt                                                                 <  3:439418/149‑1 (MQ=255)
                                                                                 ggTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGc                                                              <  3:847774/149‑1 (MQ=255)
                                                                                     gTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCaa                                                          >  3:695972/1‑149 (MQ=255)
                                                                                     gTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCaa                                                          >  4:746977/1‑149 (MQ=255)
                                                                                         gggCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                       <  2:759234/148‑1 (MQ=255)
                                                                                           gCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACcg                                                     <  2:243043/148‑1 (MQ=255)
                                                                                             ttGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTgg                                                                          <  4:165097/125‑1 (MQ=255)
                                                                                             ttGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCgcgc                                                  >  4:740220/1‑149 (MQ=255)
                                                                                                 gCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               <  4:65887/100‑1 (MQ=255)
                                                                                                 gCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                               >  3:65887/1‑100 (MQ=255)
                                                                                                  ccGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAccc                                              >  1:644215/1‑148 (MQ=255)
                                                                                                  ccGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAccc                                              <  2:730852/148‑1 (MQ=255)
                                                                                                    ggACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagt                                           >  3:517404/1‑149 (MQ=255)
                                                                                                    ggACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagt                                           <  3:924435/149‑1 (MQ=255)
                                                                                                     gACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagta                                          <  4:404714/149‑1 (MQ=255)
                                                                                                       cGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagtag                                         <  1:521410/148‑1 (MQ=255)
                                                                                                        ggCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGtt                                       <  3:696009/149‑1 (MQ=255)
                                                                                                          cAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTa                                      >  1:393810/1‑148 (MQ=255)
                                                                                                           aaCTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAc                                     >  2:883231/1‑148 (MQ=255)
                                                                                                                    ctgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTg                                                                                                      >  3:917687/1‑74 (MQ=255)
                                                                                                                    ctgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCgaga                            >  2:206468/1‑148 (MQ=255)
                                                                                                                    ctgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                           <  3:707455/149‑1 (MQ=255)
                                                                                                                     tgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                           <  2:836755/148‑1 (MQ=255)
                                                                                                                        tggtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGt                        >  2:668845/1‑148 (MQ=255)
                                                                                                                        tggtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTg                       >  3:25093/1‑149 (MQ=255)
                                                                                                                         ggtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGa                      <  3:746977/149‑1 (MQ=255)
                                                                                                                            ggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGc                   <  3:105249/149‑1 (MQ=255)
                                                                                                                             gCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                           >  2:348041/1‑140 (MQ=255)
                                                                                                                             gCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                           <  1:348041/140‑1 (MQ=255)
                                                                                                                             gCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGc                   >  1:667498/1‑148 (MQ=255)
                                                                                                                               ggTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGt                 <  1:550304/148‑1 (MQ=255)
                                                                                                                                gtgtAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtt                <  2:811454/148‑1 (MQ=255)
                                                                                                                                   tagtaACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc             >  1:684476/1‑148 (MQ=255)
                                                                                                                                    atcaaCGCATGAAATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCt           >  4:442272/4‑149 (MQ=255)
                                                                                                                                      taACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTg         >  4:431435/1‑149 (MQ=255)
                                                                                                                                      caaCGCATGAAATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCtt          >  2:724545/2‑148 (MQ=255)
                                                                                                                                      caaCGCATGAAATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCtt          >  1:472834/2‑148 (MQ=255)
                                                                                                                                       aaCGCATGAAATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                           <  3:492639/82‑1 (MQ=255)
                                                                                                                                       aaCGCATGAAATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                           >  4:492639/1‑82 (MQ=255)
                                                                                                                                        aCGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCt       >  3:446755/1‑149 (MQ=255)
                                                                                                                                        aCGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCt       >  4:786081/1‑149 (MQ=255)
                                                                                                                                         cGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCt       >  2:79754/1‑148 (MQ=255)
                                                                                                                                         cGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGGGATGCGTTTGCCTTGCTc      >  4:887183/1‑149 (MQ=255)
                                                                                                                                            aTTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCa    <  1:309713/148‑1 (MQ=255)
                                                                                                                                            aTTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCa    <  1:206468/148‑1 (MQ=255)
                                                                                                                                            aTTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAg   <  4:103280/149‑1 (MQ=255)
                                                                                                                                             tttCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAgg  <  4:767648/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGG  >  NC_000913/4515186‑4515475

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: