Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,211,469 G→A 5.5% S82N (AGC→AAC)  yehR → lipoprotein, DUF1307 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,211,4690GA 5.5% 236.7 / 2.4 72S82N (AGC→AAC) yehRlipoprotein, DUF1307 family
Reads supporting (aligned to +/- strand):  ref base G (28/40);  new base A (2/2);  total (30/42)
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.60e-01

CCAATCTGAACGGCACTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAACGTTTCTGCT  >  NC_000913/2211325‑2211608
                                                                                                                                                |                                                                                                                                           
ccAATCTGAACGGCACTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGcc                                                                                                                                         >  1:163811/1‑149 (MQ=255)
 cAATCTGAACGGCACTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAACGCCa                                                                                                                                        >  1:106903/1‑149 (MQ=255)
 cAATCTGAACGGCACTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAACGCCa                                                                                                                                        >  1:313709/1‑149 (MQ=255)
              aCTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAATCAAAATTCAATTTGCCTCCATTGGTACAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACAt                                                                                                                           <  4:255809/149‑1 (MQ=255)
               cTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATc                                                                                                                          >  3:450221/1‑149 (MQ=255)
                tGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCACCCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTCAGCGCCAAATACAAAAACATcg                                                                                                                         >  4:90826/1‑149 (MQ=255)
                tGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATcg                                                                                                                         >  2:565394/1‑149 (MQ=255)
                   aaTTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCgg                                                                                                                      <  3:396637/149‑1 (MQ=255)
                     ttGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGgt                                                                                                                    >  2:599444/1‑149 (MQ=255)
                     ttGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGgt                                                                                                                    >  4:679922/1‑149 (MQ=255)
                     ttGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGgt                                                                                                                    <  4:324648/149‑1 (MQ=255)
                     ttGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGgt                                                                                                                    <  1:43374/149‑1 (MQ=255)
                        ccATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGtt                                                                                                                 >  2:143393/1‑149 (MQ=255)
                         cATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTg                                                                                                                >  3:9020/1‑149 (MQ=255)
                          aTTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTga                                                                                                               >  4:121600/1‑149 (MQ=255)
                           ttACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTgaa                                                                                                              <  1:143393/149‑1 (MQ=255)
                           ttACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTgaa                                                                                                              <  2:201558/149‑1 (MQ=255)
                           ttACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTgaa                                                                                                              <  4:74035/149‑1 (MQ=255)
                           ttACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTACCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTgaa                                                                                                              <  1:455636/149‑1 (MQ=255)
                                 aTGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCAAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGAGTTGAAGaaaaa                                                                                                        >  4:148496/1‑149 (MQ=255)
                                    tCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTa                                                                                                     <  1:581860/149‑1 (MQ=255)
                                      tACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  3:503438/149‑1 (MQ=255)
                                      tACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  3:159488/149‑1 (MQ=255)
                                      tACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  1:635096/149‑1 (MQ=255)
                                      tACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   >  4:282407/1‑149 (MQ=255)
                                      tACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  1:221292/149‑1 (MQ=255)
                                         aaaGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCta                                                                                                <  1:318999/149‑1 (MQ=255)
                                         aaaGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCta                                                                                                <  2:414462/149‑1 (MQ=255)
                                            ggTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATAc                                                                                             <  1:408671/149‑1 (MQ=255)
                                               gACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCACAGAATATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGa                                                                                          <  3:644729/149‑1 (MQ=255)
                                                  aaGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATAc                                                                                       <  3:17663/149‑1 (MQ=255)
                                                     gTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTa                                                                                    <  1:612904/149‑1 (MQ=255)
                                                           aaGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCa                                                                              >  2:122461/1‑149 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  2:533024/126‑1 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   >  1:533024/1‑126 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   >  3:292081/1‑126 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAAc                                                                                                   <  4:292081/126‑1 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAgg                                                                            <  1:508448/149‑1 (MQ=255)
                                                             gCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAgg                                                                            >  4:201047/1‑149 (MQ=255)
                                                                 tcttctGAAACCAAAACTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCCACGCGCAGGaaaa                                                                        <  2:338151/149‑1 (MQ=255)
                                                                  cttctGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAAc                                                                       <  4:99699/149‑1 (MQ=255)
                                                                   ttctGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACg                                                                      >  4:229366/1‑149 (MQ=255)
                                                                      tGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGa                                                                   <  4:411502/149‑1 (MQ=255)
                                                                       gAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGAc                                                                  >  1:645807/1‑149 (MQ=255)
                                                                            cAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAGATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCg                                                             >  3:434907/1‑149 (MQ=255)
                                                                             aaaaTTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGa                                                            >  3:407952/1‑149 (MQ=255)
                                                                              aaaTTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGTCTATCGat                                                           >  2:573230/1‑149 (MQ=255)
                                                                                aTTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGatat                                                         <  1:558560/149‑1 (MQ=255)
                                                                                aTTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGatat                                                         <  2:163811/149‑1 (MQ=255)
                                                                                    aaTTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGaa                                                     >  4:519066/1‑149 (MQ=255)
                                                                                       ttGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGaaaaa                                                  <  1:182526/149‑1 (MQ=255)
                                                                                         gCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGt                                                >  3:514221/1‑149 (MQ=255)
                                                                                          ccTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTg                                               <  3:679922/149‑1 (MQ=255)
                                                                                            tCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGa                                             <  4:450221/149‑1 (MQ=255)
                                                                                             ccATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGAt                                            >  1:229760/1‑149 (MQ=255)
                                                                                                   gTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTaaa                                      <  3:8697/149‑1 (MQ=255)
                                                                                                   gTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAACGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTaaa                                      <  2:313709/149‑1 (MQ=255)
                                                                                                   gTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAACGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTaaa                                      <  2:106903/149‑1 (MQ=255)
                                                                                                      cAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGa                                                                   <  2:598170/117‑1 (MQ=255)
                                                                                                      cAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGa                                                                   >  1:598170/1‑117 (MQ=255)
                                                                                                        accaccAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCt                                 >  1:305107/1‑149 (MQ=255)
                                                                                                           accaAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCa                              <  1:599444/149‑1 (MQ=255)
                                                                                                            ccaAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAg                             <  3:266866/149‑1 (MQ=255)
                                                                                                            ccaAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAg                             >  2:407673/1‑149 (MQ=255)
                                                                                                             caAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAgg                            <  1:361368/149‑1 (MQ=255)
                                                                                                                agaagaTGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTa                         <  4:452713/149‑1 (MQ=255)
                                                                                                                    gaTGCTGCCAAGACACTTGAGACGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTc                     <  2:325593/149‑1 (MQ=255)
                                                                                                                             aaGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCaa            >  2:605771/1‑149 (MQ=255)
                                                                                                                              aGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAAc           <  1:373832/149‑1 (MQ=255)
                                                                                                                               gACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAACg          >  4:147031/1‑149 (MQ=255)
                                                                                                                                 cacTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAACGtt        <  2:528513/149‑1 (MQ=255)
                                                                                                                                       aGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAACGTTTctgct  <  3:598543/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
CCAATCTGAACGGCACTGAAATTGCCATTACCTATGTCTACAAAGGTGACAAGGTGCTTAAGCAATCTTCTGAAACCAAAATTCAATTTGCCTCCATTGGTGCAACCACCAAAGAAGATGCTGCCAAGACACTTGAGCCGTTAAGCGCCAAATACAAAAACATCGCGGGTGTTGAAGAAAAATTAACCTATACCGATACCTACGCGCAGGAAAACGTGACTATCGATATGGAAAAAGTGGATTTTAAAGCCCTGCAGGGTATTTCAGGAATCAACGTTTCTGCT  >  NC_000913/2211325‑2211608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: