Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,672 C→G 5.1% W63C (TGG→TGC)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6720CG 5.1% 205.1 / 2.6 78W63C (TGG→TGC) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base G (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

TAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATA  >  NC_000913/361524‑361820
                                                                                                                                                    |                                                                                                                                                    
tAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAcc                                                                                                                                                      >  4:3552/1‑149 (MQ=255)
   cGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAg                                                                                                                                                   <  8:3021/149‑1 (MQ=255)
     gAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGc                                                                                                                                                 >  1:96281/1‑149 (MQ=255)
            tAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTcc                                                                                                                                          >  7:157106/1‑149 (MQ=255)
              aCGTTGGGTGCAATCAGTATGATATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTCGAACCGTATGAGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTccac                                                                                                                                        >  8:253984/1‑147 (MQ=255)
                 ttGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTAc                                                                                                                                     >  7:181870/1‑149 (MQ=255)
                    ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGt                                                                                                                                  <  7:383303/149‑1 (MQ=255)
                        cAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCa                                                                                                                              <  6:180627/149‑1 (MQ=255)
                         aaTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGGCGTTTTCCCCTACCGTGGCa                                                                                                                              >  6:366734/1‑148 (MQ=255)
                         aaTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCa                                                                                                                              >  1:332719/1‑148 (MQ=255)
                         aaTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCa                                                                                                                              >  7:35107/1‑148 (MQ=255)
                           tCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCg                                                                                                                                   >  8:151170/1‑141 (MQ=255)
                           tCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCg                                                                                                                                   <  7:151170/141‑1 (MQ=255)
                               tACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACAttt                                                                                                                       <  6:359526/149‑1 (MQ=255)
                                 cGTTATCACCTATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGAATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATATACAGGAGGTTCTACCCAGTCGTTTTCACCTACCGTTGCAAACATTTCt                                                                                                                     <  3:285817/149‑1 (MQ=255)
                                 cGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTc                                                                                                                      <  8:35107/148‑1 (MQ=255)
                                 cGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCt                                                                                                                     <  3:3552/149‑1 (MQ=255)
                                   ttATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCtt                                                                                                                    >  8:263201/1‑148 (MQ=255)
                                      tCACCGATTGTTACCGTGTAGTCATCGACAATTGTTAAATTGATATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCATACATTTCTTTAAt                                                                                                                <  4:123441/149‑1 (MQ=255)
                                      tCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACGCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAAt                                                                                                                >  7:47530/1‑149 (MQ=255)
                                      tCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACGCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAAt                                                                                                                >  7:53995/1‑149 (MQ=255)
                                        aCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAAt                                                                                                                >  8:85502/1‑147 (MQ=255)
                                               gTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCttt                                                                                                        <  2:329532/148‑1 (MQ=255)
                                                ttACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTtct                                                                                                      <  3:296462/149‑1 (MQ=255)
                                                   ccGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCtt                                                                                                   <  8:157106/149‑1 (MQ=255)
                                                    cGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCttt                                                                                                  >  6:307128/1‑149 (MQ=255)
                                                          gTCATCGACAATGTTTAAATTGAAATTTGCATAAAAATTGCTGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCaa                                                                                            <  4:225742/149‑1 (MQ=255)
                                                          gTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCaa                                                                                            <  1:2430/149‑1 (MQ=255)
                                                           tCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAAc                                                                                           >  6:347759/1‑149 (MQ=255)
                                                             aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCTGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACt                                                                                          <  7:253984/148‑1 (MQ=255)
                                                             aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCATACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACtt                                                                                         >  7:143284/1‑149 (MQ=255)
                                                             aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACtt                                                                                         >  5:138230/1‑149 (MQ=255)
                                                             aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTAAACtt                                                                                         >  3:216724/1‑149 (MQ=255)
                                                             aTCGACAATCGTTAAATTGAAATATGCATAAAAATTGCGGCCTATATGCATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACt                                                                                          >  4:318273/1‑148 (MQ=255)
                                                              tCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTc                                                                                        >  7:18932/1‑149 (MQ=255)
                                                              tCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTc                                                                                        >  2:80103/1‑149 (MQ=255)
                                                              tCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTc                                                                                        >  5:337144/1‑149 (MQ=255)
                                                                 aCAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTga                                                                                     <  1:80103/149‑1 (MQ=255)
                                                                 aCAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTga                                                                                     <  2:285470/149‑1 (MQ=255)
                                                                  cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTgat                                                                                    <  5:108795/149‑1 (MQ=255)
                                                                     tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCATACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTgatgga                                                                                 >  3:61786/1‑149 (MQ=255)
                                                                           aaTTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                            <  7:85502/148‑1 (MQ=255)
                                                                           aaTTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAg                                                                           <  2:153274/149‑1 (MQ=255)
                                                                               gAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGa                                                                        <  8:18932/148‑1 (MQ=255)
                                                                               gAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGa                                                                        <  5:307128/148‑1 (MQ=255)
                                                                                  aTTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTaa                                                                    <  6:337144/149‑1 (MQ=255)
                                                                                   tttGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTaaa                                                                   <  3:366600/149‑1 (MQ=255)
                                                                                     tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAAc                                                                  <  7:263201/148‑1 (MQ=255)
                                                                                      gCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACt                                                                 <  4:311841/148‑1 (MQ=255)
                                                                                        aTAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCa                                                               >  2:273193/1‑148 (MQ=255)
                                                                                         tAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATa                                                             <  6:261343/149‑1 (MQ=255)
                                                                                          aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGc                                                                              <  4:79333/131‑1 (MQ=255)
                                                                                          aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGc                                                                              >  3:79333/1‑131 (MQ=255)
                                                                                          aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATa                                                             >  7:215792/1‑148 (MQ=255)
                                                                                           aaaaTTGCGGCCTATATGGATGTTGGAACCGTTAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACa                                                           <  6:129084/149‑1 (MQ=255)
                                                                                                     ccTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACGCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGtttt                                                 <  8:47530/149‑1 (MQ=255)
                                                                                                     ccTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACGCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGtttt                                                 <  8:53995/149‑1 (MQ=255)
                                                                                                     ccTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGttt                                                  >  4:81667/1‑148 (MQ=255)
                                                                                                        atatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTcc                                               >  6:96937/1‑148 (MQ=255)
                                                                                                         tatGGAAGATGGAACCGTAAGAGAAATAGACAGCCGGTAATACCCAGCCGTTTACACCTACCGTGGCAAACATATCTTTAATCAGGGTGTCACTGTTTGCAAGTTCGGCTGGATGCGAGTGATTAAACTCATACATTAACGTGTTccag                                             >  4:46637/1‑147 (MQ=255)
                                                                                                          atGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACA‑TAACGTTTTCCCAc                                            <  2:332719/148‑1 (MQ=255)
                                                                                                               tGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGt                                       >  7:131692/1‑149 (MQ=255)
                                                                                                                gTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCa                                             >  7:324433/1‑142 (MQ=255)
                                                                                                                gTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCa                                             <  8:324433/142‑1 (MQ=255)
                                                                                                                   ggAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCCAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCttt                                   >  6:15437/1‑149 (MQ=255)
                                                                                                                   ggAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTACAAATTGTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCttt                                   >  4:320296/1‑149 (MQ=255)
                                                                                                                       ccGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTc                               <  5:366734/149‑1 (MQ=255)
                                                                                                                       ccGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTc                               >  5:376856/1‑149 (MQ=255)
                                                                                                                         gTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCg                             >  4:341650/1‑149 (MQ=255)
                                                                                                                         gTAAGAGAAATAGACAGGCGGTTCTACACACGCGTTTTCCCCTACCGTGGCAAAAATTTCTTTAATCAGTCTTTCTCTTATTTCAACTTCTGATCGAAGCGAGTGATTAAACTCATACATTAACCTTTTCACACCAAGTCTTTTTTACg                             >  5:30002/1‑149 (MQ=255)
                                                                                                                          tAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCgg                            >  7:89730/1‑149 (MQ=255)
                                                                                                                                    aGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                            >  6:187757/1‑91 (MQ=255)
                                                                                                                                    aGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                            <  5:187757/91‑1 (MQ=255)
                                                                                                                                     gACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCg                  <  3:81667/148‑1 (MQ=255)
                                                                                                                                      aCAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                <  3:318273/149‑1 (MQ=255)
                                                                                                                                      aCAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                >  1:79539/1‑149 (MQ=255)
                                                                                                                                          gCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATTGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACat             >  8:378932/1‑148 (MQ=255)
                                                                                                                                             gTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCg         >  1:33503/1‑149 (MQ=255)
                                                                                                                                                 tACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTa      >  2:55186/1‑148 (MQ=255)
                                                                                                                                                   cccAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGTCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTa      <  4:216724/146‑1 (MQ=255)
                                                                                                                                                   cccAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAAt   <  8:215792/149‑1 (MQ=255)
                                                                                                                                                    ccAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTaaa    >  8:241909/1‑147 (MQ=255)
                                                                                                                                                    ccAGGCGTTTTCCCCTACCGAGGCACACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATa  <  8:181870/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATA  >  NC_000913/361524‑361820

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