Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,908 C→G 10.0% intergenic (+17/+22) yfeN → / ← yfeR putative outer membrane protein/transcriptional regulator of yefH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,9080CG10.0% 307.3 / 21.7 90intergenic (+17/+22)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (43/38);  new base G (5/4);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525768‑2526056
                                                                                                                                            |                                                                                                                                                    
tGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                                              <  2:376008/149‑1 (MQ=255)
   aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTa                                                                                                                                           >  4:726640/1‑149 (MQ=255)
   aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTa                                                                                                                                           <  4:317924/149‑1 (MQ=255)
          tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       <  4:1436687/146‑1 (MQ=255)
          tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       >  3:1436687/1‑146 (MQ=255)
          tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTc                                                                                                                                    >  4:593165/1‑149 (MQ=255)
                          acaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTAt                                                                                                                                          >  4:1310788/1‑127 (MQ=255)
                          acaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTAt                                                                                                                                          <  3:1310788/127‑1 (MQ=255)
                          acaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATAc                                                                                                                    >  2:1036582/1‑149 (MQ=255)
                            aTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  1:1156008/144‑1 (MQ=255)
                            aTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       >  2:1156008/1‑144 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGGTTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                             >  4:367897/1‑117 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGGTTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                             <  3:367897/117‑1 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       >  1:486602/1‑123 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       <  2:486602/123‑1 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAc                                                                                                                                  >  2:639645/1‑128 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAc                                                                                                                                  <  1:639645/128‑1 (MQ=255)
                                    ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       <  1:493397/120‑1 (MQ=255)
                                    ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                       >  2:493397/1‑120 (MQ=255)
                                    ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                          <  1:1036582/149‑1 (MQ=255)
                                         ttCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGAt                                                                                                             >  3:115748/1‑141 (MQ=255)
                                         ttCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGAt                                                                                                             <  4:115748/141‑1 (MQ=255)
                                         ttCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                          <  2:825822/144‑1 (MQ=255)
                                         ttCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                          >  1:825822/1‑144 (MQ=255)
                                          tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                          <  4:1101386/143‑1 (MQ=255)
                                          tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                          >  3:1101386/1‑143 (MQ=255)
                                          tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                    <  3:593165/149‑1 (MQ=255)
                                             gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                                              <  1:616525/104‑1 (MQ=255)
                                             gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                                              >  2:616525/1‑104 (MQ=255)
                                                tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                         >  2:554874/1‑122 (MQ=255)
                                                tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                         <  1:554874/122‑1 (MQ=255)
                                                                 aaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTcc                                                                                                     <  1:263519/125‑1 (MQ=255)
                                                                 aaCTGGGCGAAGATGTTTACGAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTcc                                                                                                     >  2:263519/1‑125 (MQ=255)
                                                                  aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                                 >  2:212528/1‑112 (MQ=255)
                                                                  aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                                 <  1:212528/112‑1 (MQ=255)
                                                                    tGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCgg                                                                                                                                                >  2:55140/1‑79 (MQ=255)
                                                                    tGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCgg                                                                                                                                                <  1:55140/79‑1 (MQ=255)
                                                                     gggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGt                                                                         <  3:726640/149‑1 (MQ=255)
                                                                         gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  2:645112/99‑1 (MQ=255)
                                                                         gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       >  1:645112/1‑99 (MQ=255)
                                                                         gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                                            <  3:300774/110‑1 (MQ=255)
                                                                         gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                                            >  4:300774/1‑110 (MQ=255)
                                                                          aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCgg                                                                                                                                                >  2:322696/1‑73 (MQ=255)
                                                                          aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCgg                                                                                                                                                <  1:322696/73‑1 (MQ=255)
                                                                                           gCGATAATCTACTCCATCAAAGTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACgg                                                                                                               <  3:539386/89‑1 (MQ=255)
                                                                                           gCGATAATCTACTCCATCAAAGTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACgg                                                                                                               >  4:539386/1‑89 (MQ=255)
                                                                                             gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATa                                                                                                                     >  4:17612/1‑81 (MQ=255)
                                                                                             gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATa                                                                                                                     <  3:17612/81‑1 (MQ=255)
                                                                                             gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                 >  3:1128202/1‑149 (MQ=255)
                                                                                             gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                 >  3:63941/1‑149 (MQ=255)
                                                                                                       tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGt                                       <  4:1128202/149‑1 (MQ=255)
                                                                                                           ttaaaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                             <  2:727832/105‑1 (MQ=255)
                                                                                                           ttaaaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                             >  1:727832/3‑107 (MQ=255)
                                                                                                           tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                           >  1:1267444/1‑125 (MQ=255)
                                                                                                           tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                           <  2:1267444/125‑1 (MQ=255)
                                                                                                           tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                   >  3:1367207/1‑149 (MQ=255)
                                                                                                           ccAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                        >  2:695519/2‑96 (MQ=255)
                                                                                                           ccAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                        <  1:695519/95‑1 (MQ=255)
                                                                                                            cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgc                                                                                      >  1:583889/1‑97 (MQ=255)
                                                                                                            cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgc                                                                                      <  2:583889/97‑1 (MQ=255)
                                                                                                            cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  4:1367207/149‑1 (MQ=255)
                                                                                                                tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGTTACATCGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACTTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  2:1046768/145‑1 (MQ=255)
                                                                                                                tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  >  1:1046768/1‑145 (MQ=255)
                                                                                                                tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  4:222102/145‑1 (MQ=255)
                                                                                                                tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  >  3:222102/1‑145 (MQ=255)
                                                                                                                      ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        >  3:1307471/1‑149 (MQ=255)
                                                                                                                      ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        >  1:697905/1‑149 (MQ=255)
                                                                                                                      ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        >  2:1151138/1‑149 (MQ=255)
                                                                                                                      ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        >  4:726205/1‑149 (MQ=255)
                                                                                                                          gattgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGt                                       >  3:1447027/1‑130 (MQ=255)
                                                                                                                          gattgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGt                                       <  4:1447027/130‑1 (MQ=255)
                                                                                                                           attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCCCGTACCACGTCCCATAAAGCTTCa                                               <  1:610245/121‑1 (MQ=255)
                                                                                                                           attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCCCGTACCACGTCCCATAAAGCTTCa                                               >  2:610245/1‑121 (MQ=255)
                                                                                                                            ttgattACTGCCGGGGCTCTCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt                  >  1:1268174/1‑149 (MQ=255)
                                                                                                                               attACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                                >  4:1220480/1‑68 (MQ=255)
                                                                                                                               attACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                                <  3:1220480/68‑1 (MQ=255)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                            <  2:769812/86‑1 (MQ=39)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                            >  1:769812/2‑87 (MQ=38)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                               <  2:282054/115‑1 (MQ=255)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                               <  4:700885/115‑1 (MQ=255)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                               >  1:282054/2‑116 (MQ=255)
                                                                                                                                atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                               >  3:700885/2‑116 (MQ=255)
                                                                                                                                  aCTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTGTCGTACCACGTCCCATaaa                                                     <  3:864807/108‑1 (MQ=255)
                                                                                                                                  aCTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                     >  4:864807/1‑108 (MQ=255)
                                                                                                                                  aCTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtctc                                                                      <  1:223984/91‑1 (MQ=255)
                                                                                                                                  aCTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtctc                                                                      >  2:223984/1‑91 (MQ=255)
                                                                                                                                    tGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                            >  2:489626/1‑131 (MQ=255)
                                                                                                                                        ggggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACTTCCCATAAAGCTTCAGCAGCGGTGGAAAGTGAACGGTTCTTCCTTCGCACCAGCAt      >  4:623510/1‑149 (MQ=255)
                                                                                                                                        ggggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                             <  2:166494/126‑1 (MQ=255)
                                                                                                                                        ggggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                             >  1:166494/1‑126 (MQ=255)
                                                                                                                                            cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAc  <  4:1307471/149‑1 (MQ=255)
                                                                                                                                            cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAc  <  1:1151138/149‑1 (MQ=255)
                                                                                                                                            cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAc  >  2:819202/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525768‑2526056

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