Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,124 +G intergenic (+58/+26) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,1231.G95.1% 226.6 / 7.3 103intergenic (+57/+27)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base G (44/53);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGT  >  NC_000913/3957050‑3957194
                                                                                |                                                                        
ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                           >  2:2183171/1‑76 (MQ=17)
ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                           >  1:4098605/1‑76 (MQ=17)
  gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGTGGGAGAAAACACTCAAgg                                                                        >  1:667915/1‑75 (MQ=17)
  gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                        >  2:3331788/1‑75 (MQ=17)
  gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                        >  1:2672247/1‑75 (MQ=17)
   cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                                       <  2:4829245/76‑3 (MQ=17)
   cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                                       >  1:2894758/1‑74 (MQ=17)
    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                                      >  2:4053990/1‑73 (MQ=17)
    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                                      <  1:2029390/76‑4 (MQ=17)
    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                                      >  2:3383084/1‑73 (MQ=17)
    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                                      <  2:2672247/76‑4 (MQ=17)
       tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                   <  2:5180557/76‑1 (MQ=14)
       tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                   <  2:773407/76‑1 (MQ=14)
        gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                  >  1:3556923/1‑76 (MQ=14)
        gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                  <  1:3918899/76‑1 (MQ=14)
        gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                  <  1:1389653/76‑1 (MQ=14)
         cGGGTTAAGTGCGCGCTGATGGCCCCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGGAAACACACAAAGCCGtctc                                                                 >  1:3648481/1‑76 (MQ=2)
          gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCCCCGGGGGCGGGAGGAAACACTCAAAGCCCTCTcc                                                                >  1:2448633/1‑76 (MQ=9)
          gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                >  1:3834939/1‑76 (MQ=14)
           ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                               <  1:494626/76‑1 (MQ=14)
           ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                               <  2:5838529/76‑1 (MQ=14)
              tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgta                                                            <  1:2826526/76‑3 (MQ=14)
              tAAGGGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACAGGGAGAGGGGGACAACACTCACAGCGTTCTCCTggc                                                            >  2:63938/1‑75 (MQ=2)
               aattGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                                                           <  2:2108225/73‑1 (MQ=14)
               aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGaag                                                           >  1:5059748/1‑76 (MQ=14)
                aGTGCGCGCTGATGCCCTCACCCCG‑‑ACGCTGTCGCGCGGGGGGAGGGAGAAAACACTCAGGGCCTTCTCCTGgaga                                                          <  2:453276/76‑1 (MQ=9)
                aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                                          <  2:184316/76‑1 (MQ=14)
                 gTGCGCGCTGCTGCCCTCACCCCG‑‑ACCCGCTCCCGCAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                         <  2:1926414/76‑1 (MQ=9)
                 gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                         <  1:740381/76‑1 (MQ=14)
                 gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                         <  2:5297093/76‑1 (MQ=14)
                 gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                         <  2:5185756/76‑1 (MQ=14)
                  tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgagacg                                                        >  2:1598999/1‑74 (MQ=14)
                  tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                        <  2:4408198/76‑1 (MQ=17)
                  tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                        <  2:2604769/76‑1 (MQ=14)
                   gcgcgcTGATGCCCTCACCCGT‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                       <  1:2184182/76‑1 (MQ=17)
                   gcgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                       >  2:3241600/1‑76 (MQ=17)
                    cgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                      <  1:5745061/76‑1 (MQ=17)
                     gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                     >  1:1500671/1‑76 (MQ=17)
                     gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                     >  1:4704333/1‑76 (MQ=17)
                     gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                     >  1:419565/1‑76 (MQ=17)
                     gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                     >  1:3271886/1‑76 (MQ=17)
                     gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCCTCTCCTGGAGAAGGcc                                                     >  1:3257383/1‑76 (MQ=17)
                      cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                    <  2:4455165/76‑1 (MQ=17)
                      cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                    >  2:4567663/1‑76 (MQ=17)
                       gcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCCTCTCCTGGCGAAGGccct                                                   >  2:2552597/1‑74 (MQ=14)
                       gcTGATGCCCTCACCCAG‑‑ACCCTCTCCCACACGGAGAGGGAGAAAACACACAAGGCCTTCTCCGCCAGAAGGCCta                                                   >  2:2485822/1‑75 (MQ=2)
                        cTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                  <  2:5258299/76‑1 (MQ=18)
                        cTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                  <  2:4338484/76‑1 (MQ=18)
                         tGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                 <  2:4587282/76‑1 (MQ=18)
                          gATGCCCTCACCCCGGCG‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCt                                                <  2:5579256/76‑1 (MQ=17)
                          gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                >  2:4392837/1‑76 (MQ=18)
                          gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                >  2:4500325/1‑76 (MQ=18)
                          gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                <  1:3334456/76‑1 (MQ=18)
                           aTGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCtn                                               <  1:2937085/76‑2 (MQ=18)
                              cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTa                                            >  2:5566711/1‑76 (MQ=18)
                              cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                            >  1:753237/1‑76 (MQ=21)
                              cccTCACCCCG‑‑ACCCTCGCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                            <  2:5179180/76‑1 (MQ=21)
                                cTCACCCCG‑‑ACCCTCTCCCACAGGGAGCGGGAGAAAACACTCAACGCCTTCTCCTGGAGCAGGCCTTGCTACTAGt                                          >  2:2920327/1‑76 (MQ=18)
                                cTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                          >  1:4638063/1‑76 (MQ=21)
                                cTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                          <  1:1019727/76‑1 (MQ=21)
                                 tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                         >  1:4734886/1‑76 (MQ=21)
                                 tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                         <  2:3556923/76‑1 (MQ=21)
                                           accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                 <  1:3976284/73‑1 (MQ=38)
                                           accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                 >  2:2578911/4‑76 (MQ=38)
                                            ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                <  1:2183171/74‑1 (MQ=38)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               >  1:1417336/2‑76 (MQ=255)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               <  2:2360700/75‑1 (MQ=255)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               >  2:4641426/2‑76 (MQ=255)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               <  1:2440930/75‑1 (MQ=255)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               <  1:2747566/75‑1 (MQ=255)
                                             cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                               <  1:3971592/75‑1 (MQ=255)
                                               tctcCCGCAGGGAGAGGGAGAAACCACGCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCAc                             <  1:181878/76‑1 (MQ=255)
                                                ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                            <  1:4641427/76‑1 (MQ=255)
                                                 tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                           <  2:1074753/76‑1 (MQ=255)
                                                  cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                          >  2:4331914/1‑76 (MQ=255)
                                                     acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                       >  1:682085/1‑76 (MQ=255)
                                                     acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                       <  1:2920327/76‑1 (MQ=255)
                                                         ggagagggagaAAACACGCAAGGCCTTATCCTGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCACCTTAatgatg                   >  1:1673822/1‑76 (MQ=255)
                                                          gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAGTGAtgt                  <  2:419565/76‑1 (MQ=255)
                                                            gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGGTGCGGTACAGCACCTTAATGATGTGa                >  1:4010810/1‑76 (MQ=255)
                                                             agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat               >  2:3615936/1‑76 (MQ=255)
                                                              gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata              <  1:1466827/76‑1 (MQ=255)
                                                              gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAGTAGTTGCGGTACAGCACCTTAATGATGAGata              >  2:1775617/1‑76 (MQ=255)
                                                               ggagaAAACACTCAAGGCCTACTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGACGTGatac             >  2:563763/1‑75 (MQ=255)
                                                                gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTCTTAGGTGCGGTACAGCATCTTCATGCTGTGATATc            >  1:2860064/1‑76 (MQ=255)
                                                                gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTACTGATGTGATATc            >  2:4638546/1‑76 (MQ=255)
                                                                gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc            <  2:3023906/76‑1 (MQ=255)
                                                                gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc            <  1:360313/76‑1 (MQ=255)
                                                                 agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc           <  1:5812144/76‑1 (MQ=255)
                                                                  gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg          >  2:3104970/1‑76 (MQ=255)
                                                                   aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa         >  1:4421794/1‑76 (MQ=255)
                                                                   aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa         >  2:5297681/1‑76 (MQ=255)
                                                                    aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGAGATCAGaa        >  1:1915995/1‑76 (MQ=255)
                                                                     aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc       <  1:1071117/76‑1 (MQ=255)
                                                                     aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc       <  1:3139331/76‑1 (MQ=255)
                                                                       cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCCCCTTAATGATGTGATATCCGCACTg     >  2:333848/1‑76 (MQ=255)
                                                                       cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTg     <  1:1598999/76‑1 (MQ=255)
                                                                        acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg    >  1:896231/1‑76 (MQ=255)
                                                                         cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg   >  2:5244521/1‑76 (MQ=255)
                                                                          tCAAGCCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGACATCCGCACTgcgg  >  2:784551/1‑73 (MQ=255)
                                                                                |                                                                        
TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGT  >  NC_000913/3957050‑3957194

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: