Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,456,098 A→G Q41Q (CAA→CAG yheG → YheG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,456,0980AG99.6% 845.0 / NA 262Q41Q (CAA→CAGyheGYheG
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (136/124);  total (138/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

GCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456009‑3456187
                                                                                         |                                                                                         
  gcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCa                                                                                      <  1:1381235/93‑1 (MQ=255)
   cTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAg                                                                                     <  2:2741575/93‑1 (MQ=255)
   cTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAg                                                                                     >  2:3093458/1‑93 (MQ=255)
   cTGGCACAACGAATTGGTTAACTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAg                                                                                     >  2:3059888/1‑93 (MQ=255)
    tGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAgg                                                                                    >  1:95296/1‑93 (MQ=255)
     ggCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGt                                                                                   >  1:4607467/1‑93 (MQ=255)
     ggCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGt                                                                                   <  2:1394017/93‑1 (MQ=255)
      gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGtt                                                                                  <  1:5887410/93‑1 (MQ=255)
      gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGtt                                                                                  <  2:6212038/93‑1 (MQ=255)
       cacaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTg                                                                                 >  2:2458121/1‑93 (MQ=255)
       cacaACGAATTGGTTAGCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTg                                                                                 <  2:5345607/93‑1 (MQ=255)
          aaCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGcc                                                                              <  2:4882573/93‑1 (MQ=255)
          aaCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGcc                                                                              <  1:4445875/93‑1 (MQ=255)
           aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             <  2:411894/93‑1 (MQ=255)
             gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGt                                                                           <  1:1727477/93‑1 (MQ=255)
              aaTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTg                                                                          >  1:3998178/1‑93 (MQ=255)
              aaTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTg                                                                          >  2:2473804/1‑93 (MQ=255)
               aTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGc                                                                         <  1:3821449/93‑1 (MQ=255)
                ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGcc                                                                        >  1:1676756/1‑93 (MQ=255)
                   gTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGa                                                                     >  1:3152100/1‑93 (MQ=255)
                   gTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGa                                                                     <  1:1297808/93‑1 (MQ=255)
                    ttATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGAc                                                                    <  1:4781895/93‑1 (MQ=255)
                    ttATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGAc                                                                    <  1:4912388/93‑1 (MQ=255)
                      aTCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCa                                                                  >  1:1151417/1‑93 (MQ=255)
                      aTCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCa                                                                  <  2:3533963/93‑1 (MQ=255)
                      aTCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCa                                                                  <  2:4975534/93‑1 (MQ=255)
                       tctcTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCAt                                                                 <  1:2322123/93‑1 (MQ=255)
                         tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGa                                                               >  1:58339/1‑93 (MQ=255)
                         tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGa                                                               <  1:5988376/93‑1 (MQ=255)
                         tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGa                                                               <  1:5083032/93‑1 (MQ=255)
                          ctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAc                                                              <  2:4143572/93‑1 (MQ=255)
                           tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             <  2:706941/93‑1 (MQ=255)
                           tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             <  1:1920060/93‑1 (MQ=255)
                            aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            <  2:2356546/93‑1 (MQ=255)
                            aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            >  2:2927533/1‑93 (MQ=255)
                              tCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                          >  1:2679869/1‑93 (MQ=255)
                               cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCt                                                         <  2:2611186/93‑1 (MQ=255)
                                agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        <  2:1719598/93‑1 (MQ=255)
                                agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        >  1:5519880/1‑93 (MQ=255)
                                agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        <  2:3494008/93‑1 (MQ=255)
                                   gggCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGGTACTCGATATTCAGCGTCAGGTTAGCCAGTGCCAGACCATGACCAGCTTGTc                                                     >  2:1207585/1‑93 (MQ=255)
                                     gCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcac                                                   >  1:6277324/1‑93 (MQ=255)
                                       gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAg                                                 >  2:5255891/1‑93 (MQ=255)
                                       gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAg                                                 >  2:4258063/1‑93 (MQ=255)
                                       gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAg                                                 <  1:4664460/93‑1 (MQ=255)
                                        agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAgg                                                >  2:740920/1‑93 (MQ=255)
                                         gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                               >  1:4348846/1‑93 (MQ=255)
                                         gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                               <  1:5028687/93‑1 (MQ=255)
                                          attattTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                              <  1:348341/93‑1 (MQ=255)
                                            tatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGa                                            >  2:6097052/1‑93 (MQ=255)
                                              tttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtt                                          >  1:3118169/1‑93 (MQ=255)
                                               ttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAttt                                         >  2:1514859/1‑93 (MQ=255)
                                                  tctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      <  1:2253178/93‑1 (MQ=255)
                                                  tctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      >  1:3745359/1‑93 (MQ=255)
                                                   ctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                     <  1:5739512/93‑1 (MQ=255)
                                                   ctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                     <  2:3715303/93‑1 (MQ=255)
                                                   ctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                     <  2:5568958/93‑1 (MQ=255)
                                                   ctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                     <  1:6151338/93‑1 (MQ=255)
                                                    tcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAg                                    >  2:1976282/1‑93 (MQ=255)
                                                     cGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGc                                   <  2:5500421/93‑1 (MQ=255)
                                                      gCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCt                                  >  1:4182946/1‑93 (MQ=255)
                                                         gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               >  1:3659343/1‑93 (MQ=255)
                                                          aTACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGc                              >  2:4032683/1‑93 (MQ=255)
                                                            aCGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTg                            >  1:6202474/1‑93 (MQ=255)
                                                             cGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTgg                           <  1:5414811/93‑1 (MQ=255)
                                                             cGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTgg                           >  1:3269458/1‑93 (MQ=255)
                                                              gCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGa                          <  2:1126012/93‑1 (MQ=255)
                                                               ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         <  2:4585840/93‑1 (MQ=255)
                                                               ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCCGCTTGaa                         >  1:713898/1‑93 (MQ=255)
                                                                cATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAg                        >  1:5259061/1‑93 (MQ=255)
                                                                 aTTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGGAGc                       >  2:5921375/1‑93 (MQ=255)
                                                                 aTTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGc                       <  1:5995574/93‑1 (MQ=255)
                                                                 aTTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGc                       >  2:3208365/1‑93 (MQ=255)
                                                                   tCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                     >  1:1812129/1‑93 (MQ=255)
                                                                    cGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGt                    <  2:5298943/93‑1 (MQ=255)
                                                                     gAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTTAAGCTCGGCTGGGAGCggga                   >  2:3327686/1‑91 (MQ=255)
                                                                     gAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                   <  1:1487112/93‑1 (MQ=255)
                                                                      aGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAt                  >  2:1467617/1‑93 (MQ=255)
                                                                      aGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAt                  >  2:658293/1‑93 (MQ=255)
                                                                          gTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCa              >  2:393459/1‑93 (MQ=255)
                                                                          gTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCa              <  1:5691025/93‑1 (MQ=255)
                                                                          gTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCa              >  2:5982619/1‑93 (MQ=255)
                                                                           ttACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCAt             <  2:874766/93‑1 (MQ=255)
                                                                            tACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc            >  2:4490423/1‑93 (MQ=255)
                                                                            tACTCGATATTCACCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGCTTTTGAAGCTCAGCTGGGAGCGGGATTCCATc            >  2:467649/1‑93 (MQ=255)
                                                                             aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           <  1:5773373/93‑1 (MQ=255)
                                                                             aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           <  2:1884385/93‑1 (MQ=255)
                                                                              cTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAg          <  2:5691242/93‑1 (MQ=255)
                                                                               tCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGa         >  1:1997100/1‑93 (MQ=255)
                                                                               tCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGGATTCCATCAGa         >  2:5774755/1‑93 (MQ=255)
                                                                                 gATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt       <  1:3059028/93‑1 (MQ=255)
                                                                                 gATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt       >  1:1096917/1‑93 (MQ=255)
                                                                                  atatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATa      <  1:2849738/93‑1 (MQ=255)
                                                                                   tatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAc     <  2:3980432/93‑1 (MQ=255)
                                                                                   tatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAc     <  2:3287567/93‑1 (MQ=255)
                                                                                   tatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAc     <  1:1236350/93‑1 (MQ=255)
                                                                                   tatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGATGATTTTGAAGCTCGGCTGGAAGCGGGATTCCATCAGAATAc     >  1:3242474/1‑93 (MQ=255)
                                                                                     ttCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCg   <  2:19859/93‑1 (MQ=255)
                                                                                      tCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg  <  2:3416649/93‑1 (MQ=255)
                                                                                      tCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg  >  2:1689212/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456009‑3456187

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