Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,346,318 C→T 5.1% A485V (GCG→GTG)  nrdA → ribonucleoside‑diphosphate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,346,3180CT 5.1% 266.9 / 3.6 78A485V (GCG→GTG) nrdAribonucleoside‑diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

TAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGACGG  >  NC_000913/2346172‑2346465
                                                                                                                                                  |                                                                                                                                                   
tAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCgg                                                                                                                                                   >  1:460599/1‑149 (MQ=255)
 aaCCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGt                                                                                                                                                  >  4:160151/1‑149 (MQ=255)
   ccTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTc                                                                                                                                                <  2:625414/149‑1 (MQ=255)
     tgtgCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGTGGTTCGt                                                                                                                                              >  2:466896/1‑149 (MQ=255)
     tgtgCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGTGGTTCGt                                                                                                                                              >  2:624239/1‑149 (MQ=255)
            ggAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTg                                                                                                                                       <  1:185200/149‑1 (MQ=255)
             aagaTAGCCCTGCCGACCCAACAGCTGAAGGACGTCAACGACGAGAACGGTGAAATCGCGCCGAATAAGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGa                                                                                                                                      >  4:348759/2‑149 (MQ=255)
              agaTAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGAc                                                                                                                                     <  3:160851/149‑1 (MQ=255)
              agaTAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGAc                                                                                                                                     <  2:702097/149‑1 (MQ=255)
                aTAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGAcgc                                                                                                                                   >  1:430558/1‑149 (MQ=255)
                    cccTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCgctg                                                                                                                               >  2:662970/1‑149 (MQ=255)
                      cTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCgctgct                                                                                                                             >  4:450864/1‑149 (MQ=255)
                        gCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTgg                                                                                                                           <  4:44062/149‑1 (MQ=255)
                           gACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGAtt                                                                                                                        <  2:306000/149‑1 (MQ=255)
                             ccAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTAt                                                                                                                      >  1:567518/1‑149 (MQ=255)
                              cAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATc                                                                                                                     >  4:647430/1‑149 (MQ=255)
                               aaaCCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCa                                                                                                                    <  2:378437/149‑1 (MQ=255)
                                aaCCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGAc                                                                                                                                     <  3:256763/131‑1 (MQ=255)
                                aaCCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGAc                                                                                                                                     >  4:256763/1‑131 (MQ=255)
                                 aCCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGCTTATCAgg                                                                                                                  >  2:380985/1‑149 (MQ=255)
                                  ccGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGa                                                                                                                 <  2:99334/149‑1 (MQ=255)
                                        aaCGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTAccc                                                                                                           >  1:556252/1‑149 (MQ=255)
                                          cgacgTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGa                                                                                                         <  2:27795/149‑1 (MQ=255)
                                            acgTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATc                                                                                                       >  4:689601/1‑149 (MQ=255)
                                                         gaACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGTGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAAc                                                                                          <  1:624239/149‑1 (MQ=255)
                                                         gaACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGTGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAAc                                                                                          <  1:466896/149‑1 (MQ=255)
                                                         gaACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAAc                                                                                          >  4:218094/1‑149 (MQ=255)
                                                               tGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAg                                                                                    >  2:532895/1‑149 (MQ=255)
                                                                   aTCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGAt                                                                                <  3:514708/149‑1 (MQ=255)
                                                                   aTCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGAt                                                                                <  3:148529/149‑1 (MQ=255)
                                                                          tgtgTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGa                                                                                                                                      >  4:166415/1‑88 (MQ=255)
                                                                          tgtgTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGa                                                                                                                                      <  3:166415/88‑1 (MQ=255)
                                                                          tgtgTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGgtcgt                                                                         >  1:792200/1‑149 (MQ=255)
                                                                           gtgtACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGgtcgtc                                                                        <  1:532895/149‑1 (MQ=255)
                                                                            tgtACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGa                                                                                                                 <  4:59109/107‑1 (MQ=255)
                                                                            tgtACGCTGTCTGCATTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGa                                                                                                                 >  3:59109/1‑107 (MQ=255)
                                                                              tACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTa                                                                     >  2:705715/1‑149 (MQ=255)
                                                                                    gTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTgg                                                               <  2:793873/149‑1 (MQ=255)
                                                                                    gTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTgg                                                               <  3:160151/149‑1 (MQ=255)
                                                                                    gTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTgg                                                               <  4:45311/149‑1 (MQ=255)
                                                                                    gTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTgg                                                               <  2:556252/149‑1 (MQ=255)
                                                                                     tCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCa                                                                                                                    >  2:79610/1‑95 (MQ=255)
                                                                                     tCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCa                                                                                                                    <  1:79610/95‑1 (MQ=255)
                                                                                     tCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTggg                                                              <  3:689601/149‑1 (MQ=255)
                                                                                      cTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGt                                                             <  2:712866/149‑1 (MQ=255)
                                                                                      cTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGt                                                             <  4:31872/149‑1 (MQ=255)
                                                                                       tGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCTCCAAACGTGGAGCGATGGGTCGTCGTACGCCGGGTa                                                            >  4:515547/1‑149 (MQ=255)
                                                                                       tGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTa                                                            >  1:204694/1‑149 (MQ=255)
                                                                                         cTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTAtt                                                          >  3:726579/1‑149 (MQ=255)
                                                                                            tCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGgt                                                       <  3:450864/149‑1 (MQ=255)
                                                                                               aCCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGtgtg                                                    <  1:662970/149‑1 (MQ=255)
                                                                                                  tGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATc                                                 >  1:773897/1‑149 (MQ=255)
                                                                                                       gcAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTa                                                            >  4:65079/1‑133 (MQ=255)
                                                                                                       gcAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTa                                                            <  3:65079/133‑1 (MQ=255)
                                                                                                         aaTTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCg                                          >  1:494381/1‑149 (MQ=255)
                                                                                                          aTTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCACCTTCGc                                         <  2:578219/149‑1 (MQ=255)
                                                                                                          aTTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGc                                         <  3:604626/149‑1 (MQ=255)
                                                                                                             aataaCCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTg                                                                                       >  3:545257/1‑100 (MQ=255)
                                                                                                             aataaCCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTg                                                                                       <  4:545257/100‑1 (MQ=255)
                                                                                                                 aCCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCCGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTtactac                                  >  1:620434/1‑149 (MQ=255)
                                                                                                                       aTGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCg                            <  1:380985/149‑1 (MQ=255)
                                                                                                                              ggAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACg                     <  3:647430/149‑1 (MQ=255)
                                                                                                                               gAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACgg                    <  1:140946/149‑1 (MQ=255)
                                                                                                                               gAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACgg                    >  3:388203/1‑149 (MQ=255)
                                                                                                                                    gCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAAc               <  2:567518/149‑1 (MQ=255)
                                                                                                                                    gCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAAc               <  2:204694/149‑1 (MQ=255)
                                                                                                                                      tGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGc             >  3:320097/1‑149 (MQ=255)
                                                                                                                                        gCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCa                       <  1:640224/137‑1 (MQ=255)
                                                                                                                                        gCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCa                       >  2:640224/1‑137 (MQ=255)
                                                                                                                                        gCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTa           >  3:126619/1‑149 (MQ=255)
                                                                                                                                          aaTTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACt         >  2:174804/1‑149 (MQ=255)
                                                                                                                                               tGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCCAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGAc    >  1:674055/1‑149 (MQ=255)
                                                                                                                                               tGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGAc    <  4:112514/149‑1 (MQ=255)
                                                                                                                                               tGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGAc    <  2:678826/149‑1 (MQ=255)
                                                                                                                                               tGGCGGTTCGTGCACTTGACGCGCTGCCGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGAc    <  2:620434/149‑1 (MQ=255)
                                                                                                                                                ggCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATc                                                 >  2:181300/1‑103 (MQ=255)
                                                                                                                                                ggCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATc                                                 <  1:181300/103‑1 (MQ=255)
                                                                                                                                                 gCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGACgg  >  1:133222/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTGAAATCGCGCTGTGTACGCTGTCTGCTTTCAACCTGGGCGCAATTAATAACCTGGATGAACTGGAAGAGCTGGCAATTCTGGCGGTTCGTGCACTTGACGCGCTGCTGGATTATCAGGATTACCCGATCCCGGCCGCCAAACGTGGAGCGATGGGTCGTCGTACGCTGGGTATTGGTGTGATCAACTTCGCTTACTACCTGGCGAAGCACGGTAAACGCTACTCCGACGG  >  NC_000913/2346172‑2346465

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