Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,463,325 A→G D2G (GAC→GGC)  yfdC → putative transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,463,3250AG94.1% 901.6 / 33.5 321D2G (GAC→GGC) yfdCputative transport
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base G (132/169);  total (141/180)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 330 total aligned reads displayed.

TCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACC  >  NC_000913/2463235‑2463414
                                                                                          |                                                                                         
  tACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAAc                                                                                       <  1:645996/93‑1 (MQ=255)
      ttttCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACa                                                                                   >  2:7796424/1‑93 (MQ=255)
      ttttCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACa                                                                                   >  2:2625027/1‑93 (MQ=255)
         tCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACaaaa                                                                                <  2:6686083/93‑1 (MQ=255)
          cAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAAt                                                                               <  1:7492989/93‑1 (MQ=255)
           agagTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAAtt                                                                              >  1:8682729/1‑93 (MQ=255)
           agagTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAAtt                                                                              <  1:681905/93‑1 (MQ=255)
            gagTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTg                                                                             >  1:8543570/1‑93 (MQ=255)
              gTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGAt                                                                           >  2:4566274/1‑93 (MQ=255)
                aaCTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCa                                                                         <  1:9999324/93‑1 (MQ=255)
                 aCTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCaa                                                                        >  2:368741/1‑93 (MQ=255)
                 aCTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCaa                                                                        <  1:4208616/93‑1 (MQ=255)
                  cTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAac                                                                       >  2:7317109/1‑93 (MQ=255)
                   tCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAaca                                                                      >  1:1534958/1‑93 (MQ=255)
                    cccGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAacac                                                                     <  2:884152/93‑1 (MQ=255)
                    cccGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAacac                                                                     <  2:4234913/93‑1 (MQ=255)
                      cGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAg                                                                   >  1:1089137/1‑93 (MQ=255)
                      cGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAg                                                                   >  2:9288969/1‑93 (MQ=255)
                        tCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACAGAGcg                                                                 <  2:6493144/93‑1 (MQ=255)
                        tCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGcg                                                                 >  1:4888773/1‑93 (MQ=255)
                        tCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGcg                                                                 >  2:9460263/1‑93 (MQ=255)
                         cAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGcga                                                                >  2:1824806/1‑93 (MQ=255)
                           gTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGcgacg                                                              >  1:8160520/1‑93 (MQ=255)
                              gCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGaaa                                                           <  2:8996374/93‑1 (MQ=255)
                              gCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGACTAACACAGcgacgcat                                                           >  1:8387420/1‑90 (MQ=255)
                               cGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAAt                                                          <  1:7559302/93‑1 (MQ=255)
                               cGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAAt                                                          <  1:409313/93‑1 (MQ=255)
                                gAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAAtt                                                         >  2:3002564/1‑93 (MQ=255)
                                gAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAAtt                                                         >  2:3591614/1‑93 (MQ=255)
                                gAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAAtt                                                         <  2:7039901/93‑1 (MQ=255)
                                gAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAANNGATCAACACAGCGACGAANtt                                                         >  1:781405/1‑93 (MQ=255)
                                  gCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGa                                                       <  1:2317336/93‑1 (MQ=255)
                                  gCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGa                                                       >  2:7901740/1‑93 (MQ=255)
                                    aaaaaaGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAg                                                     <  1:5448716/93‑1 (MQ=255)
                                    aaaaaaGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAg                                                     >  2:5016879/1‑93 (MQ=255)
                                     aaaaaGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGt                                                    <  1:1411179/93‑1 (MQ=255)
                                      aaaaGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGtt                                                   <  2:53954/93‑1 (MQ=255)
                                        aaGCCGCTACGCTTTAGCTATACGTGCTAATTCAATAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTga                                                 >  1:1848409/1‑93 (MQ=255)
                                        aaGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTga                                                 <  2:2463246/93‑1 (MQ=255)
                                         aGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTgag                                                <  1:2406295/93‑1 (MQ=255)
                                         aGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTgag                                                <  1:9904681/93‑1 (MQ=255)
                                          gCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGACGTTgaga                                               >  2:4813513/1‑93 (MQ=255)
                                           ccGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGagag                                              >  2:9747955/1‑93 (MQ=255)
                                            cgcgACGCTTTAGATATACGTTCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGc                                             <  1:6204704/93‑1 (MQ=255)
                                            cGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGc                                             <  2:3371107/93‑1 (MQ=255)
                                              cTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCga                                           <  1:8906460/93‑1 (MQ=255)
                                               tACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCgaa                                          <  1:3093279/93‑1 (MQ=255)
                                               tACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCgaa                                          >  1:3477588/1‑93 (MQ=255)
                                               tACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCgaa                                          >  2:6318239/1‑93 (MQ=255)
                                                  gCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCgaagaa                                       <  1:4481215/93‑1 (MQ=255)
                                                   cTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGaaa                                      <  2:1753829/93‑1 (MQ=255)
                                                    tttAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGaaaa                                     <  1:5108058/93‑1 (MQ=255)
                                                     ttAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGaaaaa                                    <  1:2171975/93‑1 (MQ=255)
                                                       aGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAaga                                  <  2:1037218/93‑1 (MQ=255)
                                                       aGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAaga                                  >  2:5931760/1‑93 (MQ=255)
                                                       aGCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAaga                                  <  1:8192087/93‑1 (MQ=255)
                                                        gCTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAagag                                 <  2:3852726/93‑1 (MQ=255)
                                                         cTATACGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAgc                                >  2:3135562/1‑93 (MQ=255)
                                                         cTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAgc                                >  1:1708399/1‑93 (MQ=255)
                                                           ataCGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAgcgc                              >  1:3201095/1‑93 (MQ=255)
                                                            taCGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGcgcg                             >  1:1181066/1‑93 (MQ=255)
                                                            taCGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGcgcg                             <  1:6042646/93‑1 (MQ=255)
                                                            taCGTGCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGcgcg                             <  1:9553083/93‑1 (MQ=255)
                                                                 gCTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaa                        >  1:7404495/1‑93 (MQ=255)
                                                                  cTAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaa                       >  1:329909/1‑93 (MQ=255)
                                                                   tAATTCAAGAGAAGAGACCATGGTCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaa                      <  2:5517380/93‑1 (MQ=255)
                                                                   tAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaa                      >  2:1455138/1‑93 (MQ=255)
                                                                   tAATTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaa                      <  2:3413757/93‑1 (MQ=255)
                                                                    aaTTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGCGAGCGAAGAAAAAGAGAGCGTCaaaaaa                     >  1:6172263/1‑93 (MQ=255)
                                                                    aaTTCAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaaa                     >  1:9019836/1‑93 (MQ=255)
                                                                     aTTCAAGAGAAGAGACCATGGGCAACGACAAAAATGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaaaa                    <  1:9726532/93‑1 (MQ=255)
                                                                        cAAGAGAAGAGACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatag                 >  2:166217/1‑93 (MQ=255)
                                                                         aagagaagagACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaataga                <  1:2642687/93‑1 (MQ=255)
                                                                         aagagaagagACCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaataga                >  2:2820991/1‑93 (MQ=255)
                                                                          agagaagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaa               >  1:9814587/1‑93 (MQ=255)
                                                                          agagaagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaa               <  1:7695792/93‑1 (MQ=255)
                                                                           gagaagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaaa              <  2:7698205/93‑1 (MQ=255)
                                                                             gaagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaaata            >  2:754277/1‑93 (MQ=255)
                                                                             gaagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaaata            <  2:5192746/93‑1 (MQ=255)
                                                                              aagagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaaatag           <  2:1129477/93‑1 (MQ=255)
                                                                               agagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAaatagaaataga          <  1:3893184/93‑1 (MQ=255)
                                                                                gagaCCATGGGCAACGACAAAATTGATCAACACAGCGACGCAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGAt         >  1:91842/1‑93 (MQ=255)
                                                                                  gaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGa       <  1:8560796/93‑1 (MQ=255)
                                                                                  gaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGa       <  1:3615611/93‑1 (MQ=255)
                                                                                  gaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGa       >  2:4364808/1‑93 (MQ=255)
                                                                                  gaCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGa       <  2:409357/93‑1 (MQ=255)
                                                                                   aCCATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGaa      <  2:1350861/93‑1 (MQ=255)
                                                                                    ccATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAg     <  1:4424979/93‑1 (MQ=255)
                                                                                    ccATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAg     <  1:5857359/93‑1 (MQ=255)
                                                                                    ataTGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAg     <  2:2214472/91‑1 (MQ=255)
                                                                                     cATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGa    <  2:8235715/93‑1 (MQ=255)
                                                                                     cATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGa    <  1:6981330/93‑1 (MQ=255)
                                                                                     cATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGa    <  1:9390475/93‑1 (MQ=255)
                                                                                     cATGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGAGGAAGa    <  2:4054209/93‑1 (MQ=255)
                                                                                      aTGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAc   <  2:8396577/93‑1 (MQ=255)
                                                                                       tGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAcc  >  1:6700509/1‑93 (MQ=255)
                                                                                       tGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAcc  <  2:253148/93‑1 (MQ=255)
                                                                                       tGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAcc  >  2:9123788/1‑93 (MQ=255)
                                                                                       tGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAcc  <  1:6895746/93‑1 (MQ=255)
                                                                                       tGGGCAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGAcc  >  2:1257999/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACC  >  NC_000913/2463235‑2463414

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