Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,727 G→A 100% G36D (GGC→GAC)  yfcG → GSH‑dependent disulfide bond oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,7270GA96.7% 280.4 / ‑0.1 91G36D (GGC→GAC) yfcGGSH‑dependent disulfide bond oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (42/46);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGC  >  NC_000913/2420631‑2420819
                                                                                                |                                                                                            
tctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACgg                                                                                           <  2:250273/100‑1 (MQ=255)
  tATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTc                                                                                         >  2:620531/1‑100 (MQ=255)
   aTTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCa                                                                                        <  1:2206693/100‑1 (MQ=255)
   aTTTCGCCACGACACCCAATGGTAACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACgg                                                                                           <  2:275552/97‑1 (MQ=255)
       cGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGttt                                                                                    <  2:2832089/100‑1 (MQ=255)
         ccccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCg                                                                                  >  2:2468386/1‑100 (MQ=255)
         ccccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCg                                                                                  >  2:1337691/1‑100 (MQ=255)
          cccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGACCTGGATTATCGCTTGATTAAGATAGACCTGAGGAAAGACGATAAGTTTCg                                                                                  >  1:1111168/1‑99 (MQ=255)
           ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                >  2:1558917/1‑100 (MQ=255)
           ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                >  1:2746694/1‑100 (MQ=255)
           ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                <  2:456506/100‑1 (MQ=255)
             gACACCCAATGGTCACAAAATTTCGCTGTTTTTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCg                                                                              <  1:909945/100‑1 (MQ=255)
             gACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCg                                                                              <  1:19436/100‑1 (MQ=255)
              acacCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCgg                                                                             >  1:2502684/1‑100 (MQ=255)
              acacCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCgg                                                                             <  1:2228894/100‑1 (MQ=255)
               cacCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGa                                                                            >  2:2557307/1‑100 (MQ=255)
               cacCCAATGATCACAAAATTACGCTTTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCgg                                                                             <  2:2308769/99‑1 (MQ=255)
                acCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGaa                                                                           >  2:2238328/1‑100 (MQ=255)
                 cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAt                                                                          >  1:2801173/1‑100 (MQ=255)
                 cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGAGAAAAACGATCAGTTTCGTCCGGAAt                                                                          >  1:2632534/1‑100 (MQ=255)
                   cAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAtt                                                                         >  1:107394/1‑99 (MQ=255)
                     aTGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAttttt                                                                      <  2:2090324/100‑1 (MQ=255)
                       ggTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgc                                                                    <  2:2859375/100‑1 (MQ=255)
                        ntCACAAACTTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTCGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAGTTTTgc                                                                    >  2:529913/2‑99 (MQ=255)
                        gTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                   <  1:1542878/100‑1 (MQ=255)
                        gTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                   <  2:1748046/100‑1 (MQ=255)
                        gTCACAAAATTACGCTGTTTCTAGAAGAAGCAGAGCTAGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                   >  2:1613519/1‑100 (MQ=255)
                         tCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcgc                                                                  <  1:309111/100‑1 (MQ=255)
                         tCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcgc                                                                  <  2:698903/100‑1 (MQ=255)
                          cacaAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                 >  2:2502165/1‑100 (MQ=255)
                          cacaAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                 >  2:3011282/1‑100 (MQ=255)
                          cacaAAATTACGCTGTTTCTAGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                 >  2:740283/1‑100 (MQ=255)
                           acaAAGTTACGCTGTTTCTCGAAGAAGCAGAGCTGGTTTATCGCTTGATTAAGTTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgc                                                                    <  2:2479711/96‑1 (MQ=255)
                           acaAAATTACGCTGTTTCAAGAAGAAGCAGAGCTGGATAATCGCTTGATTAAGGTAAACCTGGGGAAAGACGATCAGTTTCGTCCGGAATTTTTGCGCAt                                                                >  2:64005/1‑100 (MQ=255)
                               aaTTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTCAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCg                                                            >  2:920735/1‑100 (MQ=255)
                                 ttACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGcc                                                          <  2:470567/100‑1 (MQ=255)
                                 ttACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGcc                                                          <  1:1256969/100‑1 (MQ=255)
                                   aCGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCca                                                        <  2:2247436/100‑1 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                       >  2:2770032/1‑100 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                       >  1:1928213/1‑100 (MQ=255)
                                        gTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaac                                                   <  1:661841/100‑1 (MQ=255)
                                        gTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAAACGTTCAATTTCGTCCAAAATATTTGAGCATTTCgacaaaaaa                                                    >  1:343635/1‑91 (MQ=255)
                                        gTTTCTCGAAGAAACAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaac                                                   >  1:2060178/1‑100 (MQ=255)
                                          ttCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAATGTATTCCTTGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaacaa                                                 <  1:285863/100‑1 (MQ=255)
                                          ttCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaacaa                                                 <  1:2115167/100‑1 (MQ=255)
                                          ttCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaacaa                                                 >  1:2101779/1‑100 (MQ=255)
                                          ttCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACATGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaacaa                                                 >  2:1488425/1‑100 (MQ=255)
                                           tctcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACaaa                                                >  1:2531330/1‑100 (MQ=255)
                                            ctcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACaaaa                                               >  2:726861/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                             >  2:1258370/1‑100 (MQ=255)
                                                aagaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTcc                                           >  2:1076669/1‑100 (MQ=255)
                                                 agaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCg                                          >  1:2335462/1‑100 (MQ=255)
                                                  gaagCAGAGCTGGATTATCGCTTGATTAAGGTATACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGTATTTCGCCCAACAACAAAATTCCgg                                         <  2:1180046/100‑1 (MQ=255)
                                                   aagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGTGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                        <  1:1488441/100‑1 (MQ=255)
                                                   aagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                        >  1:1196583/1‑100 (MQ=255)
                                                    agCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCa                                       >  1:2837562/1‑100 (MQ=255)
                                                     gCATAGCTGTATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCgg                                         <  1:1595948/97‑1 (MQ=255)
                                                      cAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAt                                     >  2:1241333/1‑100 (MQ=255)
                                                      cAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAt                                     <  1:2805427/100‑1 (MQ=255)
                                                      cAGAGCTGGATTATAGCTTGATTAAGGTAGACCTGGGGAAAGAAGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATCCCGGCaa                                      >  2:2374365/1‑99 (MQ=255)
                                                       agagCTGGATTATCGCTTGATTAAGGTAGTTCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAt                                     <  1:396319/99‑1 (MQ=255)
                                                        gagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttg                                   >  1:1500252/1‑100 (MQ=255)
                                                          gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                 <  1:1548245/100‑1 (MQ=255)
                                                           cTGGATTATCTCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgttg                                <  1:483562/100‑1 (MQ=255)
                                                            tGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGa                               <  1:2662082/100‑1 (MQ=255)
                                                             ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                              <  1:2068599/100‑1 (MQ=255)
                                                             ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                              <  2:2060155/100‑1 (MQ=255)
                                                                   tCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtct                        >  1:503644/1‑100 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtctc                       >  2:1659620/1‑100 (MQ=255)
                                                                     gCTTGATTAAGTTAGTCCTGGGGAAATACGTTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                      <  1:2994381/100‑1 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                      >  1:1182559/1‑100 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                      <  1:608072/100‑1 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                      >  1:2034141/1‑100 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                      <  1:1857180/100‑1 (MQ=255)
                                                                      cTTGATTAAGGTAGACCTGGGGAAAGACAGTCTGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCt                     <  1:806865/100‑1 (MQ=255)
                                                                         gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                  <  1:734171/100‑1 (MQ=255)
                                                                            tAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGAt               <  2:1713639/100‑1 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg              <  2:2314129/100‑1 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg              <  2:1724367/100‑1 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGACTTTCGCCCCACAACAAAATTCCGGCCCTTGTTGCTCCTTCTCCTACCGa                >  2:3026548/1‑98 (MQ=255)
                                                                              aGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATgg             >  2:1018987/1‑100 (MQ=255)
                                                                               ggTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggc            <  1:1558935/100‑1 (MQ=255)
                                                                               ggTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAAAAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggc            >  1:110198/1‑100 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg          <  2:419599/100‑1 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg          <  1:355795/100‑1 (MQ=255)
                                                                                  aGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcggc         <  1:1978704/100‑1 (MQ=255)
                                                                                     ccTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGa       <  2:2994042/99‑1 (MQ=255)
                                                                                     ccTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTACCGATAGCAGCGaa      >  1:472637/1‑100 (MQ=255)
                                                                                       tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc    <  2:2837530/100‑1 (MQ=255)
                                                                                         nggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGc  >  2:1083516/2‑100 (MQ=255)
                                                                                         gggAAAGACGGTCATTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGc  <  2:1884708/100‑1 (MQ=255)
                                                                                                |                                                                                            
TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGC  >  NC_000913/2420631‑2420819

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