Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,836 G→A 5.5% pseudogene (33/198 nt) ypaA → pseudogene, C‑terminal fragment, Transposase_31 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,8360GA 5.5% 302.1 / 2.2 95pseudogene (33/198 nt)ypaApseudogene, C‑terminal fragment, Transposase_31 family
Reads supporting (aligned to +/- strand):  ref base G (52/35);  new base A (3/2);  total (56/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAG  >  NC_000913/2357697‑2357975
                                                                                                                                           |                                                                                                                                           
aaCAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagg                                                                                                                                      <  2:143478/147‑1 (MQ=255)
   aaGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCaggaagg                                                                                                                                  <  1:214737/148‑1 (MQ=255)
          tAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACATCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCa                                                                                                                            >  8:89883/1‑147 (MQ=255)
          tAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAGTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCa                                                                                                                            >  6:419777/1‑147 (MQ=255)
            aaaTTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaa                                                                                                                           <  8:318895/146‑1 (MQ=255)
            aaaTTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaa                                                                                                                           >  7:318895/1‑146 (MQ=255)
                tGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggc                                                                                                                     >  8:16370/1‑148 (MQ=255)
                           tGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGAAAACCCAACTCACTGATGATGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                          >  8:359330/1‑148 (MQ=255)
                               ggAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaa                                                                                                      >  2:463491/1‑148 (MQ=255)
                               ggAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTAACTGATGCTGATATTCCGCACTGTCATTAAAAGGAGAGATTAATGAGTATTGCTGAACGGCTGCGGCAGGAAGTACATCAAATTGGCTGGGAGGGAGGTaag                                                                                                      >  4:267228/1‑147 (MQ=255)
                                 aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtgg                                                                                                                      >  1:82510/1‑130 (MQ=255)
                                 aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtgg                                                                                                                      <  2:82510/130‑1 (MQ=255)
                                    aGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGa                                                                                                 <  1:384961/148‑1 (MQ=255)
                                       gCTTTGAACGTGACATTGTGCTGGCGACAACCCAACTCACTGATGCAGATATTGCAAACTGTCAGTAAAAGTAAACATTATTGAATTTTGGTGGAGGACGCAGGCAGTTAGAACATAAAATTGGATGGCCGTAAGTTAAAtatg                                                                                                  >  2:427990/1‑141 (MQ=255)
                                           tGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                            <  5:419777/146‑1 (MQ=255)
                                            gAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                          <  8:94605/147‑1 (MQ=255)
                                            gAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                          >  7:94605/1‑147 (MQ=255)
                                             aaCGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                         <  7:281424/147‑1 (MQ=255)
                                               cGTGAGATTGTGCTGGCGACAACCCAACTCACTGAGGCTGATATTCCGCACTGTCATTATAAGGTGAGATTAATGACTATTGCTGAACGGCTGCGGGAGGAAGGCCATCAAATTGGCTGGCAGGAAGGTAAAGTAGAAGGTTTGCag                                                                                       >  1:121941/1‑146 (MQ=25)
                                                gTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATg                                                                                      >  4:99020/1‑147 (MQ=255)
                                                 tGAGATTGTGCTGGCGACCACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                    >  4:225444/1‑148 (MQ=255)
                                                 tGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                    >  6:55619/1‑148 (MQ=255)
                                                 tGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                    <  5:55619/148‑1 (MQ=255)
                                                 tGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                    >  1:402501/1‑148 (MQ=255)
                                                     aTTGTGCTGGCGACCACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCATATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAg                                                                                >  1:105982/1‑148 (MQ=255)
                                                     aTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaa                                                                                                      >  4:79023/1‑126 (MQ=255)
                                                     aTTGTGCTGGCGACAACCCAACTCACTGATGCGGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaa                                                                                                      <  3:79023/126‑1 (MQ=255)
                                                              gCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTaaaa                                                                       <  7:16370/148‑1 (MQ=255)
                                                               cGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAt                                                                      >  2:231502/1‑148 (MQ=255)
                                                                gACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAtt                                                                     >  4:159012/1‑148 (MQ=255)
                                                                 aCAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGTTTAATGACTATTGCTGAACGGCTGCTGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                         <  8:467030/127‑1 (MQ=255)
                                                                 aCAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                         >  7:467030/1‑127 (MQ=255)
                                                                  cAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGc                                                                   >  1:287997/1‑148 (MQ=255)
                                                                  cAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGc                                                                   >  1:294752/1‑148 (MQ=255)
                                                                   atcccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCt                                                                  <  7:89883/146‑1 (MQ=255)
                                                                    aCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCtt                                                                 <  1:463491/148‑1 (MQ=255)
                                                                     cccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGa                                                                                     >  2:257263/1‑127 (MQ=255)
                                                                     cccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGa                                                                                     <  1:257263/127‑1 (MQ=255)
                                                                     cccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCttt                                                                >  1:249322/1‑148 (MQ=255)
                                                                       cAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgc                                                              >  3:237561/1‑148 (MQ=255)
                                                                         aCTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTCATGACTATTGCTGAAAGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                            <  2:287997/148‑1 (MQ=255)
                                                                         aCTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                             >  1:68241/1‑147 (MQ=255)
                                                                         aCTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                             >  8:247622/1‑147 (MQ=255)
                                                                          cTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTACTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCa                                                           >  6:267282/1‑148 (MQ=255)
                                                                            cACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAATGCTTTGCGCATg                                                         >  4:183036/1‑148 (MQ=255)
                                                                             aCTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                            <  2:294752/144‑1 (MQ=255)
                                                                             aCTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCAt                                                          <  3:225444/146‑1 (MQ=255)
                                                                             aCTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCAt                                                          <  2:402501/146‑1 (MQ=255)
                                                                              cTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCt                                                       >  7:436479/1‑148 (MQ=255)
                                                                                 aTGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGa                                                    <  3:159012/148‑1 (MQ=255)
                                                                                   gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                       <  3:155183/95‑1 (MQ=255)
                                                                                   gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                       >  4:155183/1‑95 (MQ=255)
                                                                                      gATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg                                               >  4:118120/1‑148 (MQ=255)
                                                                                           tCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTaaaa                                                                       <  2:347815/119‑1 (MQ=255)
                                                                                           tCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTaaaa                                                                       >  1:347815/1‑119 (MQ=255)
                                                                                              gAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAGCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                       >  3:179716/1‑148 (MQ=255)
                                                                                              gAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                       <  7:359330/148‑1 (MQ=255)
                                                                                              gAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                       <  3:118120/148‑1 (MQ=255)
                                                                                                     cATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGaa                                                                                                            >  7:468323/1‑72 (MQ=255)
                                                                                                     cATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGaa                                                                                                            <  8:468323/72‑1 (MQ=255)
                                                                                                     cATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAc                                 >  3:133747/1‑147 (MQ=255)
                                                                                                      aTTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                               <  1:231502/148‑1 (MQ=255)
                                                                                                       ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                               >  8:106666/1‑147 (MQ=255)
                                                                                                       ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                               <  7:106666/147‑1 (MQ=255)
                                                                                                       ttAAAAGAAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAg                              >  4:127726/1‑148 (MQ=255)
                                                                                                        tAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAgg                             >  3:145206/1‑148 (MQ=255)
                                                                                                        tAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAgg                             >  5:427590/1‑148 (MQ=255)
                                                                                                          aaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                       <  1:23278/72‑1 (MQ=40)
                                                                                                          aaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTg                           >  2:62328/1‑148 (MQ=255)
                                                                                                          aaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                       >  2:23278/1‑72 (MQ=14)
                                                                                                               agagaTTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGt                                                                                             <  4:184446/77‑1 (MQ=255)
                                                                                                               agagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAt                                                                                                                              >  7:419250/1‑44 (MQ=255)
                                                                                                               agagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAt                                                                                                                              <  8:419250/44‑1 (MQ=255)
                                                                                                               agagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGt                                                                                             >  3:184446/1‑77 (MQ=255)
                                                                                                                    ttAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                  <  6:218113/115‑1 (MQ=255)
                                                                                                                    ttAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                  >  5:218113/1‑115 (MQ=255)
                                                                                                                      aaTGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAATTTGCGTTGCGCATGCTGGATCAGGGCTTTGATCGTTACCCGGTGCTCGCGGCCAcc               >  2:363703/1‑148 (MQ=37)
                                                                                                                        tGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAcc               >  7:45448/1‑146 (MQ=255)
                                                                                                                        tGACTATTGCAGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAcc               >  7:262396/1‑146 (MQ=255)
                                                                                                                         gACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAg            >  7:319169/1‑148 (MQ=255)
                                                                                                                             agtGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGATTGCATGAACAAGCCATTAAACTTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa        <  3:267228/146‑1 (MQ=255)
                                                                                                                             aTTGCTGCACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa        >  6:10192/1‑148 (MQ=255)
                                                                                                                             aTTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCt          >  1:34942/1‑146 (MQ=255)
                                                                                                                             aTTGCTGAACGGCTCCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                            <  5:179064/64‑1 (MQ=255)
                                                                                                                             aTTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                            >  6:179064/1‑64 (MQ=255)
                                                                                                                             aTTGCTGAACGGCTACGGCAGGAAGGCCATCAAATTGGCTGGCAGGAAGGTAAAtt                                                                                                    <  1:253329/56‑1 (MQ=11)
                                                                                                                             aTTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAAtt                                                                                                    <  1:246280/56‑1 (MQ=255)
                                                                                                                               tGCTGAACGGCTGCGGTAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGc      >  6:266738/1‑148 (MQ=255)
                                                                                                                                   gAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                  >  1:88359/1‑100 (MQ=255)
                                                                                                                                   gAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                  <  2:88359/100‑1 (MQ=255)
                                                                                                                                   gAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg  >  4:112992/1‑148 (MQ=255)
                                                                                                                                    aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg  <  5:267282/147‑1 (MQ=255)
                                                                                                                                    aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg  <  2:68241/147‑1 (MQ=255)
                                                                                                                                     aCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg                                               >  2:215679/1‑101 (MQ=255)
                                                                                                                                     aCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg                                               <  1:215679/101‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAG  >  NC_000913/2357697‑2357975

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: