Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,595,908 T→C Q355R (CAG→CGG)  mdoB ← phosphoglycerol transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,595,9080TC85.3% 736.5 / 119.7 313Q355R (CAG→CGG) mdoBphosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (125/142);  total (148/165)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 314 total aligned reads displayed.

TCATCGCTAAATGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATCCGG  >  NC_000913/4595816‑4595999
                                                                                            |                                                                                           
 cATCGCTAAATGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTg                                                                                            <  1:548750/93‑1 (MQ=255)
   tCGCTAAATGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGc                                                                                          >  2:6287579/1‑93 (MQ=255)
   tCGCTAAATGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGc                                                                                          <  1:767963/93‑1 (MQ=255)
         aaTGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAAc                                                                                    <  1:7464982/93‑1 (MQ=255)
           tGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACtt                                                                                  <  1:5200687/93‑1 (MQ=255)
             gTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTAc                                                                                <  2:3980455/93‑1 (MQ=255)
             gTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTAc                                                                                <  1:5410504/93‑1 (MQ=255)
             gTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTAc                                                                                <  2:4222082/93‑1 (MQ=255)
              tCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTAcc                                                                               <  2:1414279/93‑1 (MQ=255)
              tCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTAcc                                                                               <  2:7430625/93‑1 (MQ=255)
                agaagaGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACcgc                                                                             <  1:2248799/93‑1 (MQ=255)
                 gaagaGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACcgcg                                                                            >  1:3036295/1‑93 (MQ=255)
                   agagACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCt                                                                          >  1:1523460/1‑93 (MQ=255)
                    gagaCGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTg                                                                         >  2:3068988/1‑93 (MQ=255)
                    gagaCGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTg                                                                         >  1:3115027/1‑93 (MQ=255)
                     agaCGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCTCGCTGa                                                                        <  2:5772086/93‑1 (MQ=255)
                      gacgacGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGaa                                                                       <  2:3603334/93‑1 (MQ=255)
                          acgaTGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGAc                                                                   >  1:4547463/1‑93 (MQ=255)
                           cgaTGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACt                                                                  <  2:6619903/93‑1 (MQ=255)
                            gaTGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTg                                                                 >  2:3561534/1‑93 (MQ=255)
                            gaTGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTg                                                                 >  1:6342014/1‑93 (MQ=255)
                             aTGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGa                                                                >  1:595212/1‑93 (MQ=255)
                              tGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGAt                                                               <  2:175344/93‑1 (MQ=255)
                              tGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGAt                                                               <  2:4414967/93‑1 (MQ=255)
                              tGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGAt                                                               >  2:3745375/1‑93 (MQ=255)
                              tGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGAt                                                               <  2:2440490/93‑1 (MQ=255)
                               gACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGAtt                                                              >  2:6709688/1‑93 (MQ=255)
                                aCGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTc                                                             <  2:414397/93‑1 (MQ=255)
                                aCGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTc                                                             <  2:6300790/93‑1 (MQ=255)
                                 cGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCg                                                            >  1:4064921/1‑93 (MQ=255)
                                  ggTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCgg                                                           <  2:210969/93‑1 (MQ=255)
                                   gTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGt                                                          >  1:3211965/1‑93 (MQ=255)
                                   gTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGt                                                          >  2:4744566/1‑93 (MQ=255)
                                    tATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGtt                                                         >  1:1277651/1‑93 (MQ=255)
                                    caTCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGtt                                                         <  1:2757677/92‑1 (MQ=255)
                                     aTCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGttt                                                        <  1:1040788/93‑1 (MQ=255)
                                     aTCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGttt                                                        >  1:2021229/1‑93 (MQ=255)
                                     aTCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGttt                                                        >  1:5881554/1‑93 (MQ=255)
                                      tCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTa                                                       >  1:1903015/1‑93 (MQ=255)
                                      tCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTa                                                       >  2:2539268/1‑93 (MQ=255)
                                      tCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTa                                                       <  2:3862318/93‑1 (MQ=255)
                                       cTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTAc                                                      <  1:1765378/93‑1 (MQ=255)
                                        tttAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTAcc                                                     <  1:7146518/93‑1 (MQ=255)
                                        tttAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTAcc                                                     <  1:3936428/93‑1 (MQ=255)
                                         ttAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCg                                                    >  1:6255512/1‑93 (MQ=255)
                                          tAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGt                                                   >  2:5541501/1‑93 (MQ=255)
                                          tAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGt                                                   >  1:5587137/1‑93 (MQ=255)
                                            aaCCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCa                                                 <  1:3881632/93‑1 (MQ=255)
                                              ccACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCaaa                                               <  1:5105152/93‑1 (MQ=255)
                                                  ggTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCa                                           <  2:1215206/93‑1 (MQ=255)
                                                     gACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCatat                                        <  1:1348416/93‑1 (MQ=255)
                                                      aCGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAtt                                       <  1:2878807/93‑1 (MQ=255)
                                                       cGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttt                                      <  1:2627462/93‑1 (MQ=255)
                                                       cGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttt                                      >  2:5632797/1‑93 (MQ=255)
                                                        gCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAtttt                                     <  1:4324972/93‑1 (MQ=255)
                                                         cTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttttt                                    <  1:359048/93‑1 (MQ=255)
                                                         cTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttttt                                    >  2:651767/1‑93 (MQ=255)
                                                         cTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttttt                                    >  2:863496/1‑93 (MQ=255)
                                                         cTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAttttt                                    >  1:1255451/1‑93 (MQ=255)
                                                          tttGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAtttttt                                   <  2:5442180/93‑1 (MQ=255)
                                                          tttGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAtttttt                                   >  2:5981210/1‑93 (MQ=255)
                                                          tttGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATAtttttt                                   >  1:921782/1‑93 (MQ=255)
                                                           ttGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTg                                  >  2:5922335/1‑93 (MQ=255)
                                                           ttGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTg                                  <  2:2622139/93‑1 (MQ=255)
                                                           ttGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTg                                  <  1:2488002/93‑1 (MQ=255)
                                                            tGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGc                                 <  2:2858374/93‑1 (MQ=255)
                                                            tGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGc                                 <  2:5086765/93‑1 (MQ=255)
                                                            tGATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGc                                 <  1:7034236/93‑1 (MQ=255)
                                                             gATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCg                                <  1:6521385/93‑1 (MQ=255)
                                                             gATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCg                                <  1:6285469/93‑1 (MQ=255)
                                                             gATTTTGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCg                                <  1:3373241/93‑1 (MQ=255)
                                                               ttttGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGt                              >  2:3266169/1‑93 (MQ=255)
                                                               ttttGTTGATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGt                              >  1:3799992/1‑93 (MQ=255)
                                                                 ttgttgATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTa                            <  2:7316369/93‑1 (MQ=255)
                                                                  tgttgATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTAc                           <  2:1932190/93‑1 (MQ=255)
                                                                  tgttgATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTAc                           <  1:6618592/93‑1 (MQ=255)
                                                                   gttgATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACa                          >  1:1609571/1‑93 (MQ=255)
                                                                    ttgATAAACGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAg                         <  2:3296434/93‑1 (MQ=255)
                                                                     tgATAAACGTCGCGATGTTCTCCCGGCTGCGACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAgg                        <  2:6112537/93‑1 (MQ=255)
                                                                          aaCGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACga                   <  2:1019171/93‑1 (MQ=255)
                                                                          aaCGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACga                   <  2:2885705/93‑1 (MQ=255)
                                                                           aCGTCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACgag                  <  2:7125187/93‑1 (MQ=255)
                                                                             gTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACgagag                >  2:5152435/1‑93 (MQ=255)
                                                                             gTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACgagag                >  2:6991028/1‑93 (MQ=255)
                                                                              tCGCGATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACgagaga               >  1:4832347/1‑93 (MQ=255)
                                                                               cgcgATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGAt              <  1:6104528/93‑1 (MQ=255)
                                                                                gcgATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATa             >  1:527975/1‑93 (MQ=255)
                                                                                gcgATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATa             >  2:6337705/1‑93 (MQ=255)
                                                                                gcgATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATa             >  1:4203417/1‑93 (MQ=255)
                                                                                 cgATGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATaa            <  1:6371125/93‑1 (MQ=255)
                                                                                  gATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATaaa           >  1:68693/1‑93 (MQ=255)
                                                                                   aTGTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATaaaa          >  1:353988/1‑93 (MQ=255)
                                                                                     gTTCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAAcc        <  1:2053308/93‑1 (MQ=255)
                                                                                      ttCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCa       >  2:3765741/1‑93 (MQ=255)
                                                                                      ttCTCCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCa       >  1:4704889/1‑93 (MQ=255)
                                                                                        ctcCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATc     <  2:719782/93‑1 (MQ=255)
                                                                                         tcCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATcc    <  2:7242057/93‑1 (MQ=255)
                                                                                         tcCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATcc    >  2:6181614/1‑93 (MQ=255)
                                                                                         tcCCGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATcc    >  1:1477116/1‑93 (MQ=255)
                                                                                           cTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATCCgg  >  1:2374393/1‑93 (MQ=255)
                                                                                            |                                                                                           
TCATCGCTAAATGGTCAGAAGAGACGACGATGACGGTATCTTTAAACCACGGTGACGCTTTGATTTTGTTGATAAACGTCGCGATGTTCTCCTGGCTGCAACTTACCGCGCTGAATGACTGATTCGGTTTACCGTCAAAATCATATTTTTTGCGGTTACAGGTACGAGAGATAAAACCATCCGG  >  NC_000913/4595816‑4595999

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