Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,307,304 G→T 5.9% L438I (CTT→ATT)  yojI ← microcin J25 efflux ABC transporter permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,307,3040GT 5.9% 228.4 / 3.2 68L438I (CTT→ATT) yojImicrocin J25 efflux ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTCGCC  >  NC_000913/2307156‑2307451
                                                                                                                                                    |                                                                                                                                                   
tCACGACGGAAGTGTAGATCATGATCCTCCGCCAATTAATCAATAATGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAg                                                                                                                                                     >  2:34978/1‑149 (MQ=255)
tCACGACGGAAGTGAGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTGTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAg                                                                                                                                                     >  8:71108/1‑149 (MQ=255)
  aagacgGAAGTGTGGATCCTGATCCGCCGCCCATTCCTCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCCCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCt                                                                                                                                                   <  6:143959/147‑1 (MQ=255)
    gacgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTa                                                                                                                                                 >  7:160466/1‑149 (MQ=255)
            tgtgGATCCTGATCCGCCGCCCATTCATCGATCTGGATAATATCTCGTTCTTCTGACAGCGCCAGCAACAGCGCCAGGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCTTTGCTTAACTCAAGCTTATGAGCCa                                                                                                                                          <  4:189588/148‑1 (MQ=255)
                 aTCCTGATCCGCCTCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTa                                                                                                                                                 <  1:269365/136‑1 (MQ=255)
                 aTCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGACACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTa                                                                                                                                                 >  2:269365/1‑136 (MQ=255)
                       aTCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTgcg                                                                                                                              >  7:310812/1‑149 (MQ=255)
                          cgccgcCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCAtttt                                                                                                                                      >  4:63139/1‑138 (MQ=255)
                          cgccgcCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCAtttt                                                                                                                                      <  3:63139/138‑1 (MQ=255)
                             cgccCATTCATCCAGCAGGATAATATCGCGTTGTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCgccagcc                                                                                                                        <  7:71108/149‑1 (MQ=255)
                                  aTTCATCCAGGATGATAATATCGCGTTATTCTTCCAGCGACATCAACAGCGCCACGCGTTTTTTCTCCCCTTTTGATAACTTCAGGTTAACACTACTCCCGTTGCTTAACTCAAGCTTCGGAGCCAATTTCAGCGGCgccagccccatc                                                                                                                   <  5:246467/149‑6 (MQ=255)
                                   ttCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAATCTTa                                                                                                                                                 <  3:101553/118‑1 (MQ=255)
                                   ttCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAATCTTa                                                                                                                                                 >  4:101553/1‑118 (MQ=255)
                                   ttCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAATCTTa                                                                                                                                                 <  3:96785/118‑1 (MQ=255)
                                   ttCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAATCTTa                                                                                                                                                 >  4:96785/1‑118 (MQ=255)
                                           gcagGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGt                                                                                                          >  2:103853/1‑149 (MQ=255)
                                            cagGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGtt                                                                                                         <  8:97615/149‑1 (MQ=255)
                                               gATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGATCCATTTTCAGCTGCGGCAGCCACTGCTCATCCAGTGGCg                                                                                                      >  3:92555/1‑149 (MQ=255)
                                                ataataTCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCgg                                                                                                     <  2:354340/149‑1 (MQ=255)
                                                ataataTCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCGTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCgg                                                                                                     >  2:253373/1‑149 (MQ=255)
                                                 taataTCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCggg                                                                                                    <  2:335172/149‑1 (MQ=255)
                                                      tCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTg                                                                                                        >  7:205025/1‑140 (MQ=255)
                                                      tCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTg                                                                                                        <  8:205025/140‑1 (MQ=255)
                                                      tCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGc                                                                                               >  7:10935/1‑149 (MQ=255)
                                                      tCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGc                                                                                               <  1:34978/149‑1 (MQ=255)
                                                      tCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCGTTTGATAACGTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTGCTCAACCAGTTGCGGGTTAGc                                                                                               >  3:334260/1‑149 (MQ=255)
                                                             cttctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTa                                                                                        >  1:76413/1‑149 (MQ=255)
                                                             cttctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTa                                                                                        >  8:16075/1‑149 (MQ=255)
                                                                ctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTAc                                                                                       <  7:169348/147‑1 (MQ=255)
                                                                   ccagcgccagcAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTcc                                                                                  >  7:341248/1‑149 (MQ=255)
                                                                      gcgccagcAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGc                                                                               <  7:16075/149‑1 (MQ=255)
                                                                               aCAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGt                                                                      >  7:263160/1‑149 (MQ=255)
                                                                                     ccACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGtt                                                                     <  7:209196/144‑1 (MQ=255)
                                                                                     ccACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGtt                                                                     >  8:209196/1‑144 (MQ=255)
                                                                                            ttttttCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGc                                                          >  2:171343/1‑148 (MQ=255)
                                                                                             tttttCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGc                                                          >  1:348046/1‑147 (MQ=255)
                                                                                              ttttCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCa                                                        <  8:160466/148‑1 (MQ=255)
                                                                                              ttttCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCa                                                        <  2:348046/148‑1 (MQ=255)
                                                                                                 tCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAACCAGCCAGAc                                                     <  1:103853/148‑1 (MQ=255)
                                                                                                          ttGATAACTTCAGGTTAACCATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTTCGGGTTAGCGGGTTTACACTCCGGCCCAAGCACTAGATCAAACAGCCAGACATCGGTaaa                                            >  3:202078/1‑148 (MQ=255)
                                                                                                            gATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTc                                                                                   >  7:284391/1‑107 (MQ=255)
                                                                                                            gATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTc                                                                                   <  8:284391/107‑1 (MQ=255)
                                                                                                            gATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAaca                                          <  8:263160/148‑1 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGttt                                                                                         <  4:123527/95‑1 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAAc                                                             >  4:277758/1‑123 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAAc                                                             <  3:277758/123‑1 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAacac                                         >  4:59499/1‑143 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAacac                                         <  3:59499/143‑1 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCg                                    >  7:209227/1‑148 (MQ=255)
                                                                                                                  ttCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCg                                    >  5:343845/1‑148 (MQ=255)
                                                                                                                     aGGTTAACACTACTCCCGAGGCTACACTCACGCTTCTGCGACATTATCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGCTGATCAAACAGCCAGACATCGGTAAACACTGCCGaaaa                                <  4:229040/149‑1 (MQ=255)
                                                                                                                      ggTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGAc                                                     >  3:237198/1‑127 (MQ=255)
                                                                                                                      ggTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGAc                                                     <  4:237198/127‑1 (MQ=255)
                                                                                                                           aCAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGttt                          >  7:5005/1‑149 (MQ=255)
                                                                                                                                  gCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGttt                                                                                         <  2:153446/79‑1 (MQ=255)
                                                                                                                                  gCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGttt                                                                                         >  1:153446/1‑79 (MQ=255)
                                                                                                                                         gCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTaaa                                            <  3:136141/117‑1 (MQ=255)
                                                                                                                                         gCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTaaa                                            >  4:136141/1‑117 (MQ=255)
                                                                                                                                              aCTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGgttgt       >  4:135626/1‑149 (MQ=255)
                                                                                                                                               cTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCCGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACGGCCGAAAACAGTTTGCGATAATCTTCCGgttgtt      >  1:189083/1‑149 (MQ=255)
                                                                                                                                               cTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCa                                                        >  4:155818/1‑99 (MQ=255)
                                                                                                                                               cTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCa                                                        <  3:155818/99‑1 (MQ=255)
                                                                                                                                                tCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCTGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTc     >  3:119183/1‑149 (MQ=255)
                                                                                                                                                tCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGaaaa                                <  4:264130/122‑1 (MQ=255)
                                                                                                                                                tCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGgtt         <  4:76397/145‑1 (MQ=255)
                                                                                                                                                tCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGgtt         >  3:76397/1‑145 (MQ=255)
                                                                                                                                                  aaGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCTGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTc     <  4:119183/147‑1 (MQ=255)
                                                                                                                                                  aaGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTcgc   >  8:341348/1‑149 (MQ=255)
                                                                                                                                                   aGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTcgcc  >  1:114551/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCGTTGCTTAACTCAAGCTTATGAGCCATTTTCAGCTGCGCCAGCCACTTCTCAACCAGTTGCGGGTTAGCGGGTTTACCCTCCGGCCCCAGCAGTTGATCAAACAGCCAGACATCGGTAAACACTGCCGAAAACAGTTTGCGATAATCTTCCGGTTGTTCGCC  >  NC_000913/2307156‑2307451

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