Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,009,002 G→T 6.5% P20T (CCA→ACA)  rhtB ← homoserine, homoserine lactone and S‑methyl‑methionine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,009,0020GT 6.5% 274.1 / 7.9 77P20T (CCA→ACA) rhtBhomoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base T (2/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

AGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGG  >  NC_000913/4008858‑4009149
                                                                                                                                                |                                                                                                                                                   
aGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCgcc                                                                                                                                                 >  2:342886/1‑147 (MQ=255)
aGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGAc                                                                                                                                                 <  2:228543/149‑1 (MQ=255)
   gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAgg                                                                                                                                              >  3:261392/1‑149 (MQ=255)
   gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATGGCACCAGAGCCTGGCGACAgg                                                                                                                                              >  4:585239/1‑149 (MQ=255)
    ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGc                                                                                                                                             <  4:527154/149‑1 (MQ=255)
        aaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCGCCAGGCtctg                                                                                                                                         <  1:342886/149‑4 (MQ=255)
         aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCGCCAGAGCCTGTCGCCAGGCtctgg                                                                                                                                        <  3:329542/149‑5 (MQ=255)
              gTCCCCAACCCCACGTCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                   <  4:163803/149‑1 (MQ=255)
              gTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                   <  2:475179/149‑1 (MQ=255)
              gTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                   <  1:530704/149‑1 (MQ=255)
               tCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATc                                                                                                                                  <  2:282306/149‑1 (MQ=255)
                 cccAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGa                                                                                                                                <  4:199312/149‑1 (MQ=255)
                  ccAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGAt                                                                                                                               <  3:108058/149‑1 (MQ=255)
                      ccccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATTGATGTca                                                                                                                           <  2:425449/149‑1 (MQ=255)
                      ccccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTca                                                                                                                           >  4:30444/1‑149 (MQ=255)
                       cccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcag                                                                                                                          >  4:213169/1‑149 (MQ=255)
                         cACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagca                                                                                                                        <  4:464363/149‑1 (MQ=255)
                          aCGCCAACCAGCACAATATGAATCGCCACTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagcag                                                                                                                       >  3:513993/1‑149 (MQ=255)
                            gCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGt                                                                                                                     >  1:98588/1‑149 (MQ=255)
                               aaCCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAgg                                                                                                                  >  3:96526/1‑149 (MQ=255)
                                  cAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCGCCAGGCtctg                                                                                                                                         >  3:525894/1‑120 (MQ=255)
                                  cAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCGCCAGGCtctg                                                                                                                                         <  4:525894/123‑4 (MQ=255)
                                      acaATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAacca                                                                                                           >  2:138311/1‑149 (MQ=255)
                                       caATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccac                                                                                                          <  3:261650/149‑1 (MQ=255)
                                        aaTATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccacc                                                                                                         >  1:20222/1‑149 (MQ=255)
                                        aaTATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccacc                                                                                                         >  2:264546/1‑149 (MQ=255)
                                        aaTATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccacc                                                                                                         >  1:588624/1‑149 (MQ=255)
                                        aaTATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccacc                                                                                                         >  1:540973/1‑149 (MQ=255)
                                         atatGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAAccacca                                                                                                        <  2:98588/149‑1 (MQ=255)
                                           atGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAC‑‑CGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCt                                                                                                    >  1:313371/1‑149 (MQ=255)
                                             gAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCt                                                                                                    >  3:132804/1‑149 (MQ=255)
                                              aaTCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTa                                                                                                   >  3:195487/1‑149 (MQ=255)
                                               aTCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTaa                                                                                                  <  3:83406/149‑1 (MQ=255)
                                                 cGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAgg                                                                                                <  1:571465/149‑1 (MQ=255)
                                                   ccAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTc                                                                                              <  2:506507/149‑1 (MQ=255)
                                                   ccAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTc                                                                                              <  1:244408/149‑1 (MQ=255)
                                                    cAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCa                                                                                             >  3:294359/1‑149 (MQ=255)
                                                            tCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAAc                                                                                     >  4:333456/1‑149 (MQ=255)
                                                              tGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTc                                                                                   >  4:407825/1‑149 (MQ=255)
                                                               gAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTcc                                                                                  <  3:257096/149‑1 (MQ=255)
                                                                aaGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTccc                                                                                 >  3:570253/1‑149 (MQ=255)
                                                                    ccAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtg                                                                             >  4:441801/1‑149 (MQ=255)
                                                                     cAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                            >  2:567854/1‑149 (MQ=255)
                                                                     cAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                            >  1:427588/1‑149 (MQ=255)
                                                                          aTAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTa                                                                       >  3:167084/1‑149 (MQ=255)
                                                                            aGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTAtt                                                                     >  2:258477/1‑149 (MQ=255)
                                                                             gACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTAttt                                                                    >  4:175046/1‑149 (MQ=255)
                                                                              aCGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTAtttt                                                                   <  2:20222/149‑1 (MQ=255)
                                                                                 ccACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtgt                                                                <  1:567854/149‑1 (MQ=255)
                                                                                   aCCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGc                                                              >  2:310714/1‑149 (MQ=255)
                                                                                   aCCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGc                                                              <  1:367045/149‑1 (MQ=255)
                                                                                               tAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTa                                                  <  3:213169/149‑1 (MQ=255)
                                                                                               tAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTa                                                  >  4:210266/1‑149 (MQ=255)
                                                                                                aaCCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTAc                                                 <  1:310714/149‑1 (MQ=255)
                                                                                                   cGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGt                                              >  1:319788/1‑149 (MQ=255)
                                                                                                           aacGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGt                                      <  4:513993/147‑1 (MQ=255)
                                                                                                           aGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGt                                      <  4:96526/149‑1 (MQ=255)
                                                                                                              gAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgt                                   >  4:100940/1‑149 (MQ=255)
                                                                                                              gAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgt                                   <  2:350548/149‑1 (MQ=255)
                                                                                                               aGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgtg                                  <  3:526825/149‑1 (MQ=255)
                                                                                                                ggtggtCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                            >  3:366358/1‑106 (MQ=255)
                                                                                                                ggtggtCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                            <  4:366358/106‑1 (MQ=255)
                                                                                                                   ggtCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAg                              <  4:602884/149‑1 (MQ=255)
                                                                                                                          gtgtTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCt                       >  3:328815/1‑149 (MQ=255)
                                                                                                                             ttGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTc                                                                                              <  2:574465/75‑1 (MQ=255)
                                                                                                                             ttGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTc                                                                                              >  1:574465/1‑75 (MQ=255)
                                                                                                                             ttGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaa                    >  1:599932/1‑149 (MQ=255)
                                                                                                                               gATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaaaa                  >  1:438043/1‑149 (MQ=255)
                                                                                                                                     aCCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGa            >  1:65953/1‑149 (MQ=255)
                                                                                                                                     aCCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACACTAAAAAAAACGa            <  4:570253/149‑1 (MQ=255)
                                                                                                                                       cAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtt          <  1:68483/149‑1 (MQ=255)
                                                                                                                                       cAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtt          >  2:478757/1‑149 (MQ=255)
                                                                                                                                       cAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtt          <  1:478757/149‑1 (MQ=255)
                                                                                                                                         gagCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtttt        >  3:451842/1‑149 (MQ=255)
                                                                                                                                           gCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTAc      <  2:313371/149‑1 (MQ=255)
                                                                                                                                             cTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGt    <  4:328815/149‑1 (MQ=255)
                                                                                                                                             cTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGt    <  1:138311/149‑1 (MQ=255)
                                                                                                                                               ggCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTgg  >  2:426434/1‑149 (MQ=255)
                                                                                                                                               ggCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTgg  >  4:507922/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGG  >  NC_000913/4008858‑4009149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑