Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,872,822 T→C A207A (GCA→GCG cysN ← ATP‑sulfurylase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,872,8220TC83.0% 665.2 / 132.2 306A207A (GCA→GCGcysNATP‑sulfurylase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base C (120/134);  total (146/160)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 310 total aligned reads displayed.

TGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGTCTC  >  NC_000913/2872730‑2872911
                                                                                            |                                                                                         
  gATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGc                                                                                         <  2:2005176/93‑1 (MQ=255)
  gATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAgcgc                                                                                         <  2:353967/93‑1 (MQ=255)
   aTCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCa                                                                                        <  1:2841364/93‑1 (MQ=255)
    tctcCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCag                                                                                       >  2:2187061/1‑93 (MQ=255)
     ctcCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCaga                                                                                      <  1:1072538/93‑1 (MQ=255)
     ctcCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCaga                                                                                      <  1:1293065/93‑1 (MQ=255)
      tcCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCagag                                                                                     >  2:686804/1‑93 (MQ=255)
        cACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCagagag                                                                                   <  2:5834504/93‑1 (MQ=255)
         aCGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGc                                                                                  >  1:4727633/1‑93 (MQ=255)
         aCGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGc                                                                                  <  1:2185879/93‑1 (MQ=255)
           ggTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCgg                                                                                <  2:2500120/93‑1 (MQ=255)
            gTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGc                                                                               >  2:2212238/1‑93 (MQ=255)
             tttCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGca                                                                              <  1:4377551/93‑1 (MQ=255)
              ttCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGcac                                                                             >  2:4189317/1‑93 (MQ=255)
              ttCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGcac                                                                             <  1:2693482/93‑1 (MQ=255)
               tCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGcaca                                                                            >  1:6799480/1‑93 (MQ=255)
                 cAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACaaa                                                                          <  2:1021979/93‑1 (MQ=255)
                  aGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAg                                                                         <  2:2075841/93‑1 (MQ=255)
                  aGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAg                                                                         <  1:2266064/93‑1 (MQ=255)
                  aGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAg                                                                         >  1:1558706/1‑93 (MQ=255)
                     aCTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCgg                                                                      <  1:4026964/93‑1 (MQ=255)
                     aCTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCgg                                                                      >  1:7015575/1‑93 (MQ=255)
                      cTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTGCCAGCGCAGAGAGCGGCACAAAGCGGa                                                                     <  2:6425546/93‑1 (MQ=255)
                        tCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGata                                                                   >  2:1867629/1‑93 (MQ=255)
                        tCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGata                                                                   <  1:5552388/93‑1 (MQ=255)
                        tCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGata                                                                   <  2:5601512/93‑1 (MQ=255)
                         cGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGatat                                                                  >  2:1413758/1‑93 (MQ=255)
                          gAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATc                                                                 >  2:522020/1‑93 (MQ=255)
                          gAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATc                                                                 <  1:1911021/93‑1 (MQ=255)
                              aGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGa                                                             >  2:5431108/1‑93 (MQ=255)
                              aGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGa                                                             <  1:6065826/93‑1 (MQ=255)
                                tgtCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGAtt                                                           <  1:2929961/93‑1 (MQ=255)
                                    ggACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTAccc                                                       <  1:2161638/93‑1 (MQ=255)
                                      aCCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCgg                                                     <  1:3372773/93‑1 (MQ=255)
                                       ccGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGc                                                    >  2:7016710/1‑93 (MQ=255)
                                        cGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                   >  1:570528/1‑93 (MQ=255)
                                        cGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                   >  2:2953424/1‑93 (MQ=255)
                                        cGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                   <  1:4938132/93‑1 (MQ=255)
                                        cGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                   <  1:4164182/93‑1 (MQ=255)
                                        cGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                   <  1:2372086/93‑1 (MQ=255)
                                          cTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCagcagc                                                 <  2:5903126/93‑1 (MQ=255)
                                          cTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGc                                                 <  2:3417830/93‑1 (MQ=255)
                                          cTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGc                                                 <  2:3321995/93‑1 (MQ=255)
                                           tGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCt                                                <  2:4068108/93‑1 (MQ=255)
                                            gTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTg                                               >  1:950656/1‑93 (MQ=255)
                                             tACCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGc                                              <  2:6559763/93‑1 (MQ=255)
                                              aCCACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGcc                                             >  2:6094769/1‑93 (MQ=255)
                                               ccACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGccc                                            <  1:2537795/93‑1 (MQ=255)
                                                cACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCg                                           >  2:4749651/1‑93 (MQ=255)
                                                cACGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCg                                           <  2:4353513/93‑1 (MQ=255)
                                                 aCGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCgg                                          >  1:299438/1‑93 (MQ=255)
                                                 aCGGCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCgg                                          >  1:1610398/1‑93 (MQ=255)
                                                  cGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGc                                         <  2:6179287/93‑1 (MQ=255)
                                                   ggCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCa                                        >  1:5875794/1‑93 (MQ=255)
                                                    gCATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCaa                                       >  1:1798464/1‑93 (MQ=255)
                                                     cATACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCaaa                                      <  2:662225/93‑1 (MQ=255)
                                                      aTACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                                     <  2:1117408/93‑1 (MQ=255)
                                                      aTACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                                     <  2:3047689/93‑1 (MQ=255)
                                                      aTACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAg                                     <  2:785300/93‑1 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    >  2:5094876/1‑93 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    >  2:2592467/1‑93 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    >  2:162365/1‑93 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    >  2:2431506/1‑93 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    <  1:851411/93‑1 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAgg                                    >  1:6882522/1‑93 (MQ=255)
                                                       tACTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAACTGCCCGGCAAAgg                                    >  2:384927/1‑93 (MQ=255)
                                                        aCTTTCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGt                                   >  2:596864/1‑93 (MQ=255)
                                                            tCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCaaa                               <  1:4639236/93‑1 (MQ=255)
                                                            tCACTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCaaa                               <  1:508091/93‑1 (MQ=255)
                                                              aCTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAaata                             >  1:7183283/1‑93 (MQ=255)
                                                               cTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAaataa                            >  2:2737963/1‑93 (MQ=255)
                                                               cTTTGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAaataa                            >  2:6705419/1‑93 (MQ=255)
                                                                tttGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAaataat                           >  2:4482043/1‑93 (MQ=255)
                                                                 ttGCGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATc                          >  1:7077685/1‑93 (MQ=255)
                                                                    cGATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTc                       <  1:5617882/93‑1 (MQ=255)
                                                                     gATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCa                      >  1:5225342/1‑93 (MQ=255)
                                                                     gATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCa                      <  2:5691837/93‑1 (MQ=255)
                                                                     gATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCa                      >  1:1707077/1‑93 (MQ=255)
                                                                     gATGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCa                      <  1:7429729/93‑1 (MQ=255)
                                                                      aTGCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCAc                     >  1:1680118/1‑93 (MQ=255)
                                                                        gCCACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGa                   >  2:832587/1‑93 (MQ=255)
                                                                         ccACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGaa                  <  1:76853/93‑1 (MQ=255)
                                                                         ccACGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGaa                  >  1:4478145/1‑93 (MQ=255)
                                                                            cGTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATAc               <  1:6496770/93‑1 (MQ=255)
                                                                             gTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACg              <  2:5714077/93‑1 (MQ=255)
                                                                             gTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACg              >  2:3652225/1‑93 (MQ=255)
                                                                             gTTGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACg              <  2:2763756/93‑1 (MQ=255)
                                                                              ttGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACgg             >  1:5982036/1‑93 (MQ=255)
                                                                              ttGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACgg             >  1:2601412/1‑93 (MQ=255)
                                                                              ttGCCGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACgg             <  1:3541005/93‑1 (MQ=255)
                                                                               tgccGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACggg            <  2:4889389/90‑1 (MQ=255)
                                                                                gccGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGt           <  2:4624265/91‑1 (MQ=255)
                                                                                gccGCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGt           <  1:3737724/91‑1 (MQ=255)
                                                                                 tCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTg          <  2:7081834/93‑1 (MQ=255)
                                                                                   gCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGaa        >  2:6217168/1‑93 (MQ=255)
                                                                                   gCCTTCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGaa        >  1:6315273/1‑93 (MQ=255)
                                                                                      ttCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGt     >  1:3288288/1‑93 (MQ=255)
                                                                                      ttCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGt     >  1:5158263/1‑93 (MQ=255)
                                                                                      ttCCAGCGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGt     <  1:1183715/93‑1 (MQ=255)
                                                                                       tCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGtc    <  1:2017959/93‑1 (MQ=255)
                                                                                            |                                                                                         
TGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGTCTC  >  NC_000913/2872730‑2872911

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