Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,988,051 A→T 7.9% intergenic (‑92/‑171) azuC ← / → yecR acid‑inducible small membrane‑associated protein/lipoprotein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,0510AT 7.9% 217.7 / 7.2 75intergenic (‑92/‑171)azuC/yecRacid‑inducible small membrane‑associated protein/lipoprotein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (3/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTT  >  NC_000913/1987908‑1988199
                                                                                                                                               |                                                                                                                                                    
aCGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATca                                                                                                                                                  >  3:136433/1‑148 (MQ=255)
  gTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACAt                                                                                                                                               <  7:15402/149‑1 (MQ=255)
     ttGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATaaa                                                                                                                                            >  4:150022/1‑149 (MQ=255)
       gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATa                                                                                                                                          >  5:306970/1‑149 (MQ=255)
        cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCTCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaaga                                                                                                                                          <  2:3399/148‑3 (MQ=255)
        cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaaga                                                                                                                                          <  5:187172/148‑3 (MQ=255)
             gTTAATCAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACGCCGACCATTTGTTATACAGGTGTAAATGTCTTTGCGCTTATAATCTGCCAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCtt                                                                                                                                    <  2:108131/149‑1 (MQ=255)
                   gAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATg                                                                                                                               <  6:200767/148‑1 (MQ=255)
                   gAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATg                                                                                                                               <  3:352109/148‑1 (MQ=255)
                    aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                             >  2:59953/1‑149 (MQ=255)
                    aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGt                                                                                                                              <  8:251719/148‑1 (MQ=255)
                     aCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                             <  6:306970/148‑1 (MQ=255)
                     aCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTta                                                                                                                            <  1:87746/149‑1 (MQ=255)
                     aCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTta                                                                                                                            >  3:179030/1‑149 (MQ=255)
                      cATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTta                                                                                                                            >  1:7745/1‑148 (MQ=255)
                               cAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTACAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGa                                                                                                                   >  1:76636/1‑148 (MQ=255)
                                aGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATa                                                                                                                 <  7:256773/149‑1 (MQ=255)
                                   aTTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATc                                                                                                              <  4:207423/149‑1 (MQ=255)
                                   aTTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATc                                                                                                              <  8:373105/149‑1 (MQ=255)
                                         cgcAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTTCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCaaaaaa                                                                                                        <  5:331261/149‑1 (MQ=255)
                                         cgcAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCaaaaaa                                                                                                        <  2:246531/149‑1 (MQ=255)
                                           cAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCaaaaaa                                                                                                        >  5:330117/1‑147 (MQ=255)
                                            aGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAg                                                                                                     <  8:293322/149‑1 (MQ=255)
                                            aGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAg                                                                                                     <  4:136433/149‑1 (MQ=255)
                                                  aTCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACCCCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGcccccc                                                                                               >  5:96174/1‑149 (MQ=255)
                                                   tCGCGGCACTCCGACCATTTGTTATACAGGTTTTATTGTCTTTGCGCTTATAATATTAGAACCATGTCAATATTCATTCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAATAATGATAATCAAAAAACAGCCCCCCt                                                                                              <  7:237287/149‑1 (MQ=255)
                                                     gcgGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                             >  3:130709/1‑116 (MQ=255)
                                                     gcgGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                             <  4:130709/116‑1 (MQ=255)
                                                         cACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCt                                                                                        <  7:14461/149‑1 (MQ=255)
                                                          aCTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCt                                                                                        >  1:32262/1‑148 (MQ=255)
                                                           cTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                      <  3:254305/149‑1 (MQ=255)
                                                            tCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                      >  6:167329/1‑148 (MQ=255)
                                                            tCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                      >  1:388814/1‑148 (MQ=255)
                                                              cGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGacac                                                                                   >  6:64607/1‑149 (MQ=255)
                                                              cGACCAGTTGTTATACAGGTGTTATTGTCTTTGAGCTTCTAATATGACAACCATAACAACAATCAATCTTTATGTGATCCACATCACATAAATACCCCTTTAATGTTATAAAAATGCTCATCAAACAACAGCCCCCCTATTTCTGaca                                                                                    <  4:76356/148‑1 (MQ=255)
                                                              agcccATTTGTTCTACATCTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAACTCACTCTTTATGTGATACACATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGaca                                                                                    <  6:276261/145‑1 (MQ=255)
                                                                  cATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTAc                                                                               <  3:326738/149‑1 (MQ=255)
                                                                    tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGc                                                                                                    >  7:11553/1‑126 (MQ=255)
                                                                    tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGc                                                                                                    <  8:11553/126‑1 (MQ=255)
                                                                        ttATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCACCCCTATTTCTGACACGTACAGATgg                                                                         >  6:255288/1‑149 (MQ=255)
                                                                            aCAGGTGTTATTGTCTTTGCGCTTATAATTTGGCATCCATCAAAAAAATCAATCTTTATGTGATACAAATCACATAAATTCACCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCGTATTTCGGACTCCGACCGATggg                                                                        >  2:308229/1‑145 (MQ=255)
                                                                            aCAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGc                                                                        >  3:254870/1‑146 (MQ=255)
                                                                                   ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAgc                                                              >  5:92232/1‑149 (MQ=255)
                                                                                   ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAgc                                                              <  4:179030/149‑1 (MQ=255)
                                                                                      ttGTCTTTGCGCTTATAATATGACAACCATCACAATATTCAATCTTTATGTGATACATCTGACATAAATACCCCTTTAATGTTATAAAAATGATAATGAAAAAAAAGCCCCCCTATTTCTGAACCCTACAGATCGCAAGAAATAGCGcc                                                           >  1:210466/1‑149 (MQ=255)
                                                                                         tCTTTGCGCTTATAATATGACAACCATCAAAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATTATAATCAAAAAACAGCCCCCCGATTTCTGACACCTACAGATGGCAAGAAATAGCGCCAGc                                                        >  8:276361/1‑149 (MQ=255)
                                                                                                     caataTGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                             <  2:76636/147‑1 (MQ=255)
                                                                                                       atatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTcc                                          >  2:111692/1‑149 (MQ=255)
                                                                                                        tatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACCGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCg                                         >  4:19400/1‑149 (MQ=255)
                                                                                                        tatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCg                                         >  5:9173/1‑149 (MQ=255)
                                                                                                        tatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCg                                         <  5:156507/149‑1 (MQ=255)
                                                                                                             cAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAt                                    >  8:356580/1‑149 (MQ=255)
                                                                                                                ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                             >  1:235619/1‑137 (MQ=255)
                                                                                                                ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCttt                                             <  2:235619/137‑1 (MQ=255)
                                                                                                                ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGt                                 >  2:58548/1‑149 (MQ=255)
                                                                                                                 cATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAAAGTTATAAAAATGATAATCAAAAAACAGCCCCAATATTTGTGACACCTACAGATGGCAAGATATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTc                                >  8:357615/1‑149 (MQ=255)
                                                                                                                  aTCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTTTCg                               >  5:220162/1‑149 (MQ=255)
                                                                                                                     acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgcag                            <  2:7745/149‑1 (MQ=255)
                                                                                                                     acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                             <  6:94265/148‑1 (MQ=255)
                                                                                                                     acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                             <  6:92232/148‑1 (MQ=255)
                                                                                                                     acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                             <  6:330117/148‑1 (MQ=255)
                                                                                                                     acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATATCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgca                             <  2:80164/148‑1 (MQ=255)
                                                                                                                       aaaaaTCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgcagca                          <  6:96174/149‑1 (MQ=255)
                                                                                                                        aaaaTCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCAc                         <  2:205648/149‑1 (MQ=255)
                                                                                                                             caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCACCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTTCCAGGCGTATTTTCCGGCCATTGTCGCAGAACTGTaa                    >  6:297751/1‑149 (MQ=255)
                                                                                                                              ggtatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAt                                    <  7:232711/128‑1 (MQ=255)
                                                                                                                              ggtatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAt                                    >  8:232711/5‑132 (MQ=255)
                                                                                                                                tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCaa                                                                      >  7:122917/1‑96 (MQ=255)
                                                                                                                                tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCaa                                                                      <  8:122917/96‑1 (MQ=255)
                                                                                                                                tatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAcgc                 >  3:29204/3‑149 (MQ=255)
                                                                                                                                  tttATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAcgcg                >  1:188081/1‑148 (MQ=255)
                                                                                                                                        tgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATa         <  1:59953/149‑1 (MQ=255)
                                                                                                                                        tgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATa         <  4:50075/149‑1 (MQ=255)
                                                                                                                                        tgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATa         <  5:64607/149‑1 (MQ=255)
                                                                                                                                           tACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTg      >  8:310295/1‑149 (MQ=255)
                                                                                                                                               aaTCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCttt  <  5:167329/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTTT  >  NC_000913/1987908‑1988199

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