Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,078,857 T→C Y143C (TAC→TGC)  sbmC ← SbmC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,078,8570TC99.4% 553.5 / NA 174Y143C (TAC→TGC) sbmCSbmC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (77/96);  total (78/96)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

AAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAG  >  NC_000913/2078767‑2078946
                                                                                          |                                                                                         
aaaaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCAc                                                                                         >  1:2516761/1‑93 (MQ=255)
    aaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCAt                                                                                     >  2:2688908/1‑93 (MQ=255)
      agagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCt                                                                                   <  1:887907/93‑1 (MQ=255)
      agagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCt                                                                                   >  1:1532309/1‑93 (MQ=255)
      agagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCt                                                                                   >  1:3092129/1‑93 (MQ=255)
        agTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTc                                                                                 <  1:4336837/93‑1 (MQ=255)
        agTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTc                                                                                 <  2:2606606/93‑1 (MQ=255)
         gTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGTTATCCCAGCACCCATCTTcc                                                                                >  1:3523291/1‑93 (MQ=255)
         gTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTcc                                                                                >  1:1134373/1‑93 (MQ=255)
         gTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTcc                                                                                <  2:353636/93‑1 (MQ=255)
         gTGGGCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCGCCCATCTTcc                                                                                <  1:899693/93‑1 (MQ=255)
          tGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCcg                                                                               >  1:4713037/1‑93 (MQ=255)
              cATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGcc                                                                           <  1:2447884/93‑1 (MQ=255)
               aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCa                                                                          <  2:377403/93‑1 (MQ=255)
               aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCa                                                                          >  2:578154/1‑93 (MQ=255)
               aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCa                                                                          >  1:367100/1‑93 (MQ=255)
                tCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAt                                                                         <  2:2727031/93‑1 (MQ=255)
                tCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAt                                                                         <  2:1773587/93‑1 (MQ=255)
                tCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAt                                                                         >  2:2491812/1‑93 (MQ=255)
                tCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAt                                                                         <  1:955120/93‑1 (MQ=255)
                  gcgTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTg                                                                       <  2:2104180/93‑1 (MQ=255)
                   cgTTAACACACCGCCCTGAGATGGATTAGTGATGGTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGt                                                                      <  1:2501314/93‑1 (MQ=255)
                   cgTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGt                                                                      <  2:804307/93‑1 (MQ=255)
                     ttAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTc                                                                    <  2:1966112/93‑1 (MQ=255)
                     ttAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTc                                                                    >  1:640539/1‑93 (MQ=255)
                     ttAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTc                                                                    <  2:1227830/93‑1 (MQ=255)
                     ttAACACACCGCACTAAGATGAATTAGTGATGTTTTGGCGGCACCGCAACATCCATTTCGATATCCCAGCGCCCATCTTCCGCGCCATTGTTc                                                                    <  2:2667247/93‑1 (MQ=255)
                     ttAACACACCCCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTc                                                                    <  2:2667246/93‑1 (MQ=255)
                      tAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCa                                                                   <  1:3707957/93‑1 (MQ=255)
                       aaCACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCaa                                                                  <  2:2031476/93‑1 (MQ=255)
                          acacCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATa                                                               >  2:1692853/1‑93 (MQ=255)
                            acCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATaaa                                                             >  2:516775/1‑93 (MQ=255)
                            acCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATaaa                                                             >  1:4112191/1‑93 (MQ=255)
                             ccGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                            >  2:981181/1‑93 (MQ=255)
                              cGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAcc                                                           <  1:337269/93‑1 (MQ=255)
                              cGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAcc                                                           <  2:3401178/93‑1 (MQ=255)
                               gCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCt                                                          <  2:2237196/93‑1 (MQ=255)
                               gCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCt                                                          <  2:413450/93‑1 (MQ=255)
                               cgccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCt                                                          <  2:1932454/91‑1 (MQ=255)
                                    gagaTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGacg                                                     >  1:1614325/1‑91 (MQ=255)
                                     agaTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGc                                                    >  2:1218448/1‑93 (MQ=255)
                                     agaTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGc                                                    <  1:148261/93‑1 (MQ=255)
                                      gaTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCa                                                   <  2:79594/93‑1 (MQ=255)
                                       aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  >  1:1746573/1‑93 (MQ=255)
                                       aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  2:3067018/93‑1 (MQ=255)
                                         gAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAggg                                                <  2:2473654/93‑1 (MQ=255)
                                          aaTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGc                                               >  1:2282108/1‑93 (MQ=255)
                                          aaTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGc                                               <  1:3121894/93‑1 (MQ=255)
                                          aaTTAGTGATGTTTTGGCTGAACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATCGTTCAAATAAACCTCGAAGCAggga                                               >  2:707601/1‑92 (MQ=255)
                                           aTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCt                                              >  1:3830449/1‑93 (MQ=255)
                                           aTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCt                                              >  2:1178624/1‑93 (MQ=255)
                                           aTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCt                                              >  2:3316321/1‑93 (MQ=255)
                                             tAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCttt                                            <  2:4269581/93‑1 (MQ=255)
                                              aGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTg                                           <  1:4563912/93‑1 (MQ=255)
                                                tGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGt                                         >  2:755572/1‑93 (MQ=255)
                                                tGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGt                                         >  2:4369096/1‑93 (MQ=255)
                                                 gATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTa                                        >  1:1941315/1‑93 (MQ=255)
                                                  aTGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTaa                                       <  2:4046401/93‑1 (MQ=255)
                                                   tGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                      <  1:2403222/93‑1 (MQ=255)
                                                    gTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAACCCTCGAAGCAGGGCTTTGGTAACa                                     <  2:3186421/93‑1 (MQ=255)
                                                        tGGCTGCCCCGCAACATTCATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTc                                 <  1:4038034/93‑1 (MQ=255)
                                                         ggCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCa                                <  2:1626330/93‑1 (MQ=255)
                                                          gCTGCCCCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCAt                               <  1:1303621/93‑1 (MQ=255)
                                                          gCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCAt                               >  2:3919000/1‑93 (MQ=255)
                                                          gCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCAt                               <  1:367037/93‑1 (MQ=255)
                                                            tGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATaa                             >  2:1732962/1‑93 (MQ=255)
                                                            tGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATaa                             >  2:120153/1‑93 (MQ=255)
                                                             gCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAg                            >  2:3090583/1‑93 (MQ=255)
                                                                  gCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTg                       >  2:3091916/1‑93 (MQ=255)
                                                                  gCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTg                       >  2:2952881/1‑93 (MQ=255)
                                                                   cAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGt                      >  1:1669227/1‑93 (MQ=255)
                                                                   cAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGt                      <  2:1101395/93‑1 (MQ=255)
                                                                   cAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATCTCATAAGCACTGt                      >  2:1863179/1‑93 (MQ=255)
                                                                    aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     <  2:1519838/93‑1 (MQ=255)
                                                                    aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     >  2:3502358/1‑93 (MQ=255)
                                                                     acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    >  1:712419/1‑93 (MQ=255)
                                                                     acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    <  2:1126698/93‑1 (MQ=255)
                                                                     acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    <  2:3890327/93‑1 (MQ=255)
                                                                     acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    <  2:1891105/93‑1 (MQ=255)
                                                                      catacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCt                   >  2:4733842/1‑93 (MQ=255)
                                                                      catacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCt                   >  1:1691139/1‑93 (MQ=255)
                                                                      catacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCt                   <  1:1524379/93‑1 (MQ=255)
                                                                       atacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTg                  <  1:2984123/93‑1 (MQ=255)
                                                                        tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 >  1:593424/1‑93 (MQ=255)
                                                                        tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 <  2:4531200/93‑1 (MQ=255)
                                                                        tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGCTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 >  2:908391/1‑93 (MQ=255)
                                                                          catTTCTATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCaa               <  1:1075270/93‑1 (MQ=255)
                                                                          catTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCaa               <  2:1546864/93‑1 (MQ=255)
                                                                            tttCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga             <  2:279014/93‑1 (MQ=255)
                                                                            tttCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga             <  2:279013/93‑1 (MQ=255)
                                                                            tttCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga             >  1:2789583/1‑93 (MQ=255)
                                                                            tttCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga             <  1:820301/93‑1 (MQ=255)
                                                                              tCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAgaga           <  2:2037538/93‑1 (MQ=255)
                                                                              tCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAgaga           <  1:4470950/93‑1 (MQ=255)
                                                                               cGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGAc          <  2:2460785/93‑1 (MQ=255)
                                                                               cGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGAc          >  1:2811866/1‑93 (MQ=255)
                                                                                gATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACt         <  1:2667675/93‑1 (MQ=255)
                                                                                     cccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaa    <  1:3891216/93‑1 (MQ=255)
                                                                                     cccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaa    <  1:394179/93‑1 (MQ=255)
                                                                                     cccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaa    >  1:62792/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAG  >  NC_000913/2078767‑2078946

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