Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,016,608 G→A G305G (GGC→GGT)  b0955 ← putative ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,016,6080GA96.6% 477.8 / 0.9 147G305G (GGC→GGT) b0955putative ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (72/70);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 152 total aligned reads displayed.

GATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCACACAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGTGCT  >  NC_000913/1016518‑1016698
                                                                                          |                                                                                          
gATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCAcc                                                                                          >  1:2795028/1‑93 (MQ=255)
gATCTCATCCTGCATACGTTCAGCGAGGAAAACTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGCAGGTGTCAcc                                                                                          >  1:3510113/1‑93 (MQ=255)
 aTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCa                                                                                         <  1:4136604/93‑1 (MQ=255)
 aTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCa                                                                                         <  2:3532240/93‑1 (MQ=255)
  tctcATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCAt                                                                                        >  1:555904/1‑93 (MQ=255)
  tctcATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCAt                                                                                        >  1:1700022/1‑93 (MQ=255)
    tcATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATca                                                                                      >  1:2421259/1‑93 (MQ=255)
      aTCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATcacg                                                                                    <  1:756655/93‑2 (MQ=255)
      aTCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATcacg                                                                                    >  2:1948872/1‑92 (MQ=255)
       tCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATcacgc                                                                                   >  2:2113648/1‑91 (MQ=255)
        ccTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATcacgca                                                                                  <  1:2976641/93‑1 (MQ=255)
         cTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAg                                                                                 >  1:2924722/1‑93 (MQ=255)
          tGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAgg                                                                                >  1:3404702/1‑93 (MQ=255)
            cATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAgggg                                                                              >  1:3799685/1‑92 (MQ=255)
            cATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGa                                                                              >  1:390397/1‑93 (MQ=255)
            cATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGa                                                                              >  1:291766/1‑93 (MQ=255)
            cATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGa                                                                              <  1:3435669/93‑1 (MQ=255)
            cATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGa                                                                              <  2:3979804/93‑1 (MQ=255)
              tACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAgg                                                                            <  2:2323989/93‑1 (MQ=255)
                cGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCg                                                                          <  2:372348/93‑1 (MQ=255)
                cGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGAGTCACCATCACGCAGGGAGGCg                                                                          >  2:4657526/1‑93 (MQ=255)
                     aGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACctc                                                                     >  2:3598252/1‑93 (MQ=255)
                      gCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACctct                                                                    >  2:4499746/1‑93 (MQ=255)
                       cGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCtt                                                                   >  2:1263338/1‑93 (MQ=255)
                        gAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCttt                                                                  >  1:1011127/1‑93 (MQ=255)
                        gAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCttt                                                                  >  1:1126655/1‑93 (MQ=255)
                         aGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTa                                                                 >  1:3569082/1‑93 (MQ=255)
                         aGGAAAACTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGGCACCATCACGCAGGGAGGGGAGctcactg                                                                 >  2:7978/1‑89 (MQ=255)
                          ggAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTac                                                                >  1:3600876/1‑93 (MQ=255)
                             aaCCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACt                                                             <  1:3307761/93‑1 (MQ=255)
                               ccTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTgg                                                           >  1:880029/1‑93 (MQ=255)
                                cTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGTAGGCGACCTCTTTACACTGGc                                                          <  1:4648890/93‑1 (MQ=255)
                                cTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGc                                                          <  1:2212018/93‑1 (MQ=255)
                                cTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGc                                                          <  1:1369932/93‑1 (MQ=255)
                                 ttCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCg                                                         <  1:4511918/93‑1 (MQ=255)
                                 ttCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCg                                                         >  2:923552/1‑93 (MQ=255)
                                  tCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCgg                                                        >  2:55430/1‑93 (MQ=255)
                                  tCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCgg                                                        <  1:1245342/93‑1 (MQ=255)
                                    gcgcCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCACGCAGGGAGGCGACCTCTTTACACTGGCggag                                                      <  2:3171040/93‑1 (MQ=255)
                                     cgcCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCggagg                                                     >  2:2957020/1‑93 (MQ=255)
                                       ccATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGAt                                                   <  1:1000465/93‑1 (MQ=255)
                                        cATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATc                                                  >  2:1544881/1‑93 (MQ=255)
                                         aTTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATcc                                                 <  2:780923/93‑1 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATccc                                                >  1:2808745/1‑92 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCa                                                >  2:4105916/1‑93 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCa                                                >  2:3286835/1‑93 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCa                                                >  1:1089387/1‑93 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCa                                                <  2:2734086/93‑1 (MQ=255)
                                          ttCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCa                                                >  1:2581984/1‑93 (MQ=255)
                                           tCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCAc                                               >  2:417511/1‑93 (MQ=255)
                                           tCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCAc                                               <  2:3844555/93‑1 (MQ=255)
                                           tCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCAc                                               <  2:3739935/93‑1 (MQ=255)
                                           tCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCAGCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCAc                                               <  2:259364/93‑1 (MQ=255)
                                            cACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACt                                              <  2:3051220/93‑1 (MQ=255)
                                             aCGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTg                                             >  1:4029957/1‑93 (MQ=255)
                                              cGCTACTGCAGCATTAAGTTTAGCTGCTCGCCGCAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGTATCCACTGc                                            <  2:172668/93‑1 (MQ=255)
                                               gctgctGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCg                                           <  2:2568304/93‑1 (MQ=255)
                                                ctgctgCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCgg                                          >  1:3776131/1‑93 (MQ=255)
                                                ctgctgCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCgg                                          >  1:1218745/1‑93 (MQ=255)
                                                ctgctgCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCgg                                          <  1:4658876/93‑1 (MQ=255)
                                                 tgctgcAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCggg                                         <  1:3850570/93‑1 (MQ=255)
                                                  gctgcAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGc                                        >  1:2003031/1‑93 (MQ=255)
                                                       agcaTTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGc                                   >  1:1738872/1‑93 (MQ=255)
                                                        gcaTTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCg                                  <  1:3249180/93‑1 (MQ=255)
                                                        gcaTTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCg                                  >  1:2296908/1‑93 (MQ=255)
                                                           ttAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGaa                               <  1:4569667/93‑1 (MQ=255)
                                                           ttAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGaa                               >  2:1237418/1‑93 (MQ=255)
                                                            tAAGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAg                              <  1:1507380/93‑1 (MQ=255)
                                                             aaGTTTATCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCTGAGGATCCACTGCGGGCTAATCGGAAgt                             <  2:2765294/93‑1 (MQ=255)
                                                             aaGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAgg                             >  1:3146607/1‑92 (MQ=255)
                                                              aGTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAgtg                            <  1:648474/93‑1 (MQ=255)
                                                               gTTTAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAgtgt                           >  2:1626982/1‑93 (MQ=255)
                                                                  tAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTc                        <  1:4387594/93‑1 (MQ=255)
                                                                  tAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTc                        >  1:2950543/1‑93 (MQ=255)
                                                                  tAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTc                        <  1:1421631/93‑1 (MQ=255)
                                                                  tAGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTc                        <  2:213924/93‑1 (MQ=255)
                                                                   aGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCt                       <  2:685076/93‑1 (MQ=255)
                                                                   aGCTGCTCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCt                       >  1:2202364/1‑93 (MQ=255)
                                                                    gctgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCtt                      <  2:2275985/93‑1 (MQ=255)
                                                                    gctgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCtt                      >  2:3467113/1‑93 (MQ=255)
                                                                    gctgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCtt                      <  2:2235306/93‑1 (MQ=255)
                                                                     ctgctCGCCGGAGAAGGTGTCACCATCACGCAGTGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTg                     <  2:2028811/93‑1 (MQ=255)
                                                                     ctgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTg                     <  2:2499163/93‑1 (MQ=255)
                                                                     ctgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTg                     >  2:1511119/1‑93 (MQ=255)
                                                                      tgctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGt                    >  2:348369/1‑93 (MQ=255)
                                                                       gctCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGGGTTTCTTGtt                   >  2:546856/1‑93 (MQ=255)
                                                                         tCGCCGGAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTGTCTTGttga                 >  2:1688788/1‑91 (MQ=255)
                                                                              ggAGAAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGgt            <  2:1707379/93‑1 (MQ=255)
                                                                               gagaAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGgtg           >  1:4036523/1‑93 (MQ=255)
                                                                                agaAGGTGTCGCCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGgtgt          <  2:2269349/93‑1 (MQ=255)
                                                                                 gaAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGgtgca         >  2:4461162/1‑91 (MQ=255)
                                                                                 gaAGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTa         <  2:4614974/93‑1 (MQ=255)
                                                                                  gaGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAg        <  1:843463/92‑1 (MQ=255)
                                                                                  aaGGTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAg        <  2:2451199/93‑1 (MQ=255)
                                                                                    ggTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCg      >  2:3319287/1‑93 (MQ=255)
                                                                                    ggTGTCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCg      >  1:2289576/1‑93 (MQ=255)
                                                                                     gtgtCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGt     >  1:3677726/1‑93 (MQ=255)
                                                                                     gtgtCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGt     <  1:4324335/93‑1 (MQ=255)
                                                                                      tgtCACCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGtg    >  2:678540/1‑93 (MQ=255)
                                                                                        tCGCCATCACGCAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGtgct  <  1:4421164/93‑1 (MQ=255)
                                                                                          |                                                                                          
GATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCACACAGGGAGGCGACCTCTTTACACTGGCGGAGGATCCACTGCGGGCTAAGCGGAAGTGTTTCTTGTTCACCGGTGTAGCGTGCT  >  NC_000913/1016518‑1016698

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