Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,647,409 T→G 5.2% Y371D (TAT→GAT)  gor → glutathione reductase (NADPH)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,4090TG 5.2% 300.9 / 5.0 92Y371D (TAT→GAT) gorglutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base T (48/39);  new base G (2/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCG  >  NC_000913/3647261‑3647556
                                                                                                                                                    |                                                                                                                                                   
ctctGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGta                                                                                                                                                    >  2:693396/1‑150 (MQ=255)
 tctGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGtat                                                                                                                                                   >  1:1130842/1‑150 (MQ=255)
    gAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTGCTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATaaa                                                                                                                                                >  2:952793/1‑150 (MQ=255)
           tGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCtt                                                                                                                                         >  2:24297/1‑150 (MQ=255)
            gTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCttt                                                                                                                                        >  1:1002147/1‑150 (MQ=255)
              ttAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCa                                                                                                                                      >  1:520973/1‑150 (MQ=255)
              ttAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCa                                                                                                                                      >  1:648973/1‑150 (MQ=255)
               taataaCAATCCGGATGAGCATCTTGATTACAGCAACATTCCGACCGTGTTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCATGCTCGCGAGCAGTATGGCGACGATCAGGTTAAAGTTTATAAATCCTCTTTCAc                                                                                                                                     <  2:1313906/150‑1 (MQ=255)
               taataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCAc                                                                                                                                     >  2:750114/1‑150 (MQ=255)
               taataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCAc                                                                                                                                     >  2:873955/1‑150 (MQ=255)
               taataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCAc                                                                                                                                     <  1:627012/150‑1 (MQ=255)
                 ataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGAGGATTTATACACTTTCAcct                                                                                                                                   <  1:363826/150‑2 (MQ=255)
                  taaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACcgc                                                                                                                                  >  2:1244116/1‑150 (MQ=255)
                   aaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGAGGATTTATACACTTTCAcctga                                                                                                                                 <  1:390821/150‑4 (MQ=255)
                   aaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGAGGATTTATACACTTTCAcctga                                                                                                                                 <  2:1521947/150‑4 (MQ=255)
                     cAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGAt                                                                                                                               >  1:1052733/1‑150 (MQ=255)
                          cGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGtata                                                                                                                          >  1:1411022/1‑150 (MQ=255)
                           ggATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATAc                                                                                                                         >  1:1488685/1‑150 (MQ=255)
                                aGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATAccgccg                                                                                                                    >  2:1198260/1‑150 (MQ=255)
                                 gCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTCGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGt                                                                                                                   >  1:1534535/1‑150 (MQ=255)
                                    tCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTcac                                                                                                                >  1:1307235/1‑150 (MQ=255)
                                     cTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTcacc                                                                                                               <  1:1025812/150‑1 (MQ=255)
                                     cTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTcacc                                                                                                               >  1:149639/1‑150 (MQ=255)
                                     cTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTcacc                                                                                                               <  1:26965/150‑1 (MQ=255)
                                       ggATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTcaccac                                                                                                             >  1:969467/1‑150 (MQ=255)
                                        gATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACt                                                                                                            <  1:990929/150‑1 (MQ=255)
                                             cAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCg                                                                                                       >  2:1108309/1‑150 (MQ=255)
                                             cAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCg                                                                                                       >  1:924348/1‑150 (MQ=255)
                                               gCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGcc                                                                                                     >  2:265715/1‑150 (MQ=255)
                                                cAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCa                                                                                                    >  2:382608/1‑150 (MQ=255)
                                                  aCATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCtt                                                                                                                                         <  1:1482407/111‑1 (MQ=255)
                                                  aCATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCtt                                                                                                                                         >  2:1482407/1‑111 (MQ=255)
                                                   cATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGcc                                                                                                 <  2:605255/150‑1 (MQ=255)
                                                        cGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGtata                                                                                                                                                  >  1:982727/1‑96 (MQ=255)
                                                        cGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGtata                                                                                                                                                  <  2:982727/96‑1 (MQ=255)
                                                        cGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGcc                                                                                            >  2:944355/1‑150 (MQ=255)
                                                            cGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCAt                                                                                        <  1:783485/150‑1 (MQ=255)
                                                            cGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCAt                                                                                        <  1:1373709/150‑1 (MQ=255)
                                                             gTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCATGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATg                                                                                       <  2:888815/150‑1 (MQ=255)
                                                              tGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGa                                                                                      <  1:402618/150‑1 (MQ=255)
                                                               ggTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGaa                                                                                     >  2:1360196/1‑150 (MQ=255)
                                                                gTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGCTGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAg                                                                                    <  1:1347917/150‑1 (MQ=255)
                                                                  cTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCGGCCGTGCCGCATGAAGCt                                                                                  >  2:801207/1‑150 (MQ=255)
                                                                  cTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGc                                                                                   >  2:1137325/1‑149 (MQ=255)
                                                                  cTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGc                                                                                   <  1:1137325/149‑1 (MQ=255)
                                                                      aGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGgtg                                                                              >  2:1397335/1‑150 (MQ=255)
                                                                            ccgccgATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGtt                                                                        <  1:898917/150‑1 (MQ=255)
                                                                             cgccgATTGGTACTCGTGGTTTCACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTTGGGTGCGTTg                                                                       <  2:1032928/150‑1 (MQ=255)
                                                                                 gATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                   <  1:1244116/150‑1 (MQ=255)
                                                                                 gATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                   >  1:1235249/1‑150 (MQ=255)
                                                                                  aTTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCg                                                                  <  1:1198260/150‑1 (MQ=255)
                                                                                    tGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGga                                                                >  1:1031696/1‑150 (MQ=255)
                                                                                     ggTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAAGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaa                                                               <  2:1145096/150‑1 (MQ=255)
                                                                                      gTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                              <  2:520973/150‑1 (MQ=255)
                                                                                      gTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                              <  2:52510/150‑1 (MQ=255)
                                                                                       tACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaaga                                                             >  1:5383/1‑150 (MQ=255)
                                                                                        aCTGTTGGTTTGACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagag                                                            <  2:794518/150‑1 (MQ=255)
                                                                                        aCTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagag                                                            <  1:944355/150‑1 (MQ=255)
                                                                                           gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaag                                                         >  1:74353/1‑150 (MQ=255)
                                                                                           gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaag                                                         >  1:720650/1‑150 (MQ=255)
                                                                                               gTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTg                                                     >  1:453321/1‑150 (MQ=255)
                                                                                                tttAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACAGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATAGGAAGAGAAGATTGt                                                    >  1:1517405/1‑150 (MQ=255)
                                                                                                   aaCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCgg                                                 <  2:1535716/150‑1 (MQ=255)
                                                                                                   aaCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCgg                                                 <  1:24297/150‑1 (MQ=255)
                                                                                                    aCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGt                                                >  2:1164078/1‑150 (MQ=255)
                                                                                                      ggAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGTATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTAt                                              <  2:1130842/150‑1 (MQ=255)
                                                                                                         aCCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCa                                           >  1:710295/1‑150 (MQ=255)
                                                                                                          ccGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCAc                                          >  1:531557/1‑150 (MQ=255)
                                                                                                            gCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACgg                                        <  1:382608/150‑1 (MQ=255)
                                                                                                            gCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACgg                                        <  1:952793/150‑1 (MQ=255)
                                                                                                             cAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGc                                       >  1:583089/1‑150 (MQ=255)
                                                                                                                 cgcgcgAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTg                                   <  1:873955/150‑1 (MQ=255)
                                                                                                                   cgcgAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGc                                 >  1:1537954/1‑150 (MQ=255)
                                                                                                                     cgAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCtt                               >  1:490275/1‑150 (MQ=255)
                                                                                                                       aGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTg                             >  2:1480743/1‑150 (MQ=255)
                                                                                                                       aGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTg                             <  1:741803/150‑1 (MQ=255)
                                                                                                                        gCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTgg                            <  2:1411022/150‑1 (MQ=255)
                                                                                                                        gCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTgg                            <  1:1189350/150‑1 (MQ=255)
                                                                                                                          aGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTa                          >  2:1218963/1‑150 (MQ=255)
                                                                                                                          aGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTa                          >  1:199410/1‑150 (MQ=255)
                                                                                                                           gTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTAt                         >  2:540492/1‑150 (MQ=255)
                                                                                                                            tATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATg                        <  1:1227885/150‑1 (MQ=255)
                                                                                                                             aTGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGTTGCATGAATTTTTTGTGCGTTGTATGTGATGAGAAGTTTTTGTGTATTCTTGGCATTTGGTTTGGTATgg                       <  2:418756/150‑1 (MQ=255)
                                                                                                                             aTGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATgg                       <  2:1523569/150‑1 (MQ=255)
                                                                                                                             aTGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATgg                       >  1:257007/1‑150 (MQ=255)
                                                                                                                               ggCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGAc                     >  1:69158/1‑150 (MQ=255)
                                                                                                                                gCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACg                    >  2:1234319/1‑150 (MQ=255)
                                                                                                                                 cgacgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGa                   <  2:1031696/150‑1 (MQ=255)
                                                                                                                                   acgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGaaa                 <  1:1108309/150‑1 (MQ=255)
                                                                                                                                     gaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATg               <  1:1234319/150‑1 (MQ=255)
                                                                                                                                        cAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTg            >  1:1367516/1‑150 (MQ=255)
                                                                                                                                         cggTGAAAGAGGATTTATACACTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc           >  1:964201/2‑150 (MQ=255)
                                                                                                                                         cggTGAAAGAGGATTTATACACTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc           >  1:393963/2‑150 (MQ=255)
                                                                                                                                         aGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc           <  1:265715/150‑1 (MQ=255)
                                                                                                                                          ggTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCa          >  1:698565/1‑150 (MQ=255)
                                                                                                                                            tGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAgg        <  1:591375/150‑1 (MQ=255)
                                                                                                                                                  tgtATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcg  >  2:1333806/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCG  >  NC_000913/3647261‑3647556

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