Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,186 C→T 6.8% A100A (GCG→GCA)  yfcF ← glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,1860CT 6.8% 185.7 / 4.1 59A100A (GCG→GCA) yfcFglutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/2);  total (29/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 = 1.26e-01

GCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAACCTT  >  NC_000913/2420040‑2420334
                                                                                                                                                  |                                                                                                                                                    
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:611467/1‑149 (MQ=255)
                         tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCGTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaa                                                                                                                           >  3:243452/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                                                                                                                                  <  1:23825/127‑1 (MQ=255)
                                        ccGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAg                                                                                                                                  >  2:23825/1‑127 (MQ=255)
                                                 gcgcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCa                                                                                                   >  1:610328/1‑149 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                               >  2:540473/1‑149 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                               >  1:414921/1‑149 (MQ=255)
                                                     ttttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATc                                                                                               <  1:316254/149‑1 (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:485119/1‑149 (MQ=255)
                                                                                                             gATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAtt                                       <  4:330994/149‑1 (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            <  3:273317/149‑1 (MQ=255)
                                                                                                                                        gTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGa            <  1:459254/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)
                                                                                                                                                  cgcACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAACCtt  <  4:461641/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAACCTT  >  NC_000913/2420040‑2420334

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: ‑