Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,257,754 G→A 6.3% V5I (GTC→ATC)  dinF → oxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,257,7540GA 6.3% 403.0 / 10.3 127V5I (GTC→ATC) dinFoxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base G (56/63);  new base A (3/5);  total (59/68)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGG  >  NC_000913/4257616‑4257895
                                                                                                                                          |                                                                                                                                              
ttGCCAGAAAATAGCGAGGTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTCTTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGgtttg                                                                                                                                      <  3:80110/149‑1 (MQ=255)
   ccAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGgtttgttt                                                                                                                                   <  5:337538/149‑1 (MQ=255)
    cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGtttt                                                                                                                                  >  3:143182/1‑149 (MQ=255)
    cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGtttt                                                                                                                                  >  8:176796/1‑149 (MQ=255)
      gAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGttttt                                                                                                                                 >  4:292213/1‑148 (MQ=255)
                aGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGGTTTTCATCTCTCtt                                                                                                                      >  3:63765/1‑148 (MQ=255)
                  tttAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtcat                                                                                                                   <  8:316848/149‑1 (MQ=255)
                     aaaCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGCTGATTCGCAACGGCGACTGGCTGTAACATATCTCAGAGACCGCGATTCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAg                                                                                                                >  6:330030/1‑149 (MQ=255)
                      aaCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtcatca                                                                                                                 <  5:384107/147‑1 (MQ=255)
                       aCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtc                                                                                                                     >  1:418662/1‑142 (MQ=255)
                       aCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGc                                                                                                              <  4:143182/149‑1 (MQ=255)
                       aCCAATTGTCGTTGACCTTCGTCAGCAGAGCTACACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGc                                                                                                              >  6:234938/1‑149 (MQ=255)
                              gTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCt                                                                                                        <  6:98833/148‑1 (MQ=255)
                                 gTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCAt                                                                                                    <  5:234938/149‑1 (MQ=255)
                                 gTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCAt                                                                                                    >  3:397517/1‑149 (MQ=255)
                                  ttGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCa                                                                                                                              <  3:23321/122‑1 (MQ=255)
                                  ttGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCa                                                                                                                              >  4:23321/1‑122 (MQ=255)
                                  ttGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATg                                                                                                   >  7:12020/1‑149 (MQ=255)
                                  ttGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATg                                                                                                   <  1:296927/149‑1 (MQ=255)
                                    gACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAAACCGCGACGCCAGGCGGCATCGCGgtct                                                                                                                                       <  5:10793/112‑3 (MQ=38)
                                    gACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAAACCGCGACGCCAGGCGGCATCGCGgtct                                                                                                                                       >  6:10793/1‑110 (MQ=38)
                                         tCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCa                                                                                                <  8:380239/145‑1 (MQ=255)
                                          cGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTATCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTc                                                                                             <  6:448645/147‑1 (MQ=255)
                                           gTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCt                                                                                           <  2:304685/148‑1 (MQ=255)
                                           gTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCt                                                                                           >  5:5784/1‑148 (MQ=255)
                                           gTCAGCAGAGCTTCACCATTGAAGGGCTGGGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtcatca                                                                                                                 >  3:283160/1‑125 (MQ=255)
                                           gTCAGCAGAGCTTCACCATTGAAGGGCTGGGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtcatca                                                                                                                 <  4:283160/125‑1 (MQ=255)
                                               gcagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACtt                                                                                      <  3:84283/149‑1 (MQ=255)
                                                cagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGgt                                                                                                                                         <  8:324310/98‑1 (MQ=255)
                                                cagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGgt                                                                                                                                         >  7:324310/1‑98 (MQ=255)
                                                  gagCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCtt                                                                                                            >  8:260447/1‑124 (MQ=255)
                                                   agCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATc                                                                                  <  3:63211/149‑1 (MQ=255)
                                                     cTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTTTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTc                                                                                     >  4:217091/1‑144 (MQ=255)
                                                     cTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTc                                                                                     <  3:217091/144‑1 (MQ=255)
                                                        cACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGAt                                                                              >  6:260526/1‑148 (MQ=255)
                                                         aCCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATa                                                                             <  4:397517/148‑1 (MQ=255)
                                                          ccATTGAAGGGCAGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATa                                                                             >  6:334405/1‑147 (MQ=255)
                                                             ttGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCa                                                                        >  8:444432/1‑149 (MQ=255)
                                                               gAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCAct                                                                      >  5:166834/1‑149 (MQ=255)
                                                                    gCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGc                                                                 <  5:20712/149‑1 (MQ=255)
                                                                      gggCGGTGGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCa                                                                                         <  6:342303/122‑1 (MQ=255)
                                                                        gCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCAtct                                                             <  7:426753/149‑1 (MQ=255)
                                                                            ttGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGcc                                                         >  1:273493/1‑149 (MQ=255)
                                                                                    aTTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCAt                                                 >  3:314234/1‑149 (MQ=255)
                                                                                       cGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGa                                               <  3:292213/148‑1 (MQ=255)
                                                                                        gCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGAtt                                             <  2:284867/149‑1 (MQ=255)
                                                                                             ggCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtctc                                        <  7:444432/149‑1 (MQ=255)
                                                                                             ggCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtctc                                        <  8:12020/149‑1 (MQ=255)
                                                                                             ggCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtctc                                        <  7:449145/149‑1 (MQ=255)
                                                                                              gcctgatgaagagagatgaaaaacaaACCGCGACGCCAGGCGGCATCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATaaa                                                                           <  8:202706‑M1/87‑1 (MQ=255)
                                                                                              gcctgatgaagagagatgaaaaacaaACCGCGACGCCAGGCGGCATCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATaaa                                                                           >  7:202706‑M1/27‑113 (MQ=255)
                                                                                              gCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTcc                                       <  5:260526/149‑1 (MQ=255)
                                                                                                gACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCaa                                     <  6:225570/149‑1 (MQ=255)
                                                                                                gACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCaa                                     <  7:255574/149‑1 (MQ=255)
                                                                                                 tgatgaagagagatgaaaaacaaACCGCGACGCCAGGCGGCATCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtct                                         <  4:34163‑M1/121‑1 (MQ=255)
                                                                                                 tgatgaagagagatgaaaaacaaACCGCGACGCCAGGCGGCATCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtct                                         >  3:34163‑M1/24‑144 (MQ=255)
                                                                                                     gCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCt                                                                                 <  2:174919/100‑1 (MQ=255)
                                                                                                     gCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCt                                                                                 >  1:174919/1‑100 (MQ=255)
                                                                                                     gCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAatat                                  <  6:5784/147‑1 (MQ=255)
                                                                                                      cTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCAt                                                 >  2:131005/1‑131 (MQ=255)
                                                                                                      cTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCAt                                                 <  1:131005/131‑1 (MQ=255)
                                                                                                       tGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCtt                                                                                                                      >  1:163353/1‑61 (MQ=255)
                                                                                                       tGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCt                                                                                                        <  4:63478/75‑1 (MQ=255)
                                                                                                       tGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCt                                                                                                        >  3:63478/1‑75 (MQ=255)
                                                                                                       tGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATaaa                                                                           >  8:437812/1‑104 (MQ=255)
                                                                                                       tGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATaaa                                                                           <  7:437812/104‑1 (MQ=255)
                                                                                                           aCATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAAAATAAccgtta                          >  4:304674/1‑148 (MQ=255)
                                                                                                            gatgaaaaacaaACCGCGACGCCAGGCGGCATCGCTGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCGTCACTTCATCTGATAAAGCAATCTGGCATCTCGCCTTACACATGATTTTCTCCATTATcaccgttgc                         <  1:338503‑M1/137‑3 (MQ=255)
                                                                                                            cATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGa                                               <  6:106862/127‑1 (MQ=255)
                                                                                                            cATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGa                                               >  5:106862/1‑127 (MQ=255)
                                                                                                            cATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAccgttcc                         <  2:273493/149‑1 (MQ=255)
                                                                                                             atatCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAccgttccg                        >  2:359561/1‑149 (MQ=255)
                                                                                                             atatCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGCCATTCCTCACTTCAGCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAccgttccg                        >  5:120879/1‑149 (MQ=255)
                                                                                                                   ctGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAat                                    >  5:449744/1‑131 (MQ=255)
                                                                                                                   ctGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAat                                    <  6:449744/131‑1 (MQ=255)
                                                                                                                   ctGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTg                  >  3:335608/1‑149 (MQ=255)
                                                                                                                    tGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTg                     <  8:282045/145‑1 (MQ=255)
                                                                                                                    tGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTg                     <  8:275045/145‑1 (MQ=255)
                                                                                                                    tGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTg                     >  7:275045/1‑145 (MQ=255)
                                                                                                                    tGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTg                     >  7:282045/1‑145 (MQ=255)
                                                                                                                    tGAGACCGCGATGCCGCATGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 >  8:472886/1‑149 (MQ=255)
                                                                                                                     gagaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 >  3:271980/1‑148 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCCTGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 <  6:166834/147‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTGCGTTGCTgg                 <  6:189852/147‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 <  3:139237/147‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 <  6:400044/147‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTg                  <  5:225227/146‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGa               <  5:334405/149‑1 (MQ=255)
                                                                                                                      agaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGa               >  2:419514/1‑149 (MQ=255)
                                                                                                                       gaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGa               >  5:443986/1‑148 (MQ=255)
                                                                                                                       gaCCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGa               <  6:443986/148‑1 (MQ=255)
                                                                                                                          CGCGATGCCGCCTGGCATCGCGGTTTGTTTTTTCATCACTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTg            <  2:374223‑M1/149‑1 (MQ=255)
                                                                                                                            CGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTa                                                      <  6:76981‑M1/104‑1 (MQ=255)
                                                                                                                            CGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTg            <  5:455216‑M1/146‑1 (MQ=255)
                                                                                                                             GATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                 <  7:207999‑M1/140‑1 (MQ=255)
                                                                                                                             GATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTc         <  8:327180‑M1/148‑1 (MQ=255)
                                                                                                                              ATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTAc                                                     >  7:295761‑M1/1‑103 (MQ=255)
                                                                                                                              ATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTAc                                                     <  8:295761‑M1/103‑1 (MQ=255)
                                                                                                                                  cgccTGGCGTCGCGGTTTGTTTATCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTGTCCAATATCACCGATCCGTTGCAGGGACTGGGCGATCcg   >  2:229721/1‑149 (MQ=255)
                                                                                                                                   gccGGGCGGCGCGGGTTGTTTTTAATCTCTCTTCAGCATGCTTGTCTGCATGGCATTCCTCACTGCATCACCGAAAGCACCCGCGCATCTCCCCTTACCCCTGATTTTCTCCAAAATCACCGTGACGTTGCTGGGACTGGTCGATAcgg  <  6:310096/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGG  >  NC_000913/4257616‑4257895

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