Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,077 G→A 6.5% intergenic (‑113/+15) mfd ← / ← ycfT transcription‑repair coupling factor/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,0770GA 6.5% 793.9 / 28.8 245intergenic (‑113/+15)mfd/ycfTtranscription‑repair coupling factor/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (120/109);  new base A (9/7);  total (129/116)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 249 total aligned reads displayed.

TTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCA  >  NC_000913/1173929‑1174225
                                                                                                                                                    |                                                                                                                                                    
tttGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCg                                                                                                                                                      >  2:1212230/1‑149 (MQ=255)
     cGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                 <  1:201255/149‑1 (MQ=255)
      gggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTggg                                                                                                                                                <  3:1455469/149‑1 (MQ=255)
      gggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATCGTTTAGCAAATTGAATCGCCCCGGCt                                                                                                                                                   >  4:759747/1‑146 (MQ=255)
           gCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGAt                                                                                                                                           <  1:761295/149‑1 (MQ=255)
              tataACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGc                                                                                                                                        <  4:969556/149‑1 (MQ=255)
                 aaCGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgggg                                                                                                                                               <  4:782636/139‑1 (MQ=255)
                        ttGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAg                                                                                                                              <  3:1669124/149‑1 (MQ=255)
                          gTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTA                                                                                                                            <  1:980575‑M1/149‑1 (MQ=255)
                           ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAA                                                                                                                           >  1:477677‑M1/1‑149 (MQ=255)
                           ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAA                                                                                                                           >  1:391009‑M1/1‑149 (MQ=255)
                             cAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                         >  2:946156/1‑133 (MQ=255)
                             cAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACT                                                                                                                         >  2:917936‑M1/1‑149 (MQ=255)
                              aGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTg                                                                                                                        >  2:1399141/1‑149 (MQ=255)
                                  aTAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACgg                                                                                                                    <  1:947965/149‑1 (MQ=255)
                                    aaCGGTTGTCAGAATCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCCTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCTGGGGCGATGGCTTAGCGGGAGTAA                                                                                                                           <  2:450546‑M1/140‑1 (MQ=255)
                                     aCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAg                                                                                                                 <  4:1225935/149‑1 (MQ=255)
                                       ggTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCg                                                                                                               >  2:554890/1‑149 (MQ=255)
                                        gTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCgg                                                                                                              <  2:533363/149‑1 (MQ=255)
                                           gTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCa                                                                                                           <  4:231524/149‑1 (MQ=255)
                                            tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAg                                                                                                          >  2:8759/1‑149 (MQ=255)
                                            tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAg                                                                                                          >  3:977880/1‑149 (MQ=255)
                                             cAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGt                                                                                                         >  1:127483/1‑149 (MQ=255)
                                               gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTg                                                                                                                                                  >  4:1401127/1‑106 (MQ=255)
                                                 ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                     >  1:1665332/1‑149 (MQ=255)
                                                 ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                     >  3:982809/1‑149 (MQ=255)
                                                 ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                     >  4:667227/1‑149 (MQ=255)
                                                 ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTCTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                     <  1:295554/149‑1 (MQ=255)
                                                       ttAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgggg                                                                                               <  1:1170919/149‑1 (MQ=255)
                                                              atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGa                                                                                                                               >  3:619219/1‑110 (MQ=255)
                                                              atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGa                                                                                                                               <  4:619219/110‑1 (MQ=255)
                                                              atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                              <  2:365036/143‑1 (MQ=255)
                                                              atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                              >  1:365036/1‑143 (MQ=255)
                                                              atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                              <  3:761003/143‑1 (MQ=255)
                                                                 ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGaa                                                                                                      <  3:1185672/132‑1 (MQ=255)
                                                                  cTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                >  2:1124252/1‑137 (MQ=255)
                                                                  cTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGa                                                                                    <  1:194018/149‑1 (MQ=255)
                                                                   tCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                   >  4:103850/1‑149 (MQ=255)
                                                                   tCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                   <  3:667227/149‑1 (MQ=255)
                                                                     aaCATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                         <  3:1287919/93‑1 (MQ=255)
                                                                       cATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCt                                                                               <  4:352592/149‑1 (MQ=255)
                                                                       cATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGAGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCCGGGAGAACAACAGATCGCt                                                                               <  4:1431967/149‑1 (MQ=255)
                                                                          atGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGaa                                                                            <  4:977880/149‑1 (MQ=255)
                                                                            gggggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATg                                                                          >  3:747758/1‑149 (MQ=255)
                                                                              gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                <  3:42695/125‑1 (MQ=255)
                                                                              gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                >  4:42695/1‑125 (MQ=255)
                                                                               ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCAc                                                                       >  3:1167200/1‑149 (MQ=255)
                                                                                  atatCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGc                                                                                                                <  2:1659562/105‑1 (MQ=255)
                                                                                     tctcTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCt                                                                               <  3:560172/135‑1 (MQ=255)
                                                                                     tctcTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCt                                                                               >  4:560172/1‑135 (MQ=255)
                                                                                       tctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGaa                                                               >  1:1314400/1‑149 (MQ=255)
                                                                                        cTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAa                                                              >  4:206695‑M2/1‑148 (MQ=255)
                                                                                         tgctAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCgg                                                                                                                                 >  4:129370/4‑81 (MQ=255)
                                                                                         tGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAag                                                             <  3:950744‑M2/149‑3 (MQ=255)
                                                                                         tGATAAGCAATGCTGTTCATTATCCGTAAGCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAag                                                             <  1:1365592‑M2/149‑3 (MQ=255)
                                                                                             aaGCACTGCTGTTCATTATCCGTAATCTACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                         <  4:1457558/149‑1 (MQ=255)
                                                                                               gCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCggg                                                                                                                                                >  2:1674730/1‑60 (MQ=255)
                                                                                                  aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAc                                                    <  4:1564120/149‑1 (MQ=255)
                                                                                                  aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAc                                                    >  1:912735/1‑149 (MQ=255)
                                                                                                       gTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCa                                               >  4:241137/1‑149 (MQ=255)
                                                                                                        ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAagtt                                                           >  3:1302277‑M2/1‑132 (MQ=255)
                                                                                                        ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCaa                                              <  1:779612/149‑1 (MQ=255)
                                                                                                        ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCaa                                              <  3:450630/149‑1 (MQ=255)
                                                                                                           aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCgg                                                                                                                                 >  1:531991/1‑63 (MQ=255)
                                                                                                           aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                           >  2:1307302/1‑149 (MQ=255)
                                                                                                           aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                           <  3:1593274/149‑1 (MQ=255)
                                                                                                            ttATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGccc                                          >  3:1254431/1‑149 (MQ=255)
                                                                                                               tCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCaa                                              >  1:431072/1‑142 (MQ=255)
                                                                                                                 cgcaaTCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc                                     >  3:905405/4‑149 (MQ=255)
                                                                                                                  gTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                 <  2:961191/88‑1 (MQ=255)
                                                                                                                     aTCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                 >  1:1012655/1‑149 (MQ=255)
                                                                                                                         aCATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGa                             >  4:335658/1‑149 (MQ=255)
                                                                                                                          cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                >  1:142768/1‑81 (MQ=255)
                                                                                                                          cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGAt                            <  2:477677/149‑1 (MQ=255)
                                                                                                                          cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGAt                            <  2:225050/149‑1 (MQ=255)
                                                                                                                             ggTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGaa                         <  3:122854/149‑1 (MQ=255)
                                                                                                                             ggTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGaa                         <  3:241137/149‑1 (MQ=255)
                                                                                                                                  aGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAAc                    <  4:1254431/149‑1 (MQ=255)
                                                                                                                                   gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACt                   >  2:1028570/1‑149 (MQ=255)
                                                                                                                                   gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACt                   >  3:1445186/1‑149 (MQ=255)
                                                                                                                                     aaaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                 <  2:1610646/149‑1 (MQ=255)
                                                                                                                                     aaaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                 >  4:213176/1‑149 (MQ=255)
                                                                                                                                      aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAATAAACTGAc                >  2:1479995/1‑149 (MQ=255)
                                                                                                                                       aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa               >  1:1084362/1‑149 (MQ=255)
                                                                                                                                         tGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACaaa             <  4:11779/149‑1 (MQ=255)
                                                                                                                                          gAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGGACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAg            <  3:366941/149‑1 (MQ=255)
                                                                                                                                          ccaTCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAg            >  4:1618567/3‑149 (MQ=255)
                                                                                                                                            aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg          <  2:1721674/149‑1 (MQ=255)
                                                                                                                                            aTCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg          >  1:890290/1‑149 (MQ=255)
                                                                                                                                             tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa         <  1:700610/149‑1 (MQ=255)
                                                                                                                                             tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa         >  2:274701/1‑149 (MQ=255)
                                                                                                                                             tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa         <  4:103414/149‑1 (MQ=255)
                                                                                                                                             tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa         >  2:1237450/1‑149 (MQ=255)
                                                                                                                                               gCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATa       <  3:213176/149‑1 (MQ=255)
                                                                                                                                                ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATaa      >  1:1608405/1‑149 (MQ=255)
                                                                                                                                                ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATaa      >  1:1502723/1‑149 (MQ=255)
                                                                                                                                                ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCACTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGAACAGATGAACAAACTGACAAAGGGATaa      >  3:1059026/1‑149 (MQ=255)
                                                                                                                                                  ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAAc    <  1:1692021/149‑1 (MQ=255)
                                                                                                                                                  ccGGCTGGGGCGATGGCTTAACGGGAGTAACTGACGGCAGCGAGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAGGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAAc    >  4:116184/1‑149 (MQ=255)
                                                                                                                                                    ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACca  <  1:554890/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCA  >  NC_000913/1173929‑1174225

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: