Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,068 +C intergenic (‑178/‑242) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,0681.C99.5% 687.3 / ‑3.8 188intergenic (‑178/‑242)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (30/157);  total (30/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 207 total aligned reads displayed.

GTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCT  >  NC_000913/4476976‑4477153
                                                                                              |                                                                                     
  gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCATGcgacgcat                                                                                     >  2:4225764/1‑90 (MQ=255)
  gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcat                                                                                     >  2:3810202/1‑90 (MQ=255)
  gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcat                                                                                     <  2:5744246/93‑4 (MQ=255)
   tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcatt                                                                                    <  1:1023975/93‑5 (MQ=255)
   tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcatt                                                                                    >  2:4240224/1‑89 (MQ=255)
      ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                 <  2:2066009/93‑1 (MQ=255)
      ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                 <  2:4432730/93‑1 (MQ=255)
      ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                 <  2:4846359/93‑1 (MQ=255)
       cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                                <  1:751096/93‑1 (MQ=255)
         ttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                                                              <  1:1515051/93‑1 (MQ=255)
          tATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACACGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTg                                                                             <  2:704784/93‑1 (MQ=255)
           aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                                            <  2:1293318/93‑1 (MQ=255)
            tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                                           <  1:2688678/93‑1 (MQ=255)
             cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                                                                          <  2:5831133/93‑1 (MQ=255)
                  gTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                                                                     <  1:347992/93‑1 (MQ=255)
                    aaGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATGACAGGTGATTGCGgcgc                                                                   >  2:5496684/1‑93 (MQ=255)
                      gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                                 >  2:1356006/1‑93 (MQ=255)
                       aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                                                <  2:1407206/93‑1 (MQ=255)
                       aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                                                <  2:1562670/93‑1 (MQ=255)
                       aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                                                >  1:4675604/1‑93 (MQ=255)
                        ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                                                               >  1:2138838/1‑93 (MQ=255)
                        ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                                                               >  1:5804288/1‑93 (MQ=255)
                        ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                                                               <  2:3325921/93‑1 (MQ=255)
                         cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                                                              >  2:2010157/1‑93 (MQ=255)
                          gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGtt                                                             <  2:3404605/93‑1 (MQ=255)
                          ctATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGtt                                                             <  2:4126717/92‑1 (MQ=255)
                           tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                                            <  2:4364537/93‑1 (MQ=255)
                           tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                                            <  2:3202802/93‑1 (MQ=255)
                            aTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAt                                                           <  2:117645/93‑1 (MQ=255)
                             ttCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCGGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                                          <  2:4690631/93‑1 (MQ=255)
                             ttCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGCGATTGCGGCGCAATGTTAtt                                                          >  2:4320729/1‑93 (MQ=255)
                             tcgggcggTTACTGAGGAATGGCCGGGGCGAAAGCAGGAATATTCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                                                         <  1:1940756/90‑1 (MQ=255)
                              tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                                                         <  2:2615909/93‑1 (MQ=255)
                              tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc                                                         <  2:5248705/93‑1 (MQ=255)
                               cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                        >  2:4160923/1‑93 (MQ=255)
                               cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                        <  1:3729666/93‑1 (MQ=255)
                               cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                        <  2:6131499/93‑1 (MQ=255)
                               cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                        <  1:5970376/93‑1 (MQ=255)
                                ggcggATACCTGACGAGTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGCCGCAATGTTATTCAt                                                       <  1:888173/93‑1 (MQ=255)
                                ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                                                       <  2:3641742/93‑1 (MQ=255)
                                 gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa                                                      <  1:5858835/93‑1 (MQ=255)
                                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                                    <  1:1151426/93‑1 (MQ=255)
                                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                                    <  1:3110409/93‑1 (MQ=255)
                                   ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                                    >  1:2926901/1‑93 (MQ=255)
                                    gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                                                   <  2:632219/93‑1 (MQ=255)
                                    gATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGCCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt                                                   <  2:1071149/93‑1 (MQ=255)
                                     aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTTATTGCGGCGCAATGTTATTCATAGttt                                                  <  1:2218652/93‑1 (MQ=255)
                                     aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt                                                  <  1:5950478/93‑1 (MQ=255)
                                     aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt                                                  <  1:403221/93‑1 (MQ=255)
                                     aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt                                                  <  2:5644526/93‑1 (MQ=255)
                                      tACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                                 >  2:1853583/1‑93 (MQ=255)
                                      tACCTGACGAAGGACGCGGGCGATAGCATGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                                 <  2:1707721/93‑1 (MQ=255)
                                       aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                                                <  1:5595089/93‑1 (MQ=255)
                                            aCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaa                                           <  1:4461823/93‑1 (MQ=255)
                                             cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaat                                          <  1:2837071/93‑1 (MQ=255)
                                              gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                                         <  1:2439663/93‑1 (MQ=255)
                                               aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa                                        >  2:2501019/1‑93 (MQ=255)
                                               aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa                                        <  1:1921312/93‑1 (MQ=255)
                                                 tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaa                                      >  2:1336139/1‑93 (MQ=255)
                                                 tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGCTTAACCCaataaaa                                      >  2:2698536/1‑93 (MQ=255)
                                                 agACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaa                                      <  2:1916095/92‑1 (MQ=255)
                                                  gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaa                                     <  2:5929031/93‑1 (MQ=255)
                                                  gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaa                                     <  2:3691452/93‑1 (MQ=255)
                                                    cgcgGGCGATAGCAGGATGATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaata                                   <  2:281692/93‑1 (MQ=255)
                                                      cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaataaa                                 <  2:3825205/93‑1 (MQ=255)
                                                       gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAg                                <  2:2646220/93‑1 (MQ=255)
                                                        ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGt                               <  2:3267947/93‑1 (MQ=255)
                                                         gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa                              <  2:5350005/93‑1 (MQ=255)
                                                           gATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAt                            <  2:845638/93‑1 (MQ=255)
                                                              tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                          <  2:1724160/89‑1 (MQ=255)
                                                              tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                          <  1:639372/89‑1 (MQ=255)
                                                              tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                          <  1:2716913/89‑1 (MQ=255)
                                                              tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                          <  2:5419841/89‑1 (MQ=255)
                                                              tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                          <  1:5294890/89‑1 (MQ=255)
                                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata                        <  1:575632/91‑1 (MQ=255)
                                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata                        <  2:2327506/91‑1 (MQ=255)
                                                                gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata                        >  2:2426114/3‑93 (MQ=255)
                                                                  aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa                      <  1:1766222/93‑1 (MQ=255)
                                                                  aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa                      <  2:2255182/93‑1 (MQ=255)
                                                                   agATTATCCTGCAGCGTCAGGCGACACATTACAGGTGATTGCGGCGCAAAGTAATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAt                     <  1:5465086/92‑1 (MQ=255)
                                                                     aTTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAttt                   <  1:6141098/93‑1 (MQ=255)
                                                                      ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg                  <  1:5715440/93‑1 (MQ=255)
                                                                      ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg                  <  1:133227/93‑1 (MQ=255)
                                                                      ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg                  <  1:407111/93‑1 (MQ=255)
                                                                      ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg                  <  1:4318612/93‑1 (MQ=255)
                                                                       tATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGt                 <  1:5340436/93‑1 (MQ=255)
                                                                       tATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGt                 <  1:251756/93‑1 (MQ=255)
                                                                        aTCCTGCAGCGTCAGGCGACGCATTTCAGGAGATTGCGGCGGACTGTTATGCATAGTTTAACCCAATAAAAATAAAGTAATTAGACCtttttt                <  1:4324528/93‑1 (MQ=255)
                                                                        aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTACCCCAATAAAAATAAAGTAATTATACATTTGtt                <  1:5155837/93‑1 (MQ=255)
                                                                        aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGtt                <  2:454959/93‑1 (MQ=255)
                                                                         tCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGttt               <  1:2050485/93‑1 (MQ=255)
                                                                         tCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGttt               <  2:1965344/93‑1 (MQ=255)
                                                                            tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAAt            >  2:4671240/1‑93 (MQ=255)
                                                                            tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAAt            <  1:5911039/93‑1 (MQ=255)
                                                                             gCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATa           <  2:2752499/93‑1 (MQ=255)
                                                                             gCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATa           <  2:4504035/93‑1 (MQ=255)
                                                                                gCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCa        <  2:4902451/93‑1 (MQ=255)
                                                                                gCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCa        >  1:493928/1‑93 (MQ=255)
                                                                                 cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc       <  2:5128758/93‑1 (MQ=255)
                                                                                     aGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCActctc   <  1:3257064/93‑1 (MQ=255)
                                                                                              |                                                                                     
GTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCT  >  NC_000913/4476976‑4477153

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