Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,793,897 Δ1 bp 100% coding (213/552 nt) btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,8970T.98.4% 284.4 / ‑2.9 63coding (213/552 nt)btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (33/29);  minor base C (0/1);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793750‑1794040
                                                                                                                                                   |                                                                                                                                               
cgcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtg                                                                                                                                                <  1:1226463/149‑2 (MQ=255)
         ggTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                      >  1:1039504/1‑149 (MQ=255)
          gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGGCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCCTCGCTGCCCGGTTC‑TGTTCCAGAAAc                                                                                                                                     <  1:564632/149‑1 (MQ=255)
            cGGCGCTGCGGCAATCAATTCTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAAGCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTg                                                                                                                                   <  1:584369/149‑1 (MQ=255)
                gcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTg                                                                                                                               >  1:608346/1‑149 (MQ=255)
                 cTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTGC‑TGTTCCAGAAACTGGTTGc                                                                                                                              >  1:507825/1‑149 (MQ=255)
                   gCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCATTCACAACGTACATTGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCAc                                                                                                                            <  1:339373/149‑1 (MQ=255)
                     ggCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                          >  1:437538/1‑149 (MQ=255)
                      gCAATCAATTTTTGATACAGCGGATGGCGTCCTTAGCCATTACCTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGGGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCTAACAACTGGTTGTATggg                                                                                                                         >  2:509558/1‑149 (MQ=255)
                          tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                     >  2:920664/1‑149 (MQ=255)
                             atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGGTCCCGAAACTGGTTGCACGGGAATCCCa                                                                                                                  >  2:989203/1‑149 (MQ=255)
                                        aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                       <  1:843174/149‑1 (MQ=255)
                                            gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTCCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                   >  1:1287050/1‑149 (MQ=255)
                                              tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                                                                 <  1:1202250/149‑1 (MQ=255)
                                              tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCGCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCGTAAACCCtcga                                                                                                 <  2:1336319/149‑1 (MQ=255)
                                                gCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatc                                                                                               <  1:1260684/149‑1 (MQ=255)
                                                gCGTCCTTCCCCATTAACTTCAACCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatc                                                                                               <  1:1260652/149‑1 (MQ=255)
                                                   tCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAc                                                                                            <  1:613074/149‑1 (MQ=255)
                                                             ttAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTc                                                                                  >  1:312858/1‑149 (MQ=255)
                                                             ttAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTc                                                                                  <  1:285120/149‑1 (MQ=255)
                                                              tAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCCTGGGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                                 <  2:787545/149‑1 (MQ=255)
                                                              tAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                                 >  2:659194/1‑149 (MQ=255)
                                                                   tCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAata                                                                            >  1:1244449/1‑149 (MQ=255)
                                                                       tCTTACTGAACGTCGGGAACGTCCCCCCCCATGTGGCGGTAGAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtct                                                                        <  2:390297/149‑1 (MQ=255)
                                                                       tCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAGACCTCGATCGACCCAGGCTTTCTGAATATtct                                                                        <  2:390087/149‑1 (MQ=255)
                                                                            cTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAAc                                                                   >  1:1008295/1‑149 (MQ=255)
                                                                             tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTGCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                                  <  2:881124/149‑1 (MQ=255)
                                                                               aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                                >  2:1087827/1‑149 (MQ=255)
                                                                               aaCATCGGGAACGTCACCCCCCATGGGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                                <  1:69032/149‑1 (MQ=255)
                                                                                aCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGATACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCTATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                                               <  2:877478/149‑1 (MQ=255)
                                                                                 cATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                              <  2:628665/149‑1 (MQ=255)
                                                                                 cATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAACATAAAACCTAGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                              >  2:305396/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCCTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCACCTGCGCa                                                             >  1:284529/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGACACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGACTATTCTCCAACTGCTCa                                                             >  1:222193/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                                             >  1:450051/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                                             >  2:839943/1‑149 (MQ=255)
                                                                                   tCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                                                            >  1:904561/1‑149 (MQ=255)
                                                                                   tCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                                                            >  1:904575/1‑149 (MQ=255)
                                                                                      ggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGCATATTCTCCAACTGCTCATAtt                                                         >  2:511932/1‑149 (MQ=255)
                                                                                            tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggcg                                                   >  1:1403592/1‑149 (MQ=255)
                                                                                                cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                               <  2:437538/149‑1 (MQ=255)
                                                                                                 ccccATGTGGTGGTACAGTAAGTGTTAAGCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCCCGGGAATCCCAGCACCATAAATTCTCGATCGTTCCATTCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                              <  1:1224961/149‑1 (MQ=255)
                                                                                                  ccccTGTGGTGGTACAGCAAGTTCTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                             <  1:395530/149‑1 (MQ=255)
                                                                                                  cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGG‑TCCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                             <  2:507825/149‑1 (MQ=255)
                                                                                                          gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                                     >  1:21890/1‑149 (MQ=255)
                                                                                                           tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGACTATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTg                                    >  2:912025/1‑149 (MQ=255)
                                                                                                           tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAGGCCACACTTTg                                    >  1:388590/1‑149 (MQ=255)
                                                                                                           tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTg                                    <  1:836699/149‑1 (MQ=255)
                                                                                                            ggtACCGTAAGTTGTACTCTCTCCATCGCCGCCCGGTTC‑TGTTCCAGAAACTGGGTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                   <  2:1370087/149‑1 (MQ=255)
                                                                                                            ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGGTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                   <  1:479011/149‑1 (MQ=255)
                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  <  1:920664/149‑1 (MQ=255)
                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATAATCTCCAACTGCTCATATTGCGGCGTTCAGACACACTTTGAg                                  >  2:1130541/1‑149 (MQ=255)
                                                                                                                 aGTAAGTTTTAATCCCTCCCTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                              <  1:198115/149‑1 (MQ=255)
                                                                                                                  gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTTCGGCGTTACGCCACACTTTCAGGCGAc                             >  1:487566/1‑149 (MQ=255)
                                                                                                                    aaGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACAt                           <  2:1268067/149‑1 (MQ=255)
                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTgc                        >  2:832827/1‑148 (MQ=255)
                                                                                                                        tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGGATCCCAGCGCGATAAAAACTCGATCGACCCAGGCTTGCTGAATATTCGCCACCTCCTCATATTGCGGCATTAAGCCACACTTTGAGGCGCCATTCAc                       >  2:376468/1‑149 (MQ=255)
                                                                                                                        tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc                       >  1:1302680/1‑149 (MQ=255)
                                                                                                                              ctctGCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                 <  1:95976/149‑1 (MQ=255)
                                                                                                                                ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               >  1:122453/1‑149 (MQ=255)
                                                                                                                                ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               <  1:1102866/149‑1 (MQ=255)
                                                                                                                                 ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGCTCGACCCAGGCTGTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg              >  1:202607/1‑149 (MQ=255)
                                                                                                                                 ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg              >  2:1447686/1‑149 (MQ=255)
                                                                                                                                             cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc  <  2:1039504/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793750‑1794040

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