Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,815,038 A→G V163A (GTG→GCG)  yqaB ← putative phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,815,0380AG86.0% 731.6 / 107.6 308V163A (GTG→GCG) yqaBputative phosphatase
Reads supporting (aligned to +/- strand):  ref base A (27/16);  new base G (153/112);  total (180/128)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 328 total aligned reads displayed.

GATAACGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGAC  >  NC_000913/2814946‑2815129
                                                                                            |                                                                                           
  tAACGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACcgca                                                                                           <  2:4888529/93‑1 (MQ=255)
  tAACGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACcgca                                                                                           >  1:2581439/1‑93 (MQ=255)
   aaCGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACcgcac                                                                                          <  2:2334974/93‑1 (MQ=255)
   aaCGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACcgcac                                                                                          >  1:2249012/1‑93 (MQ=255)
     cGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCCAAATCGGCATCTTCAAAGACCgcacac                                                                                        <  2:6420136/93‑1 (MQ=255)
      gCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACt                                                                                       >  1:5297067/1‑93 (MQ=255)
      gCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACt                                                                                       <  1:6242815/93‑1 (MQ=255)
       cTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTg                                                                                      >  2:6534582/1‑93 (MQ=255)
        ttCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGc                                                                                     <  1:1088779/93‑1 (MQ=255)
        ttCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGc                                                                                     >  2:1686061/1‑93 (MQ=255)
         tcactcacAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCg                                                                                    <  1:5713071/93‑1 (MQ=255)
          cactcacAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGt                                                                                   <  2:5417801/93‑1 (MQ=255)
             tcacAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCgg                                                                                <  2:2022315/93‑1 (MQ=255)
             tcacAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCgg                                                                                >  2:6798174/1‑93 (MQ=255)
              cacaGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGt                                                                               <  1:4800167/93‑1 (MQ=255)
               acaGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGtt                                                                              >  1:174130/1‑93 (MQ=255)
               acaGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGtt                                                                              >  2:1435402/1‑93 (MQ=255)
               acaGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGtt                                                                              <  1:5461616/93‑1 (MQ=255)
                caGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTg                                                                             <  1:2510572/93‑1 (MQ=255)
                  gCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCa                                                                           <  1:5822757/93‑1 (MQ=255)
                  gCAAGCGAACATCCACGGCGTCCACGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCa                                                                           >  1:1494621/1‑93 (MQ=255)
                    aaGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACg                                                                         >  1:3978456/1‑93 (MQ=255)
                    aaGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACg                                                                         >  2:1610224/1‑93 (MQ=255)
                     aGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGc                                                                        >  1:4881074/1‑93 (MQ=255)
                      gCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGcc                                                                       <  2:5273716/93‑1 (MQ=255)
                      gCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGcc                                                                       <  1:1548269/93‑1 (MQ=255)
                      gCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGcc                                                                       >  2:2592327/1‑93 (MQ=255)
                      gCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGcc                                                                       <  1:4124036/93‑1 (MQ=255)
                       cGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGccc                                                                      <  1:264031/93‑1 (MQ=255)
                       cGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGccc                                                                      <  1:7297367/93‑1 (MQ=255)
                          aCATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATa                                                                   <  2:6021494/93‑1 (MQ=255)
                            aTCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACg                                                                 >  2:3246430/1‑93 (MQ=255)
                             tCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGc                                                                <  2:4709659/93‑1 (MQ=255)
                             tCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGc                                                                <  1:2342530/93‑1 (MQ=255)
                               cACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTg                                                              >  2:726231/1‑93 (MQ=255)
                               cACGGCGTCCATGCCTGCTGCACCGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTg                                                              >  1:203958/1‑93 (MQ=255)
                                aCGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTgc                                                             <  1:6558507/93‑1 (MQ=255)
                                     gTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGca                                                        <  2:6906798/93‑1 (MQ=255)
                                        cATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCacaa                                                     >  2:5718663/1‑93 (MQ=255)
                                         aTGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCacaac                                                    >  2:2872019/1‑93 (MQ=255)
                                         aTGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCacaac                                                    <  2:3025440/93‑1 (MQ=255)
                                          tGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCacaaca                                                   >  1:4472964/1‑93 (MQ=255)
                                          tGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCacaaca                                                   >  1:3491861/1‑93 (MQ=255)
                                            ccTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACaaa                                                 <  1:1934819/93‑1 (MQ=255)
                                            ccTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACaaa                                                 >  2:1288338/1‑93 (MQ=255)
                                             ctgctgCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACaaaa                                                >  2:7448774/1‑93 (MQ=255)
                                              tgctgcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACaaaaa                                               <  2:4262145/93‑1 (MQ=255)
                                               gctgcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAt                                              <  2:4403300/93‑1 (MQ=255)
                                                 tgcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCCCGCCCATACGCTGCGCGCACAACAAAAAtgt                                            >  1:6672848/1‑93 (MQ=255)
                                                    aCGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTc                                         >  1:3065336/1‑93 (MQ=255)
                                                    aCGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTc                                         >  2:6619544/1‑93 (MQ=255)
                                                    aCGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTc                                         <  2:1955925/93‑1 (MQ=255)
                                                          cgccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGgcgc                                   <  2:4976201/93‑1 (MQ=255)
                                                          cgccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGgcgc                                   >  2:5046511/1‑93 (MQ=255)
                                                           gccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGcgcg                                  >  1:6375918/1‑93 (MQ=255)
                                                           gccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGcgcg                                  >  1:2146951/1‑93 (MQ=255)
                                                            ccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCgg                                 <  2:3816510/93‑1 (MQ=255)
                                                            ccTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCgg                                 <  1:4049628/93‑1 (MQ=255)
                                                             cTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCggg                                >  2:7523448/1‑93 (MQ=255)
                                                             cTGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCggg                                >  1:3766873/1‑93 (MQ=255)
                                                              tGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGt                               <  2:1652846/93‑1 (MQ=255)
                                                              tGAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGt                               >  2:3467277/1‑93 (MQ=255)
                                                               gAATACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGtt                              >  2:945218/1‑93 (MQ=255)
                                                               gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGtt                              >  2:2555543/1‑93 (MQ=255)
                                                                aaTACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGttt                             >  1:6849192/1‑93 (MQ=255)
                                                                  tACCGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAt                           <  1:7056201/93‑1 (MQ=255)
                                                                    ccGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATgg                         <  1:2232690/93‑1 (MQ=255)
                                                                     cGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGgt                        <  2:2404061/93‑1 (MQ=255)
                                                                     cGAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGgt                        >  1:3725839/1‑93 (MQ=255)
                                                                      gAAATCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGgtg                       <  2:35367/93‑1 (MQ=255)
                                                                       aaaTCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGgtgt                      >  2:5801191/1‑93 (MQ=255)
                                                                       aaaTCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGgtgt                      <  2:6180788/93‑1 (MQ=255)
                                                                        aaTCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGtt                     <  1:3153562/93‑1 (MQ=255)
                                                                         aTCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGttt                    >  1:1674450/1‑93 (MQ=255)
                                                                         aTCGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGttt                    <  1:2862104/93‑1 (MQ=255)
                                                                          tCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTg                   >  2:3330849/1‑93 (MQ=255)
                                                                           cGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGa                  <  2:1348869/93‑1 (MQ=255)
                                                                           cGGCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGa                  <  1:7492202/93‑1 (MQ=255)
                                                                             gCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACg                <  2:6263409/93‑1 (MQ=255)
                                                                             gCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACg                >  2:474582/1‑93 (MQ=255)
                                                                             gCATCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACg                >  1:5162197/1‑93 (MQ=255)
                                                                                tCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGa             <  1:6934794/93‑1 (MQ=255)
                                                                                tCTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGa             <  2:857242/93‑1 (MQ=255)
                                                                                 cTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGAt            <  2:2835644/93‑1 (MQ=255)
                                                                                 cTTCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGAt            >  2:1032377/1‑93 (MQ=255)
                                                                                  ttCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATc           >  2:3666378/1‑93 (MQ=255)
                                                                                  ttCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATc           <  1:5996129/93‑1 (MQ=255)
                                                                                  ttCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATc           >  1:4619025/1‑93 (MQ=255)
                                                                                  ttCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATc           >  1:4305104/1‑93 (MQ=255)
                                                                                   tCAAAGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCg          >  1:1340003/1‑93 (MQ=255)
                                                                                     aaaGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGc        >  2:4599232/1‑93 (MQ=255)
                                                                                     aaaGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGc        >  2:1012688/1‑93 (MQ=255)
                                                                                     aaaGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGc        <  1:3323264/93‑1 (MQ=255)
                                                                                      aaGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGGGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCa       <  1:4733865/93‑1 (MQ=255)
                                                                                      aaGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCa       <  2:5941815/93‑1 (MQ=255)
                                                                                       aGACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAg      >  1:5939398/1‑93 (MQ=255)
                                                                                        gACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGc     >  2:2187887/1‑93 (MQ=255)
                                                                                        gACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGc     >  1:7388743/1‑93 (MQ=255)
                                                                                        gACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGc     >  1:5008653/1‑93 (MQ=255)
                                                                                        gACCGCACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGc     <  1:4550361/93‑1 (MQ=255)
                                                                                            |                                                                                           
GATAACGCTTCACTCACAGCAAGCGAACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGAC  >  NC_000913/2814946‑2815129

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