Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,413 T→G 12.8% intergenic (‑12/+102) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4130TG12.8% 434.8 / 59.0 148intergenic (‑12/+102)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (57/72);  new base G (9/10);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 150 total aligned reads displayed.

CCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070265‑4070559
                                                                                                                                                    |                                                                                                                                                   
 ccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                    >  4:813416/1‑149 (MQ=255)
  cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGa                                                                                                                                                   <  2:290702/149‑1 (MQ=255)
   tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                  >  2:451219/1‑149 (MQ=255)
   tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                  >  1:581373/1‑149 (MQ=255)
    gAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                 >  2:1143807/1‑149 (MQ=255)
    gAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGc                                                                                                                                                <  3:983658/149‑1 (MQ=255)
        ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                                                                            <  2:581373/149‑1 (MQ=255)
         gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                           <  1:1000629/149‑1 (MQ=255)
         gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                           <  3:813416/149‑1 (MQ=255)
         gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                           >  3:279210/1‑149 (MQ=255)
          ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATg                                                                                                                                          >  1:122079/1‑149 (MQ=255)
          ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtat                                                                                                                                           <  1:1148933/149‑4 (MQ=255)
            aaTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTc                                                                                                                                        <  2:1623362/149‑1 (MQ=255)
             aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCa                                                                                                                                       >  2:609071/1‑149 (MQ=255)
               tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                                                                            >  3:1575583/1‑142 (MQ=255)
               tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                                                                            <  4:1575583/142‑1 (MQ=255)
               tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                     <  1:451219/149‑1 (MQ=255)
                  gCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                      <  3:1652001/145‑1 (MQ=255)
                  gCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                      >  4:1652001/1‑145 (MQ=255)
                     aGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtatgt                                                                                                                                         >  3:722497/1‑135 (MQ=255)
                     aGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtatgt                                                                                                                                         <  4:722497/140‑6 (MQ=255)
                       aCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaaaa                                                                                                                             <  3:1452769/149‑1 (MQ=255)
                         cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACt                                                                                                                           <  3:563184/149‑1 (MQ=255)
                         cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtat                                                                                                                                           <  2:809780/134‑4 (MQ=255)
                             aCCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCACACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCAc                                                                                                                       <  4:1098616/149‑1 (MQ=255)
                                 aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCAGATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGa                                                                                                                   <  3:1112659/149‑1 (MQ=255)
                                     tCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGt                                                                                                               >  2:446679/1‑149 (MQ=255)
                                       acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                                                                             >  4:1189906/1‑149 (MQ=255)
                                                   tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                 <  3:177940/149‑1 (MQ=255)
                                                    gATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                      >  3:895819/1‑111 (MQ=255)
                                                      tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  1:1554519/149‑1 (MQ=255)
                                                        aTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATCCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                            <  4:1456943/149‑1 (MQ=255)
                                                        aTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               <  3:1860520/146‑1 (MQ=255)
                                                        aTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               >  4:1860520/1‑146 (MQ=255)
                                                         tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCt                                                                                                                                               <  2:1719270/98‑1 (MQ=255)
                                                         tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCt                                                                                                                                               >  4:1325601/1‑98 (MQ=255)
                                                         tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCt                                                                                                                                               <  3:1325601/98‑1 (MQ=255)
                                                          cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                          <  4:279210/149‑1 (MQ=255)
                                                           ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                         <  1:609071/149‑1 (MQ=255)
                                                           ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                         <  3:649794/149‑1 (MQ=255)
                                                            cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        <  1:446679/149‑1 (MQ=255)
                                                            cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        <  3:386659/149‑1 (MQ=255)
                                                             ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                     <  3:1077290/103‑1 (MQ=255)
                                                             ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               <  2:842308/141‑1 (MQ=255)
                                                             ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                       >  3:976341/1‑149 (MQ=255)
                                                                  gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             <  3:273154/138‑1 (MQ=255)
                                                                  gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             >  4:273154/1‑138 (MQ=255)
                                                                   cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               >  1:1697478/1‑135 (MQ=255)
                                                                    aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              >  3:1158170/1‑135 (MQ=255)
                                                                    aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  4:1158170/135‑1 (MQ=255)
                                                                              gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                          >  4:1021678/1‑129 (MQ=255)
                                                                                       cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACg                                                                                                                 >  2:294/1‑97 (MQ=255)
                                                                                          cAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                                          <  3:46107/149‑1 (MQ=255)
                                                                                              gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGt                                                                                                                                         >  4:1842008/1‑66 (MQ=255)
                                                                                                 tatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                   <  3:1435263/117‑1 (MQ=255)
                                                                                                   ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                 <  4:467643/101‑1 (MQ=255)
                                                                                                   ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                 >  3:467643/1‑101 (MQ=255)
                                                                                                   ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                   >  4:931616/1‑115 (MQ=255)
                                                                                                   ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                   <  3:931616/115‑1 (MQ=255)
                                                                                                       tGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  >  4:496231/1‑112 (MQ=255)
                                                                                                       tGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  3:496231/112‑1 (MQ=255)
                                                                                                         aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                   <  3:108441/61‑1 (MQ=255)
                                                                                                         aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                   >  4:108441/1‑61 (MQ=255)
                                                                                                         aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               >  1:70102/1‑97 (MQ=255)
                                                                                                         aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               <  2:70102/97‑1 (MQ=255)
                                                                                                         aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  >  2:658289/1‑110 (MQ=255)
                                                                                                           ttGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        <  2:1478631/102‑1 (MQ=255)
                                                                                                                tACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                     >  1:1081823/1‑52 (MQ=255)
                                                                                                                tACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                     <  2:1081823/52‑1 (MQ=255)
                                                                                                                 aCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                                   >  4:405861/1‑149 (MQ=255)
                                                                                                                         ccTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                            <  1:225893/100‑1 (MQ=255)
                                                                                                                          cTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                                                                              <  1:832324/97‑1 (MQ=255)
                                                                                                                          cTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTAAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                          <  3:405861/149‑1 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              >  1:391941/1‑80 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  2:420197/80‑1 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  2:391941/80‑1 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  4:1390591/80‑1 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              >  3:1390591/1‑80 (MQ=255)
                                                                                                                           ttaattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              >  1:420197/1‑80 (MQ=255)
                                                                                                                                     aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                             >  4:1421546/1‑87 (MQ=255)
                                                                                                                                       catcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        >  1:762484/4‑75 (MQ=25)
                                                                                                                                       aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCATCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  2:1139351/1‑149 (MQ=255)
                                                                                                                                       aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  4:799041/1‑149 (MQ=255)
                                                                                                                                       aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  3:1008587/1‑149 (MQ=255)
                                                                                                                                        atcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                                                                              <  1:187951/82‑1 (MQ=25)
                                                                                                                                        atcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                                                                              >  2:187951/3‑84 (MQ=25)
                                                                                                                                        atcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTaa                                 <  4:937423/127‑1 (MQ=38)
                                                                                                                                        atcAGCCTCACGGGAGGCTGATGTCAAACGCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGAAGAGAGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa            >  3:1787479/3‑149 (MQ=39)
                                                                                                                                         tcAGCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg               >  3:1002004/2‑145 (MQ=255)
                                                                                                                                         tcAGCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg               <  4:1002004/144‑1 (MQ=255)
                                                                                                                                           aGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCg         >  1:623167/1‑149 (MQ=255)
                                                                                                                                            gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg                                                                          >  1:964406/1‑83 (MQ=255)
                                                                                                                                            gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg                                                                          <  2:964406/83‑1 (MQ=255)
                                                                                                                                            gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      >  1:654424/1‑87 (MQ=255)
                                                                                                                                             ccTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       >  4:665864/1‑149 (MQ=255)
                                                                                                                                              cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                      <  4:339713/69‑1 (MQ=255)
                                                                                                                                              cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTg      >  3:187089/1‑149 (MQ=255)
                                                                                                                                               tCCCGTGGGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTgg     <  3:799041/149‑1 (MQ=255)
                                                                                                                                                cccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                        >  1:824842/1‑129 (MQ=255)
                                                                                                                                                cccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagca                        <  2:824842/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070265‑4070559

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