Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,369,793 G→A E35E (GAG→GAA)  pspE → thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,7930GA99.2% 407.7 / ‑4.7 124E35E (GAG→GAA) pspEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (56/65);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 127 total aligned reads displayed.

AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCA  >  NC_000913/1369696‑1369890
                                                                                                 |                                                                                                 
aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACa                                                                                                 >  1:3036318/1‑100 (MQ=255)
   aaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGt                                                                                              <  1:3442058/100‑1 (MQ=255)
     ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTc                                                                                            >  2:2842215/1‑100 (MQ=255)
      gCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCa                                                                                           >  2:3283721/1‑100 (MQ=255)
      gCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCa                                                                                           <  2:2107439/100‑1 (MQ=255)
        ttACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAgg                                                                                         <  1:1761018/100‑1 (MQ=255)
           cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAggggg                                                                                      <  2:3437543/100‑1 (MQ=255)
             tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGcc                                                                                    >  2:1585506/1‑100 (MQ=255)
             tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGAACGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGc                                                                                     <  1:2626293/99‑1 (MQ=255)
              gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCa                                                                                   <  2:272645/100‑1 (MQ=255)
                 ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCa                                                                                <  2:1820919/100‑1 (MQ=255)
                   ggCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCaan                                                                              <  2:32554/100‑2 (MQ=255)
                   ggCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAat                                                                              >  1:3113394/1‑100 (MQ=255)
                   ggCACTGGTGTTTTCACTGCCCGTTTGCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAat                                                                              <  1:3301374/100‑1 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                       <  1:2814780/100‑1 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                       <  2:2638753/100‑1 (MQ=255)
                            gtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTg                                                                     >  1:1534004/1‑100 (MQ=255)
                            gtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTg                                                                     <  1:937999/100‑1 (MQ=255)
                            gtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTg                                                                     >  2:2069738/1‑100 (MQ=255)
                             ttttCGCTGCCCGTTTTCGCCGCTGAATCCTGCATCGATGTTCGGGTTCCAGAGCAGTAGCAGCAAGAACACGTTCAGGGGGCCATCAATATGCCCCTga                                                                    <  1:2138043/100‑1 (MQ=255)
                             ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAAAATACGTTCAGGGGGCCATCAATATTCCCCTga                                                                    <  1:2278509/100‑1 (MQ=255)
                              tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                   >  2:2565820/1‑100 (MQ=255)
                              tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                   >  2:2404451/1‑100 (MQ=255)
                               ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                  >  2:89071/1‑100 (MQ=255)
                                  aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaa                                                               <  2:1457440/100‑1 (MQ=255)
                                   cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                              >  1:2333608/1‑100 (MQ=255)
                                    tGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGt                                                             >  2:147912/1‑100 (MQ=255)
                                     gCCCTTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                            <  1:3149557/100‑1 (MQ=255)
                                     gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                            >  2:1192676/1‑100 (MQ=255)
                                     gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                            <  2:2149407/100‑1 (MQ=255)
                                     gCCCGTTTTCGCCGCTGAACACTGGAGCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaa                                                               <  1:421196/97‑1 (MQ=255)
                                        cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa                                                         >  1:358828/1‑100 (MQ=255)
                                         gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAag                                                        <  1:3423750/100‑1 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                   >  1:630859/1‑100 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                   >  1:492152/1‑100 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                   >  1:2422551/1‑100 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                   >  1:1401211/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                  >  1:263600/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                  >  2:377171/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTACCCTGAAGGAAGTGAAAGAGCGCa                                                  >  1:1177760/1‑100 (MQ=255)
                                                ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                                 <  1:1120387/100‑1 (MQ=255)
                                                ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                                 >  2:3022429/1‑100 (MQ=255)
                                                     naaaacTGGATCGATGTTCGTGGTCCAGAGAAGTATCAGCATGGCCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGGAGTGAAAGAGCGCAt                                                 >  1:1883580/5‑95 (MQ=255)
                                                     gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                            >  2:346898/1‑100 (MQ=255)
                                                     gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                            >  1:2868281/1‑100 (MQ=255)
                                                     gAACACTGGATCGATGCTCGTGTTCCAGAGCAGTATCAGCAAGAACACGGTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                            >  2:2385961/1‑100 (MQ=255)
                                                       acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                          <  2:2379435/100‑1 (MQ=255)
                                                       acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                          >  1:2019904/1‑100 (MQ=255)
                                                       acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                          >  2:2388996/1‑100 (MQ=255)
                                                        cacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                         >  2:640568/1‑100 (MQ=255)
                                                         acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc                                        <  2:1900163/100‑1 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGCACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       >  1:3018851/1‑100 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       <  2:1732475/100‑1 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       <  2:1113585/100‑1 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       >  1:1819359/1‑100 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       >  2:2032038/1‑100 (MQ=255)
                                                          cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                       <  2:2505752/100‑1 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt                                     >  1:457333/1‑100 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATAAATATTCCCGTGAAGGAAGTGAAAGAGCGCATTGCAACCGCCGt                                     >  2:2773302/1‑100 (MQ=255)
                                                             gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt                                    >  2:569442/1‑100 (MQ=255)
                                                               tCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTcc                                  >  1:2728616/1‑100 (MQ=255)
                                                               tCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTcc                                  <  2:2838584/100‑1 (MQ=255)
                                                                cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                 <  2:1471684/100‑1 (MQ=255)
                                                                 gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                                <  1:843051/100‑1 (MQ=255)
                                                                 gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                 <  1:3277612/99‑1 (MQ=255)
                                                                   tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                              <  1:718821/100‑1 (MQ=255)
                                                                   tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                              <  1:3370380/100‑1 (MQ=255)
                                                                    gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAtn                             <  2:1060452/100‑2 (MQ=255)
                                                                    gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                             >  1:642073/1‑100 (MQ=255)
                                                                    gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                             <  1:2984295/100‑1 (MQ=255)
                                                                     ttCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaa                            >  2:738504/1‑100 (MQ=255)
                                                                      tCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaca                           <  2:3215969/100‑3 (MQ=255)
                                                                      tCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaa                           >  2:313978/1‑100 (MQ=255)
                                                                        gtgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaaaa                         <  1:2709033/100‑1 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                        <  2:51786/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                       <  2:65705/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                       >  2:649900/1‑100 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                       >  2:1980261/1‑100 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                       >  2:918764/1‑100 (MQ=255)
                                                                            tCCAGAGCAGTATCAGCAAGATCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                     <  2:1273862/100‑1 (MQ=255)
                                                                            tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                     <  1:2871283/100‑1 (MQ=255)
                                                                            tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCAGAtcggaagcgc                     >  2:308005/1‑91 (MQ=255)
                                                                             ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                    <  2:668002/100‑1 (MQ=255)
                                                                              cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                   >  1:726753/1‑100 (MQ=255)
                                                                              cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                   <  2:492157/100‑1 (MQ=255)
                                                                                  gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg               <  2:166602/100‑1 (MQ=255)
                                                                                    aGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa             <  2:2085896/100‑1 (MQ=255)
                                                                                     gTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa            <  2:2973509/100‑1 (MQ=255)
                                                                                     gTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa            >  2:1420727/1‑100 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGATCGCATTGCCACCGCCGTGCCGGATAAAAACGACTCTGCGaaac           <  1:1205947/100‑2 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg           <  1:1177225/100‑1 (MQ=255)
                                                                                       aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt          >  2:2610542/1‑100 (MQ=255)
                                                                                       aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGACACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt          >  2:2437443/1‑100 (MQ=255)
                                                                                         cAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgaga        >  2:2147141/1‑97 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCGATCAATATTCCCCGGAAGGAAGTAAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACCCCGTGAAAGTGTa       <  1:2132300/100‑1 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCGCGATAAAAACGACACCGTGAAAGAGTa       >  1:183185/1‑100 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt      <  1:1490453/100‑1 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGAGCATTGCCAGCGCCGCTCCGGATATAAACGACACCGTGAAAGAGTAt      >  1:2349767/1‑100 (MQ=255)
                                                                                              agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc   >  1:2980156/1‑100 (MQ=255)
                                                                                               gagCACGTTCAGGGGGCCATCAATATTGTCCTTAAGGAAGTGAAAGCGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCa  <  1:1958681/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCA  >  NC_000913/1369696‑1369890

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