Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,182,139 T→G Q90P (CAG→CCG)  ygiD ← 4,5‑ DOPA‑extradiol‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,182,1390TG98.4% 227.2 / ‑3.0 64Q90P (CAG→CCG) ygiD4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (35/28);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

AGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACA  >  NC_000913/3181991‑3182282
                                                                                                                                                    |                                                                                                                                               
aGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCggt                                                                                                                                               >  1:75401/1‑151 (MQ=255)
     tGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCa                                                                                                                                           <  1:519960/150‑1 (MQ=255)
      gCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAgc                                                                                                                                         <  2:914226/151‑1 (MQ=255)
       cTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACAACTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAgcg                                                                                                                                        >  1:752080/1‑151 (MQ=255)
          cTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAgcgc                                                                                                                                       >  2:1103481/1‑149 (MQ=255)
            gTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAgg                                                                                                                                    <  1:80079/150‑1 (MQ=255)
            gTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAg                                                                                                                                     <  1:1011662/149‑1 (MQ=255)
              cGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGa                                                                                                                                 <  1:767259/151‑1 (MQ=255)
                  aCTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAAccc                                                                                                                             >  2:684904/1‑151 (MQ=255)
                  aCTCAACTGCACCATCGGCATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAAccc                                                                                                                             <  1:120690/151‑1 (MQ=255)
                     cAACTGCACCATCGGGATATCAGCTTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGa                                                                                                                          <  1:638631/151‑1 (MQ=255)
                     cAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCgg                                                                                                                           <  1:552040/150‑1 (MQ=255)
                     cAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGa                                                                                                                          <  2:600663/151‑1 (MQ=255)
                      aaCTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGa                                                                                                                          <  1:899431/150‑1 (MQ=255)
                         tGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGc                                                                                                                        >  2:406283/1‑149 (MQ=255)
                         tGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGc                                                                                                                        >  1:1037053/1‑149 (MQ=255)
                          gCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAgg                                                                                                                     <  2:196739/151‑1 (MQ=255)
                          gCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAgg                                                                                                                     >  2:437947/1‑151 (MQ=255)
                            aCCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGAt                                                                                                                   >  2:947795/1‑151 (MQ=255)
                             ccATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATa                                                                                                                  <  2:119540/151‑1 (MQ=255)
                                 cGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGc                                                                                                              >  1:1158169/1‑151 (MQ=255)
                                   ggATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGt                                                                                                            <  2:245121/151‑1 (MQ=255)
                                     atatCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTAt                                                                                                          <  2:1023103/151‑1 (MQ=255)
                                     atatCAGCGTCAGGATACATATTAAACAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTAt                                                                                                          >  1:1269621/1‑151 (MQ=255)
                                       atCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGCGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATc                                                                                                         >  1:1013511/1‑150 (MQ=255)
                                        tCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCg                                                                                                        >  1:658791/1‑150 (MQ=255)
                                          aGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTAc                                                                                                     >  2:851932/1‑151 (MQ=255)
                                           gCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGGACCCGGAGCAGGATAATGCGTATCGTACa                                                                                                    <  2:232986/151‑1 (MQ=255)
                                            cGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAg                                                                                                   <  1:209838/151‑1 (MQ=255)
                                             gTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAgc                                                                                                  >  2:797292/1‑151 (MQ=255)
                                                     aCATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCg                                                                                           <  1:294249/150‑1 (MQ=255)
                                                       aTCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCgg                                                                                          >  1:1178245/1‑149 (MQ=255)
                                                       aTCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAAGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCggg                                                                                         >  1:185385/1‑150 (MQ=255)
                                                           tAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGaa                                                                                       <  1:167660/148‑1 (MQ=255)
                                                           tAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGcc                                                                                    >  2:56950/1‑151 (MQ=255)
                                                             aTCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccac                                                                                  <  1:41752/151‑1 (MQ=255)
                                                               cAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccacc                                                                                 >  2:202582/1‑150 (MQ=255)
                                                                   aCGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGt                                                                            <  2:1235416/151‑1 (MQ=255)
                                                                       cccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATg                                                                        >  2:1128147/1‑151 (MQ=255)
                                                                         cACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGaa                                                                      <  1:215309/151‑1 (MQ=255)
                                                                          aCGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAAt                                                                     >  1:484597/1‑151 (MQ=255)
                                                                               ccGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTgg                                                                >  1:979683/1‑151 (MQ=255)
                                                                                cGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTggg                                                               <  1:33572/151‑1 (MQ=255)
                                                                                          gCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTcc                                                     >  1:263770/1‑151 (MQ=255)
                                                                                              cAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATcg                                                 >  2:433623/1‑151 (MQ=255)
                                                                                                ggCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATcgcg                                               >  1:102741/1‑151 (MQ=255)
                                                                                                   ttcttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTc                                            >  2:331045/1‑151 (MQ=255)
                                                                                                   ttcttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGAACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGt                                             >  1:667996/1‑150 (MQ=255)
                                                                                                    tcttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGCACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCAACGCgg                                              >  2:135828/1‑148 (MQ=255)
                                                                                                           ccAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGt                                    <  1:54287/151‑1 (MQ=255)
                                                                                                           ccAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGt                                    <  2:485488/151‑1 (MQ=255)
                                                                                                             aGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGTGGTCACTCCTGTTc                                  >  2:1292287/1‑151 (MQ=255)
                                                                                                                gTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCAc                               >  1:268265/1‑151 (MQ=255)
                                                                                                                         aTCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAAcc                      >  2:409327/1‑151 (MQ=255)
                                                                                                                         aTCGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAAcc                      <  2:1118151/151‑1 (MQ=255)
                                                                                                                           cGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAAcc                      >  1:401280/1‑149 (MQ=255)
                                                                                                                           cGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAAcc                      >  2:291000/1‑149 (MQ=255)
                                                                                                                           cGGCGCTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCa                     <  1:461130/150‑1 (MQ=255)
                                                                                                                             gcgcTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTg                  <  1:736603/151‑1 (MQ=255)
                                                                                                                              cgcTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTg                  >  1:1188759/1‑150 (MQ=255)
                                                                                                                               gcTAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCAGTCACTCCTGTTCCACGGGTAAACCAGTGa                 <  1:947737/150‑1 (MQ=255)
                                                                                                                                 tAACAGCTCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGGACAGCGCCTGCGGGAAGACACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGcc              >  1:550248/1‑151 (MQ=255)
                                                                                                                                        tCAACCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAacc        >  1:1259924/1‑150 (MQ=255)
                                                                                                                                           aCCAGACGCGGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTAAGCCGAAaccacc     >  2:462553/1‑150 (MQ=255)
                                                                                                                                              aGACGCGGGGCCAGCGCAGGCGAACCCTGAGCAGGATAATGCGCTTCGTACAGCTCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACcaca  >  1:564773/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACA  >  NC_000913/3181991‑3182282

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑