Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,183,048 T→G L191F (TTA→TTC potC ← spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,183,0480TG96.9% 415.8 / 0.2 128L191F (TTA→TTCpotCspermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base G (63/61);  total (63/65)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

CTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTCGCCGCTTCC  >  NC_000913/1182951‑1183144
                                                                                                 |                                                                                                
cTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGaa                                                                                                <  1:1388564/100‑1 (MQ=255)
 tttGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAac                                                                                               <  1:833043/100‑1 (MQ=255)
 tttGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAac                                                                                               <  2:2967534/100‑1 (MQ=255)
 tttGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAac                                                                                               <  2:2107374/100‑1 (MQ=255)
  ttGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAaca                                                                                              >  2:1975639/1‑100 (MQ=255)
  ttGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAaca                                                                                              >  2:3783149/1‑100 (MQ=255)
   tGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAacac                                                                                             <  1:1069041/100‑1 (MQ=255)
        aTCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCCTCGACAGGGTAAAGCTGAACACCCAGc                                                                                        >  1:3540179/1‑100 (MQ=255)
        aTCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGc                                                                                        <  1:2487466/100‑1 (MQ=255)
        aTCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGc                                                                                        >  1:2160933/1‑100 (MQ=255)
         tCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGcc                                                                                       <  1:1248628/100‑1 (MQ=255)
          cGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGccc                                                                                      >  1:2355587/1‑100 (MQ=255)
           gAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCccg                                                                                     <  1:1740192/100‑1 (MQ=255)
             aTAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCccgcc                                                                                   >  2:3672566/1‑100 (MQ=255)
              tAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAATCTTAACACCCAGCccgccg                                                                                  <  2:3095363/100‑1 (MQ=255)
              tAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCccgccg                                                                                  >  1:3457481/1‑100 (MQ=255)
                aaTTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCgccgcc                                                                                <  2:2630504/100‑1 (MQ=255)
                aaTTTTTAACGGCAGAATTGCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCgccgcc                                                                                <  2:859376/100‑1 (MQ=255)
                  tttttAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCAc                                                                              <  2:370591/100‑1 (MQ=255)
                   ttttAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCAcc                                                                             <  2:927954/100‑1 (MQ=255)
                     ttAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCCCACCGTACATCGACAGGGTAAAGCGGAACACGCAGCCCGCCGCCACCGc                                                                           >  2:3358078/1‑100 (MQ=255)
                        aCGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTgg                                                                        >  1:3870459/1‑100 (MQ=255)
                         cGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAACACCCAGCTCGCCTCCACCGCTGGc                                                                       <  2:3876225/100‑1 (MQ=255)
                         cGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGc                                                                       >  1:615969/1‑100 (MQ=255)
                         cGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCGGGc                                                                       >  2:2303783/1‑100 (MQ=255)
                         cGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCCTCCATCAACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGc                                                                       >  2:2763975/1‑100 (MQ=255)
                          ggCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCa                                                                      >  2:2398980/1‑100 (MQ=255)
                          ggCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCa                                                                      <  1:1447208/100‑1 (MQ=255)
                           gCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCAt                                                                     <  2:3423317/100‑1 (MQ=255)
                           gCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCAt                                                                     >  1:510791/1‑100 (MQ=255)
                           gCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCAt                                                                     >  2:2328144/1‑100 (MQ=255)
                            cAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCAtt                                                                    >  2:261652/1‑100 (MQ=255)
                             aGAATTGCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTg                                                                   <  2:3454379/100‑1 (MQ=255)
                              gAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGc                                                                  >  2:2367824/1‑100 (MQ=255)
                               aaTTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGcc                                                                 <  2:1797043/100‑1 (MQ=255)
                                aTTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCa                                                                >  1:2697062/1‑100 (MQ=255)
                                aTTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCa                                                                <  2:1261226/100‑1 (MQ=255)
                                 tttCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAg                                                               >  1:3019998/1‑100 (MQ=255)
                                    cATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTgg                                                            >  1:579353/1‑100 (MQ=255)
                                     ntAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGc                                                           >  1:225533/2‑100 (MQ=255)
                                     aTAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGc                                                           <  2:1941080/100‑1 (MQ=255)
                                     aTAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGc                                                           <  1:2492262/100‑1 (MQ=255)
                                      tAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCa                                                          >  2:3109822/1‑100 (MQ=255)
                                      tAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCa                                                          <  1:849276/100‑1 (MQ=255)
                                      tAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCa                                                          <  1:3157618/100‑1 (MQ=255)
                                       aaCTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAg                                                         <  1:2765966/100‑1 (MQ=255)
                                       aaCTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGGGGCAg                                                         >  2:2118668/1‑100 (MQ=255)
                                          ttGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAAt                                                      <  1:1711692/100‑1 (MQ=255)
                                           tGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATg                                                     >  1:2178806/1‑100 (MQ=255)
                                            ggTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACCGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGa                                                    >  2:1297175/1‑100 (MQ=255)
                                            ggTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGa                                                    <  2:1988723/100‑1 (MQ=255)
                                             gTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGAt                                                   >  1:4033217/1‑100 (MQ=255)
                                             gTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGAt                                                   >  2:357796/1‑100 (MQ=255)
                                             gTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGAt                                                   <  1:4059997/100‑1 (MQ=255)
                                              tCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGAtt                                                  <  2:481566/100‑1 (MQ=255)
                                               ccGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGAttt                                                 >  2:482729/1‑100 (MQ=255)
                                                  gTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTc                                              >  1:1786404/1‑100 (MQ=255)
                                                   tGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCg                                             <  2:645979/100‑1 (MQ=255)
                                                   tGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCg                                             <  1:2444733/100‑1 (MQ=255)
                                                      cgaacgaaGACACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAg                                          <  1:233375/100‑1 (MQ=255)
                                                      cgaacgaaGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAg                                          >  2:2296374/1‑100 (MQ=255)
                                                        aacgaaGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGaa                                        <  2:3546286/100‑1 (MQ=255)
                                                         acgaaGAAACGACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCGGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAAt                                       <  1:2799971/100‑1 (MQ=255)
                                                          cgaaGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATg                                      <  2:2319362/100‑1 (MQ=255)
                                                          cgaaGAAACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATg                                      >  2:2654946/1‑100 (MQ=255)
                                                           gaagaaACTACCACATCGTCCATCGACAGGGTAAAGCTGAACCCCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATgg                                     <  2:567630/100‑1 (MQ=255)
                                                             agaaACTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTa                                   <  2:3485404/100‑1 (MQ=255)
                                                                aaCTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAAt                                <  2:2017159/100‑1 (MQ=255)
                                                                aaCTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAAt                                <  2:3081980/100‑1 (MQ=255)
                                                                aaCTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAAt                                >  2:2836664/1‑100 (MQ=255)
                                                                  cTACCACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTc                              <  2:2065595/100‑1 (MQ=255)
                                                                     ccACATCGTCCATCGACAGGGTAAAGCTTAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCt                           <  2:3387720/100‑1 (MQ=255)
                                                                     ccACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCt                           <  2:1798155/100‑1 (MQ=255)
                                                                     ccACATCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCt                           >  1:3695842/1‑100 (MQ=255)
                                                                        caTCGTCCATCGACAGGGTAAAGCTGAACACGCAGCCCGCCGCCACCGCTGGCATTGCCGGTGGCAGAATGAGTTTACGCAGAATGGTAAATTCGCTGGc                        <  2:1799672/100‑1 (MQ=255)
                                                                        caTCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGc                        >  2:1107513/1‑100 (MQ=255)
                                                                        caTCGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGc                        <  2:2288636/100‑1 (MQ=255)
                                                                           cGTCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCAcc                     <  2:3622394/100‑1 (MQ=255)
                                                                             tCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCga                   >  1:1147440/1‑100 (MQ=255)
                                                                             tCCATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCga                   <  2:3768145/100‑1 (MQ=255)
                                                                               cATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCgaga                 >  1:1626464/1‑100 (MQ=255)
                                                                               cATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCgaga                 >  1:2118786/1‑100 (MQ=255)
                                                                               cATCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGAGGCAGAATGATTTTTCCCAGAATGGTAAATTCGCTGGCACAgaga                 >  1:1626685/1‑100 (MQ=255)
                                                                                aTCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGAt                <  1:1094767/100‑1 (MQ=255)
                                                                                 tCGACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATc               <  1:2322545/100‑1 (MQ=255)
                                                                                  ngACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCACAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCt              >  2:477690/2‑100 (MQ=255)
                                                                                  cGACAGTGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGAGCt              >  2:1195693/1‑100 (MQ=255)
                                                                                   gACAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCtt             >  2:2999889/1‑100 (MQ=255)
                                                                                     cAGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTc           <  2:1176942/100‑1 (MQ=255)
                                                                                      agtgTAAAGCTGAACACCCAGCCCGCCTCCACCGCTGGCATTGCCAGTGGCATAATGATTTTTCGCTGAATGGTAAATTCGCTGGCACCGAGATCTTTcg          <  2:3320584/97‑1 (MQ=255)
                                                                                      aGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcg          <  1:2880006/100‑1 (MQ=255)
                                                                                      aGGGTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcg          <  1:2331981/100‑1 (MQ=255)
                                                                                       gggTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcgc         >  1:3285083/1‑100 (MQ=255)
                                                                                        ggTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcgcc        <  2:622728/100‑1 (MQ=255)
                                                                                         gTAAAGCTGCACACCCAGCCCGCCGCCACCGCTGGCATTGCAAGTGGCAGAATGATTTTTCGCAGATTAGTAAATTCGCTAGCACCGAGATCTTTcgccg       <  2:221895/100‑1 (MQ=255)
                                                                                         gTAAAGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcgccg       >  2:3782852/1‑100 (MQ=255)
                                                                                          tAAAGCTGAACATCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATGTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTcgccgc      >  2:543898/1‑100 (MQ=255)
                                                                                             aGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCTTTTTTTTTTAAATTCGCTGGCACCGAGATCTTTCGCCGCCTc   >  2:901525/1‑100 (MQ=255)
                                                                                             aGCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTCGCCGCTTc   >  1:2948536/1‑100 (MQ=255)
                                                                                              gCTGAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTCGCCGCTTcc  >  2:1016919/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAACACCCAGCCCGCCGCCACCGCTGGCATTGCCAGTGGCAGAATGATTTTTCGCAGAATGGTAAATTCGCTGGCACCGAGATCTTTCGCCGCTTCC  >  NC_000913/1182951‑1183144

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