Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,311,616 C→T 5.2% intergenic (‑82/+30) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,311,6160CT 5.2% 211.1 / 2.4 77intergenic (‑82/+30)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (31/42);  new base T (2/2);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGT  >  NC_000913/2311468‑2311761
                                                                                                                                                    |                                                                                                                                                 
tttGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGccc                                                                                                                                                   >  1:1246309/1‑149 (MQ=255)
   gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                                <  1:623006/149‑1 (MQ=255)
     tCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAgg                                                                                                                                              >  1:704914/1‑149 (MQ=255)
         aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                          <  1:1108644/149‑1 (MQ=255)
           ccaaaAAAGCAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCCTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCttt                                                                                                                                        <  2:1384331/149‑1 (MQ=255)
             aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGt                                                                                                                                      <  2:863703/149‑1 (MQ=255)
             aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGt                                                                                                                                      >  1:1290685/1‑149 (MQ=255)
              acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                     <  1:386050/149‑1 (MQ=255)
              acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                     <  2:738928/149‑1 (MQ=255)
              acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                     <  2:1534177/149‑1 (MQ=255)
              acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                     <  1:94077/149‑1 (MQ=255)
              acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGAGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCACGCAGGCCCTTTGtt                                                                                                                                     >  1:1021241/1‑149 (MQ=255)
              acaaAGAAGAGAGCGTCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATAATCCTCATGCGAAAGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAAGAATGATAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                     <  2:686536/149‑1 (MQ=255)
               caaTGACGAGCGCGGCCGACATCAGTGCCACGCGGGTAAAGCTTATATCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                    <  2:537932/149‑1 (MQ=255)
               caaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTc                                                                                                                                    <  2:187854/149‑1 (MQ=255)
               caaAGAAGAGCGCGGCAGCCAGCACTGCCACACGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCACGCACGCCCTttcttc                                                                                                                                    >  2:996500/1‑149 (MQ=255)
               aaaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                    <  1:920057/148‑1 (MQ=255)
               aaaaGAAGAGCGCGGCAGACAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCTCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                    <  2:1025570/148‑1 (MQ=255)
                 aagaagAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGa                                                                                                                                  <  2:1340419/149‑1 (MQ=255)
                  agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                                 <  1:1500455/149‑1 (MQ=255)
                  agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTGTGTACGCTGAAAAGAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                                 >  2:297574/1‑149 (MQ=255)
                   tacgagCGCGGCAGCCAGCAGTGCCACGCGTGAAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGGCGCAAGAGTACACCAAAAATCTGAGTTTGTACGCTGAAAACAATGAACAAAGGGCCTGCGGGCCCTTTTTTcattg                                                                                                                                <  1:137211/146‑5 (MQ=255)
                      gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCa                                                                                                                             <  1:813286/149‑1 (MQ=255)
                            ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCgag                                                                                                                       <  1:1313628/149‑1 (MQ=255)
                                gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCTATATCAATCGAGATTa                                                                                                                   >  1:1417082/1‑149 (MQ=255)
                                 ccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAg                                                                                                                  <  2:602689/149‑1 (MQ=255)
                                       gTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCTTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTAGATATCAATCGAGATTAGAACTgg                                                                                                            <  2:1158817/149‑1 (MQ=255)
                                           cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaaa                                                                                                        <  2:3557/149‑1 (MQ=255)
                                                gggTATAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAATCTGACTTTGTACGCTGAAAACTATGAAATAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGTTTAGAACTGGTAAACCAGa                                                                                                   <  2:279254/149‑1 (MQ=255)
                                                gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                          >  2:804738/1‑126 (MQ=255)
                                                gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                          <  1:804738/126‑1 (MQ=255)
                                                        cTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGCCCCAgcgc                                                                                           >  2:275155/1‑149 (MQ=255)
                                                            tttCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACg                                                                                       >  1:1397978/1‑149 (MQ=255)
                                                             ttCCATACGTGATTTTCCTCATGCGAAGGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGa                                                                                      <  1:200517/149‑1 (MQ=255)
                                                             ttCCATACGGGATTATCCTCATGGGATCGGTCGCAAGTGTACACCTAAAAACTGAGTTTGTAGGCTGAAATCATTGAAAAAAGGGCCCGGAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGa                                                                                      <  1:206591/149‑1 (MQ=255)
                                                               ccATACGTGATTATCCTCATGCGAACGGGCGCAAGAGTACACCAAAAAAGTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCACTCGAGATTAGAACTGGGAAACCAGACCCAGAGCTACGATg                                                                                    <  1:278547/149‑1 (MQ=255)
                                                                     gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAGGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                              >  2:1109149/1‑149 (MQ=255)
                                                                     gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                              <  1:267726/149‑1 (MQ=255)
                                                                     gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                              <  2:1230932/149‑1 (MQ=255)
                                                                       gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                            >  2:330009/1‑149 (MQ=255)
                                                                       gATTATCCTCATGCGAACGGTCGCAAGAGTACAACAAAAAACTTAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                            <  1:872626/149‑1 (MQ=255)
                                                                         ttATCCTCAAGCGAACGGGCGCAAGAGTACACCAAAAAACTGAGTTTGAACGCTGAAACCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTATACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                          <  2:1021241/149‑1 (MQ=255)
                                                                          tATCCTCATGAGAACGTGCGCAATATTAAACAACACAACTGAGATTGTACGCTGAAAAAAAAGAAAAATTGGCCCGCAGTCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                         <  2:1085788/149‑1 (MQ=255)
                                                                              cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATg                                                                     <  2:1528970/149‑1 (MQ=255)
                                                                               acATGCCAACGGTCGCCAGAGTACCCCAAAAACCTGAGTATCTACGCTGAAAACCATGCAAACAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                    <  1:1144584/148‑1 (MQ=255)
                                                                                 aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCa                                                                  >  2:201703/1‑149 (MQ=255)
                                                                                 aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCa                                                                  <  2:1285200/149‑1 (MQ=255)
                                                                                  tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                 >  1:46711/1‑149 (MQ=255)
                                                                                      aaCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTc                                                             <  2:704914/149‑1 (MQ=255)
                                                                                       aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTATGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTcc                                                            >  1:1319771/1‑148 (MQ=255)
                                                                                           tCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACTATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAg                                                        <  1:1212679/149‑1 (MQ=255)
                                                                                              cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGa                                                     <  2:111178/149‑1 (MQ=255)
                                                                                               aaGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                    >  1:687722/1‑149 (MQ=255)
                                                                                                   gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                <  2:293131/149‑1 (MQ=255)
                                                                                                    tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGTTTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                               <  1:1406832/149‑1 (MQ=255)
                                                                                                     acacCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGCGCCCGCAGGCCCGTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                              <  2:450470/149‑1 (MQ=255)
                                                                                                       acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCACGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                            >  1:578184/1‑149 (MQ=255)
                                                                                                         cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                          >  2:1537331/1‑149 (MQ=255)
                                                                                                          aaaaaaCTGAGTTTGTATGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCAAGCGTCACGAGTGAACTGGTTgtcgt                                         >  1:692632/1‑149 (MQ=255)
                                                                                                                 tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagcag                                  >  1:1302339/1‑149 (MQ=255)
                                                                                                                    gaatGTACGGTGAAACCAATCAAACAAGGGCCCGCAGGCCAATAGTTAGATATCAATCGAGAGTAGAAATGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtg                               <  2:1040095/146‑2 (MQ=255)
                                                                                                                    gTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAACAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCAAGCAGGtt                               >  1:714188/1‑149 (MQ=255)
                                                                                                                          aCGCTTAAACCAAAGAAAACAGGTCCCGCAGGCCCTATGTTCGATATCAATCGACATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGGGTTGATGCCAGCGTCACGAGGGAACTGGTTGTCGTCCAGCAGGTTGAtttt                         <  2:241649/149‑1 (MQ=255)
                                                                                                                            gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGt                       >  1:110587/1‑149 (MQ=255)
                                                                                                                             cTGAAAACAATGAAACAAGGGCCCGCAGTCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGATCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                      <  2:1417082/149‑1 (MQ=255)
                                                                                                                             cTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTGTGTTCGATATCAATCGAGATTAGAACTGTTAACCCAGACCCAGAGCTACGATGTGATGAGTGTTGATGCGAGGGTCACGAGTGACCGGGTTGGCGTCCAGCAGGGTGAGTTTGtt                      >  1:1564630/1‑148 (MQ=255)
                                                                                                                                caaaCAATGATAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGCTTAGAACTGGTATACCAGCCCCAGCTCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCCGGTTGATTTTGTAGTc                   >  2:1024969/2‑149 (MQ=255)
                                                                                                                                  aaCAATGAAAATAGGGCCCGAAGGCCCTTTGTTCGATATCAATCGAGATTATAACTTGTAAACCATACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTTACGAGGGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTaaa                 <  2:1100458/149‑1 (MQ=255)
                                                                                                                                     aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGt              >  2:594271/1‑149 (MQ=255)
                                                                                                                                      aTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTa             >  1:369050/1‑149 (MQ=255)
                                                                                                                                        gAACAAAGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg           >  1:688661/1‑149 (MQ=255)
                                                                                                                                        gAACAAAGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg           >  1:1470413/1‑149 (MQ=255)
                                                                                                                                         aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt          >  1:37927/1‑149 (MQ=255)
                                                                                                                                         aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt          >  1:179298/1‑149 (MQ=255)
                                                                                                                                         aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt          >  1:889449/1‑149 (MQ=255)
                                                                                                                                         aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt          >  1:889448/1‑149 (MQ=255)
                                                                                                                                                ggCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATg   >  2:1309843/1‑149 (MQ=255)
                                                                                                                                                 gCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAACCCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGt  >  1:396402/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGT  >  NC_000913/2311468‑2311761

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