Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,183 A→G 6.6% R101R (CGT→CGC)  yfcF ← glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,1830AG 6.6% 196.7 / 4.2 61R101R (CGT→CGC) yfcFglutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

GCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAC  >  NC_000913/2420040‑2420331
                                                                                                                                               |                                                                                                                                                    
gCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACgcgc                                                                                                                                                 >  3:461641/1‑149 (MQ=255)
 ccATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGAAGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCa                                                                                                                                                >  2:222047/1‑149 (MQ=255)
    tCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGa                                                                                                                                             >  4:256342/1‑149 (MQ=255)
    tCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGa                                                                                                                                             <  3:510820/149‑1 (MQ=255)
     cgcgAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGac                                                                                                                                            <  3:187273/149‑2 (MQ=255)
      gcgAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtt                                                                                                                                           <  2:180233/149‑1 (MQ=255)
         aaCAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTc                                                                                                                                        <  1:197200/149‑1 (MQ=255)
            aGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTaa                                                                                                                                     >  3:169971/1‑149 (MQ=255)
              tttCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAAt                                                                                                                                   >  1:206800/1‑149 (MQ=255)
                         tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGGCGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGATTCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaa                                                                                                                        >  1:238919/1‑149 (MQ=255)
                         tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCGTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaa                                                                                                                        >  3:243452/1‑149 (MQ=255)
                         tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCGTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaa                                                                                                                        >  3:611467/1‑149 (MQ=255)
                            cTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATa                                                                                                                     <  3:256342/149‑1 (MQ=255)
                             ttttCCATAGGGCGTTAGTGGCGCATTTTTTGACCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGAGCATCAGTCCGCTGCGCCGCCAGGCCTGCATCTGTCGTGCACGCGCACGATTTTCTAAATCATGAGGATAAAtcc                                                                                                                    >  2:75525/1‑147 (MQ=255)
                                 cccTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTc                                                                                                                >  3:330994/1‑149 (MQ=255)
                                    tCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCa                                                                                                             <  4:538252/149‑1 (MQ=255)
                                    tCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCa                                                                                                             <  1:257424/149‑1 (MQ=255)
                                        ccGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAg                                                                                                                               >  2:23825/1‑127 (MQ=255)
                                        ccGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAg                                                                                                                               <  1:23825/127‑1 (MQ=255)
                                          gTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCAc                                                                                                                                                     <  2:163621/103‑1 (MQ=255)
                                          gTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCAc                                                                                                                                                     >  1:163621/1‑103 (MQ=255)
                                                 gcgcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCa                                                                                                >  1:610328/1‑149 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                            >  1:414921/1‑149 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                            <  1:316254/149‑1 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                            >  2:540473/1‑149 (MQ=255)
                                                        tttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGAt                                                                                         <  2:210489/149‑1 (MQ=255)
                                                            gCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTc                                                                                     <  3:584668/149‑1 (MQ=255)
                                                              ccccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCa                                                                                   >  1:98794/1‑149 (MQ=255)
                                                                ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGa                                                                                 <  3:368582/149‑1 (MQ=255)
                                                                       gACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCg                                                                          >  2:567662/1‑149 (MQ=255)
                                                                          aacaTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCaa                                                                       <  3:239241/149‑1 (MQ=255)
                                                                          aacaTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCaa                                                                       >  1:133111/1‑149 (MQ=255)
                                                                             aTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATgg                                                                    >  2:522978/1‑149 (MQ=255)
                                                                               ccGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCa                                                                  <  2:441339/149‑1 (MQ=255)
                                                                                  tCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaa                                                               <  2:610328/149‑1 (MQ=255)
                                                                                      aCGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                      >  3:297357/1‑90 (MQ=255)
                                                                                      aCGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                      <  4:297357/90‑1 (MQ=255)
                                                                                          tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                       >  1:107024/1‑149 (MQ=255)
                                                                                                 ggATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGAtt                                                           >  4:264100/1‑138 (MQ=255)
                                                                                                 ggATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGAtt                                                           <  3:264100/138‑1 (MQ=255)
                                                                                                    tGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCaaaa                                             <  1:64017/149‑1 (MQ=255)
                                                                                                         aTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcatc                                        <  2:182494/149‑1 (MQ=255)
                                                                                                             gATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAtt                                    <  4:330994/149‑1 (MQ=255)
                                                                                                             gATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAtt                                    >  4:485119/1‑149 (MQ=255)
                                                                                                             gATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATAGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAtt                                    >  2:459254/1‑149 (MQ=255)
                                                                                                              aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                   <  2:452724/149‑1 (MQ=255)
                                                                                                                 gCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGta                                >  2:370723/1‑149 (MQ=255)
                                                                                                                  cTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaa                               <  1:540473/149‑1 (MQ=255)
                                                                                                                    gcgcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaata                             >  4:164807/1‑149 (MQ=255)
                                                                                                                        aGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACgg                         <  2:600829/149‑1 (MQ=255)
                                                                                                                         gCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGc                        <  1:522978/149‑1 (MQ=255)
                                                                                                                           cAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGcac                      >  2:615843/1‑149 (MQ=255)
                                                                                                                               ccTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACgg                  >  4:49570/1‑149 (MQ=255)
                                                                                                                                    aTCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgtgt             <  2:206800/149‑1 (MQ=255)
                                                                                                                                     aaaatCGTGCGCGTGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTc            >  2:57377/5‑149 (MQ=255)
                                                                                                                                       aatCGTGCGCGTGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGcaca                     >  4:164863/3‑138 (MQ=255)
                                                                                                                                       aatCGTGCGCGTGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGcaca                     <  3:164863/136‑1 (MQ=255)
                                                                                                                                        gTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGa         <  1:459254/149‑1 (MQ=255)
                                                                                                                                        gTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGa         <  3:273317/149‑1 (MQ=255)
                                                                                                                                        gTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGa         >  4:80858/1‑149 (MQ=255)
                                                                                                                                          cGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGAcc       <  2:238919/149‑1 (MQ=255)
                                                                                                                                               aCGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAc  <  2:107024/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAC  >  NC_000913/2420040‑2420331

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