Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,450,844 T→C T22A (ACC→GCC)  rpsJ ← 30S ribosomal subunit protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,450,8440TC90.8% 1192.2 / 86.4 457T22A (ACC→GCC) rpsJ30S ribosomal subunit protein S10
Reads supporting (aligned to +/- strand):  ref base T (12/30);  new base C (211/203);  total (223/234)
Fisher's exact test for biased strand distribution p-value = 5.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 481 total aligned reads displayed.

CAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAAC  >  NC_000913/3450752‑3450936
                                                                                            |                                                                                            
cAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGt                                                                                              <  1:6963657/93‑1 (MQ=255)
  gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTg                                                                                            >  2:6697247/1‑93 (MQ=255)
  gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTg                                                                                            >  2:5645714/1‑93 (MQ=255)
      aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTg                                                                                        >  2:9825075/1‑93 (MQ=255)
       gTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                       <  1:6450272/93‑1 (MQ=255)
       gTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTga                                                                                       <  1:8085682/93‑1 (MQ=255)
       gTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTga                                                                                       <  1:8202019/93‑1 (MQ=255)
        tGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgat                                                                                      >  2:312490/1‑93 (MQ=255)
        tGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGGCTCGACGATTTCCGCGGCtgattgat                                                                                      >  1:421377/1‑93 (MQ=255)
         gAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatc                                                                                     <  2:5928058/93‑1 (MQ=255)
         gAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatc                                                                                     <  2:1675259/93‑1 (MQ=255)
           aGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGGTGCTTgatcga                                                                                   <  2:7485396/93‑1 (MQ=255)
           aGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcga                                                                                   <  1:2521455/93‑1 (MQ=255)
           aGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcga                                                                                   <  1:2412191/93‑1 (MQ=255)
           aGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcga                                                                                   >  2:4519127/1‑93 (MQ=255)
            gcgcTCTGTGCGTGTCGGCATCGGGATCGGACCACGGACCTGCGCACGAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcgat                                                                                  <  2:920804/93‑1 (MQ=255)
             cgcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcgatc                                                                                 <  2:7045436/93‑1 (MQ=255)
             cgcTCTTTGCGTGTCGGCAGCGGCATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTgatcgatc                                                                                 <  2:7224077/93‑1 (MQ=255)
              gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCa                                                                                <  2:4300019/93‑1 (MQ=255)
              gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCa                                                                                <  1:2090496/93‑1 (MQ=255)
                tctTTGCGTGTAGGCAGCGGGATCGAACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGa                                                                              <  2:713593/93‑1 (MQ=255)
                 ctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGAc                                                                             >  2:540185/1‑93 (MQ=255)
                 ctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGAGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGAACGATCAGAc                                                                             >  2:9419833/1‑93 (MQ=255)
                  tttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACg                                                                            <  1:1524772/93‑1 (MQ=255)
                   ttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACga                                                                           <  2:3581266/93‑1 (MQ=255)
                    tGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACgat                                                                          <  1:3185386/93‑1 (MQ=255)
                    tGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACgat                                                                          <  2:2998426/93‑1 (MQ=255)
                      cGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACgatga                                                                        <  2:1892347/93‑1 (MQ=255)
                       gtgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACgataat                                                                       >  1:7331502/1‑93 (MQ=255)
                         gtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCa                                                                     >  2:9738816/1‑93 (MQ=255)
                         gtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCGGGTTAATGGATAAGACGATGATCa                                                                     >  2:9539387/1‑93 (MQ=255)
                           cGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCaaa                                                                   <  2:8423726/93‑1 (MQ=255)
                           cGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCaaa                                                                   <  2:5286371/93‑1 (MQ=255)
                            ggCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAAc                                                                  <  1:6014457/93‑1 (MQ=255)
                            ggCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAAc                                                                  <  2:1536298/93‑1 (MQ=255)
                               aGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCt                                                               <  2:825100/93‑1 (MQ=255)
                                gctggATCGGACCACGGACCTGCGCACCAGTGCGCGTGGCAGTCTCGACGATTTCCGAGGCTGCTTGATCGATCAGACGATGATCAAACGCtt                                                              <  1:10013723/90‑1 (MQ=255)
                                 cGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCttt                                                             >  1:8498199/1‑93 (MQ=255)
                                  gggATCGGACCGCGGACCTGCGCACCAGTGCGCTGGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTc                                                            <  2:2313650/93‑1 (MQ=255)
                                   ggATCGGACCACGGGCCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                           <  1:3482641/93‑1 (MQ=255)
                                   ggATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                           <  2:849295/93‑1 (MQ=255)
                                   ggATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                           >  2:6101097/1‑93 (MQ=255)
                                    gATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAg                                                          >  1:2788076/1‑93 (MQ=255)
                                    gATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTCCTTGATCGATCAGACGATGATCAAACGCTTTCAg                                                          <  1:3324116/93‑1 (MQ=255)
                                     aTCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                         <  2:4838976/93‑1 (MQ=255)
                                     aTCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                         <  1:562566/93‑1 (MQ=255)
                                     aTCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                         <  1:3605011/93‑1 (MQ=255)
                                      tCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                        >  2:5822484/1‑93 (MQ=255)
                                      tCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                        <  1:9972838/93‑1 (MQ=255)
                                      tCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                        <  2:3471425/93‑1 (MQ=255)
                                         gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                     <  2:9134466/93‑1 (MQ=255)
                                          aCCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGAt                                                    >  2:7389271/1‑93 (MQ=255)
                                           ccACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATa                                                   <  1:7574533/93‑1 (MQ=255)
                                           ccACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATa                                                   <  2:2789332/93‑1 (MQ=255)
                                           ccACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATa                                                   <  1:1964233/93‑1 (MQ=255)
                                             aCGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCTAACGCTTTCAGGCGGATACg                                                 >  2:1992346/1‑93 (MQ=255)
                                               ggACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                               <  1:9289683/93‑1 (MQ=255)
                                                 gccTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGAtt                                             <  2:2720702/92‑1 (MQ=255)
                                                  ccTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTc                                            <  1:1764140/93‑1 (MQ=255)
                                                    tGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCtt                                          <  2:8694987/93‑1 (MQ=255)
                                                    tGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCtt                                          >  2:1326161/1‑93 (MQ=255)
                                                    tGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCtt                                          <  2:5124665/93‑1 (MQ=255)
                                                       gcACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                       <  1:5255139/93‑1 (MQ=255)
                                                       gcACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                       <  1:4297487/93‑1 (MQ=255)
                                                       gcACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                       >  1:5777584/1‑93 (MQ=255)
                                                         aCCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                     >  1:1190260/1‑93 (MQ=255)
                                                         aCCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                     >  1:292013/1‑93 (MQ=255)
                                                         aCCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                     >  1:7810251/1‑93 (MQ=255)
                                                          ccAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTc                                    <  1:50659/93‑1 (MQ=255)
                                                           cAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCt                                   >  1:861671/1‑93 (MQ=255)
                                                             gTGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGc                                 >  1:5397842/1‑93 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                <  2:8871951/93‑1 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCCGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTTGGTCTGCa                                >  1:5477449/1‑93 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                >  1:4939494/1‑93 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                >  1:4563838/1‑93 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                >  1:9486626/1‑93 (MQ=255)
                                                              tGCGCTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                >  1:9091813/1‑93 (MQ=255)
                                                               gcgcTTGGCAGTCTCGGCGATTTCCGAGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGAGGATACGGATTCTTTGGATCTGCaa                               <  2:8992870/93‑2 (MQ=255)
                                                               gcgcTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCAt                               >  2:7905991/1‑93 (MQ=255)
                                                                  cTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgag                            <  2:2555558/93‑1 (MQ=255)
                                                                  cTTGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgag                            <  1:1326226/93‑1 (MQ=255)
                                                                   ttGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgaga                           <  1:3861877/93‑1 (MQ=255)
                                                                   ttGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgaga                           <  1:2965327/93‑1 (MQ=255)
                                                                   ttGGCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgaga                           >  2:2351751/1‑93 (MQ=255)
                                                                     ggCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAcc                         >  2:3946750/1‑93 (MQ=255)
                                                                      gCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                        <  1:613643/93‑1 (MQ=255)
                                                                      gCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                        >  2:816990/1‑93 (MQ=255)
                                                                      gCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                        <  1:9003921/93‑1 (MQ=255)
                                                                      gCAGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                        <  2:110755/93‑1 (MQ=255)
                                                                        aGTCTCGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCaga                      <  2:4066624/93‑1 (MQ=255)
                                                                          tctcGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGc                    >  1:6309042/1‑93 (MQ=255)
                                                                          tctcGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGc                    <  2:6311254/93‑1 (MQ=255)
                                                                            tcGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTc                  >  1:7989950/1‑93 (MQ=255)
                                                                            tcGACGATTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTc                  <  1:9749707/93‑1 (MQ=255)
                                                                            tcGACGATTTCCGCCGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGAGTCTTTGGTTCTGCATGAGAGCAGAGCTc                  >  2:3844319/1‑93 (MQ=255)
                                                                             cgacgaTTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTcc                 <  1:4815905/93‑1 (MQ=255)
                                                                               acgaTTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa               >  1:7481436/1‑93 (MQ=255)
                                                                                 gaTTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatt             <  1:5951426/93‑1 (MQ=255)
                                                                                  aTTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta            >  1:6750964/1‑93 (MQ=255)
                                                                                  aTTTCCGCGGCTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta            >  2:209332/1‑93 (MQ=255)
                                                                                            |                                                                                            
CAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAAC  >  NC_000913/3450752‑3450936

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