Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,117,238 G→A 44.2% M144I (ATG→ATA yceA → putative rhodanese‑related sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,117,2380GA44.2% 51.7 / 174.0 145M144I (ATG→ATAyceAputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (34/30);  total (75/70)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 147 total aligned reads displayed.

GAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTC  >  NC_000913/1117139‑1117337
                                                                                                   |                                                                                                   
 aaGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGc                                                                                                    <  1:1306898/100‑1 (MQ=255)
   ggTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCtt                                                                                                  >  2:2022351/1‑100 (MQ=255)
   ggTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACtt                                                                                                  >  2:1213376/1‑100 (MQ=255)
    gTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTg                                                                                                 <  1:1497334/100‑1 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCCATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGa                                                                                                >  2:2330013/1‑99 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGAc                                                                                               <  1:2192985/100‑1 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGAc                                                                                               <  1:1829217/100‑1 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGAc                                                                                               <  2:1118683/100‑1 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGAc                                                                                               >  1:4069115/1‑100 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGAc                                                                                               <  1:3449807/100‑1 (MQ=255)
       cgcgATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACg                                                                                              <  2:3321479/100‑1 (MQ=255)
        gcgATCACATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGCAAGTGAACGCCATGCTTGACGa                                                                                             >  1:3860312/1‑100 (MQ=255)
         cgATCGCATTGTTGCCGACGGTGTTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGAt                                                                                            >  1:2483132/1‑100 (MQ=255)
         cgATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGAt                                                                                            >  1:2490740/1‑100 (MQ=255)
         cgATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGAt                                                                                            >  1:3128113/1‑100 (MQ=255)
            tCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATccc                                                                                         <  2:1572761/100‑1 (MQ=255)
            tCGCATTATCGCCGACGATATTGACGATCCTCACTTTGATGCCAGCCATGTTGGTGAGTATCCGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATccc                                                                                         >  1:235759/1‑100 (MQ=255)
             cGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCg                                                                                        <  2:1449751/100‑1 (MQ=255)
              gCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGa                                                                                       >  2:526529/1‑100 (MQ=255)
              gCATTGTTGCAGACGGTATTGACGTTCCTCACTCTGACGTCAGCCACGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGa                                                                                       <  1:1050767/100‑1 (MQ=255)
               cATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGAt                                                                                      >  1:4018115/1‑100 (MQ=255)
                 ttgttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGc                                                                                    <  2:93622/100‑1 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCAc                                                                                  >  2:2484142/1‑100 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCa                                                                                   >  2:935708/1‑99 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCAc                                                                                  <  2:3428460/100‑1 (MQ=255)
                    ttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACt                                                                                 <  2:1711023/100‑1 (MQ=255)
                        cGACGGTATTGACGATCCTCACTTTGATGCCAGCACTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTAttt                                                                             >  2:1761130/1‑100 (MQ=255)
                         gACGGTATTGACGATCCTCACATTGATGCCAGTAATTTTAGTGAGTAACTGCAAGCGGCGGATATGAACGCCATACTTGACGATACCGATGCACTAttt                                                                             <  1:570057/99‑1 (MQ=255)
                           cGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATTCACTCTTTATc                                                                          >  1:2837830/1‑100 (MQ=255)
                           cGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATc                                                                          <  1:934051/100‑1 (MQ=255)
                           cGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATc                                                                          <  1:1025855/100‑1 (MQ=255)
                           cGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATc                                                                          <  2:4168841/100‑1 (MQ=255)
                           cGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATc                                                                          <  1:3981051/100‑1 (MQ=255)
                             gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGa                                                                        >  2:1787391/1‑100 (MQ=255)
                             gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGa                                                                        <  1:1837212/100‑1 (MQ=255)
                              tATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGAc                                                                       <  1:1163726/100‑1 (MQ=255)
                                 tGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATg                                                                    >  1:3814071/1‑100 (MQ=255)
                                  gACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGc                                                                   >  1:1779352/1‑100 (MQ=255)
                                   aCGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCg                                                                  <  1:569193/100‑1 (MQ=255)
                                   aCGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCg                                                                  <  2:1140167/100‑1 (MQ=255)
                                     gATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTa                                                                >  1:237674/1‑100 (MQ=255)
                                        ccTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGTCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAAcc                                                             >  1:335721/1‑100 (MQ=255)
                                         cTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCa                                                            >  1:3599069/1‑100 (MQ=255)
                                           cACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACt                                                          <  2:169355/100‑1 (MQ=255)
                                             cTTTGATGCCAGCAATGTTTGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTAt                                                        >  1:4196943/1‑100 (MQ=255)
                                              tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATg                                                       <  1:3285000/100‑1 (MQ=255)
                                              tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATg                                                       <  2:3249645/100‑1 (MQ=255)
                                                tGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAg                                                     <  1:2676799/100‑1 (MQ=255)
                                                  atgtcagataTGTTGGTGAGTATCTGCAATCGTCGGATGTGAACGCCATGCTTTACCATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAg                                                     <  1:1585595/89‑1 (MQ=255)
                                                  aTGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTa                                                   >  1:3357191/1‑100 (MQ=255)
                                                   tGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTAt                                                  >  2:2938853/1‑100 (MQ=255)
                                                   tGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTAt                                                  >  2:2939568/1‑100 (MQ=255)
                                                   tGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAATCACTATGACTAt                                                  >  1:2551069/1‑100 (MQ=255)
                                                   tGCCAGCAAAGTTGGTGAGTATCTGCAAGCGGCGGACGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATTCGTCACCACTACGAGTAt                                                  >  1:1411910/1‑100 (MQ=255)
                                                     ccAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGa                                                <  2:3142585/100‑1 (MQ=255)
                                                     ccAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGa                                                >  1:4224422/1‑100 (MQ=255)
                                                     ccAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGa                                                >  1:1929355/1‑100 (MQ=255)
                                                      cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGaa                                               >  2:314659/1‑100 (MQ=255)
                                                       aGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCTCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAg                                              <  1:2101908/100‑1 (MQ=255)
                                                           aTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                          <  1:2972657/100‑1 (MQ=255)
                                                           aTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                          >  2:3098770/1‑100 (MQ=255)
                                                           aTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                          <  2:743028/100‑1 (MQ=255)
                                                           aTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                          <  1:526527/100‑1 (MQ=255)
                                                            tGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgggg                                         <  2:1162212/100‑1 (MQ=255)
                                                            tGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgggg                                         <  2:3092096/100‑1 (MQ=255)
                                                            tGTTGGTGAGTATCCGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgggg                                         >  2:3975702/1‑100 (MQ=255)
                                                             gTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGc                                        <  2:3493327/100‑1 (MQ=255)
                                                             gTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGc                                        >  2:1954957/1‑100 (MQ=255)
                                                                 gTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACtt                                    >  2:1834888/1‑100 (MQ=255)
                                                                  tGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACttt                                   >  2:3867092/1‑100 (MQ=255)
                                                                   gAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTg                                  >  2:1638946/1‑100 (MQ=255)
                                                                   gAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTg                                  >  2:1828768/1‑100 (MQ=255)
                                                                    aGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGa                                 >  1:287679/1‑100 (MQ=255)
                                                                     gTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaa                                <  1:662862/100‑1 (MQ=255)
                                                                       aTCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaaaa                              <  2:394756/100‑1 (MQ=255)
                                                                        tCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTCTGAGTATGAAGTGGGGCACTTTGAAAAc                             >  1:457098/1‑100 (MQ=255)
                                                                         ctacAAGCGGCGGAAGTGCACGCCATGCGTGACGATCCCGATGCACTATTTATCGACACGCGTAACCACTCTGCCTATGAAGTGAGGCACTTTGaaaa                              >  1:2969136/4‑98 (MQ=255)
                                                                           gCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgcg                          >  2:2184878/1‑100 (MQ=255)
                                                                              aGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTc                       >  2:738140/1‑100 (MQ=255)
                                                                                 ggcggAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaa                    <  1:1796485/100‑1 (MQ=255)
                                                                                  gcggAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaaa                   <  2:103348/100‑1 (MQ=255)
                                                                                   cggAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAAt                  >  2:331401/1‑100 (MQ=255)
                                                                                       aGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCg              <  2:569308/100‑1 (MQ=255)
                                                                                       aGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCg              >  1:3400320/1‑100 (MQ=255)
                                                                                       aGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCg              >  1:141687/1‑100 (MQ=255)
                                                                                       aGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCg              >  1:3714210/1‑100 (MQ=255)
                                                                                        gTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCgg             <  1:1055829/100‑1 (MQ=255)
                                                                                         tGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGc            <  2:1308803/100‑1 (MQ=255)
                                                                                         tGAACCCCATGCTTGACGCTCCCGATGCAGTATTTATCGACATGCGTAACCCCTATGAGTATGAAGCGGGGCACTTTGAAAACGCGCTCGAAATTCCGGc            <  1:2566946/100‑1 (MQ=255)
                                                                                          gAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa           <  1:2143447/100‑1 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa           >  2:2404949/1‑100 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa           <  2:1178399/100‑1 (MQ=255)
                                                                                           aaCGCCATGCTTGACGATCCCGATGCATTGTTTATCGACGTGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTGGAAATTCCGGCAg          <  1:3002885/100‑1 (MQ=255)
                                                                                              gCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATa       >  1:2759652/1‑100 (MQ=255)
                                                                                               ccATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc      >  2:1601975/1‑100 (MQ=255)
                                                                                               ccATGCTTGACGAGCCCGATGCACAATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc      <  2:4072674/100‑1 (MQ=255)
                                                                                               ccATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc      >  1:4167580/1‑100 (MQ=255)
                                                                                                 aTGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCt    <  2:3998616/100‑1 (MQ=255)
                                                                                                  tGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCtt   <  1:448857/100‑1 (MQ=255)
                                                                                                   gCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTc  <  1:2022347/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTC  >  NC_000913/1117139‑1117337

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: