Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,803 C→T 6.9% Q33Q (CAG→CAA ygfF ← putative NAD(P)‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,8030CT 6.9% 277.7 / 7.0 89Q33Q (CAG→CAAygfFputative NAD(P)‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base T (3/3);  total (44/45)
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.72e-01

CTGCTAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTG  >  NC_000913/3045658‑3045947
                                                                                                                                                 |                                                                                                                                                
ctgctAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGa                                                                                                                                               <  4:370020/149‑1 (MQ=255)
   ctAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCTCCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATaa                                                                                                                                            <  3:393646/149‑1 (MQ=255)
   ctAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATaa                                                                                                                                            >  1:175988/1‑149 (MQ=255)
   ctAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATAAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATaa                                                                                                                                            <  8:322148/149‑1 (MQ=255)
         gTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAAcc                                                                                                                                      <  8:429242/149‑1 (MQ=255)
             aTCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCa                                                                                                                                  >  3:196132/1‑149 (MQ=255)
                    tGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGtata                                                                                                                           >  6:204398/1‑149 (MQ=255)
                     gATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATAc                                                                                                                          >  3:192258/1‑149 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAcc                                                                                                                         <  2:2885/149‑1 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAcc                                                                                                                         >  3:233214/1‑149 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAcc                                                                                                                         >  5:25301/1‑149 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAcc                                                                                                                         >  6:108373/1‑149 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGGTTTGCTGATAATTAACCGCCACCGtata                                                                                                                           <  7:206651/147‑1 (MQ=255)
                      aTCGATTGCTGTAAACATCGCAACGACCTGGTATTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAAcc                                                                                                                                      <  6:286025/136‑1 (MQ=255)
                       tCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCATTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAccc                                                                                                                        >  3:128394/1‑149 (MQ=255)
                              ccttAAACATCGCACCGACCTGGTATTCGGCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCAGCACTTCCTTCGCCGCGTTGAGGTTTTGCTGATAATTAACCTCCAACGTATACCCTTCTTg                                                                                                                  <  6:285576/145‑1 (MQ=255)
                              cTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGtata                                                                                                                           <  8:111414/139‑1 (MQ=255)
                              cTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGtata                                                                                                                           >  7:111414/1‑139 (MQ=255)
                              cTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTgc                                                                                                                 >  3:306169/1‑149 (MQ=255)
                                 tAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGTCTTGCGTTATTAAGTTCATAACTTCCTGCGCCGCGTGGAGGTTTTGCTTATAATTAACCGACACCGTATCCACTTCTTgcgc                                                                                                               >  6:403307/1‑148 (MQ=255)
                                 tAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTgcgc                                                                                                               >  8:231110/1‑148 (MQ=255)
                                   aaCATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCa                                                                                                             <  4:250282/148‑1 (MQ=255)
                                      aTCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         >  1:70887/1‑149 (MQ=255)
                                      aTCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         <  1:137982/149‑1 (MQ=255)
                                       tCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAccc                                                                                                                        >  3:218295/1‑133 (MQ=255)
                                       tCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATAccc                                                                                                                        <  4:218295/133‑1 (MQ=255)
                                         gCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTaa                                                                                                      <  2:70887/149‑1 (MQ=255)
                                         gCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTaa                                                                                                      <  5:194194/149‑1 (MQ=255)
                                         gCAACGACCTGGTTTTCGTAGCTGATATCCGCCTGGAGCAAGAATGCTTTGCCACCGGCTTGAGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTTATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTa                                                                                                       <  2:404304/148‑1 (MQ=255)
                                           aaCGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATg                                                                                                    <  4:192258/149‑1 (MQ=255)
                                           aaCGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTACCCGCCACCGGATACCCTTCTTGCGCCAACAGTAATg                                                                                                    >  3:47240/1‑149 (MQ=255)
                                              gACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAg                                                                                                 <  5:128541/149‑1 (MQ=255)
                                               aCCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATaa                                                                                                                                            <  4:248961/105‑1 (MQ=255)
                                               aCCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATaa                                                                                                                                            >  3:248961/1‑105 (MQ=255)
                                               aCCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGt                                                                                                <  4:233214/149‑1 (MQ=255)
                                                ccTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGAGTGATTAAGTTCCTCACTTCCTGCGCCGCGTGGAGGTTTTGGTGATAATTACAAGACCCCGTATACCCTTCTTGCGCCAACAGTAATGCAggt                                                                                               >  6:404077/1‑147 (MQ=255)
                                                 cTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         >  8:272280/1‑138 (MQ=255)
                                                 cTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGTTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         <  7:272280/138‑1 (MQ=255)
                                                 cTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         >  6:137013/1‑138 (MQ=255)
                                                 cTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAg                                                                                                         <  5:137013/138‑1 (MQ=255)
                                                  tGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTgc                                                                                             >  5:200653/1‑149 (MQ=255)
                                                   ggTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTgcc                                                                                            >  8:319383/1‑149 (MQ=255)
                                                    gTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCTTATACCCTTCTTGCGCCAACAGTAATGCAGTTgccc                                                                                           <  4:306169/149‑1 (MQ=255)
                                                    gTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTgc                                                                                             <  8:237161/147‑1 (MQ=255)
                                                     ttttCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTgcccg                                                                                          >  5:203160/1‑149 (MQ=255)
                                                          gtcgCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGa                                                                                     <  8:118731/149‑1 (MQ=255)
                                                             gCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGAt                                                                                                                                              <  6:333222/89‑1 (MQ=255)
                                                             gCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGAt                                                                                                                                              >  5:333222/1‑89 (MQ=255)
                                                               tGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTACCCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCcg                                                                                >  6:413679/1‑149 (MQ=255)
                                                               tGATATCCGCCTGGAGCACGAATGCTTAGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGtt                                                                                               >  7:190098/1‑134 (MQ=255)
                                                                 atatCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCcgcg                                                                              <  7:231110/149‑1 (MQ=255)
                                                                  tatCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGa                                                                             <  8:235087/149‑1 (MQ=255)
                                                                      cGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACt                                                                           <  1:387665/147‑1 (MQ=255)
                                                                      cGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACt                                                                           >  2:387665/1‑147 (MQ=255)
                                                                       gCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGc                                                                         >  2:270729/1‑148 (MQ=255)
                                                                            gAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGccacca                                                                    <  2:175988/148‑1 (MQ=255)
                                                                             aGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGt                                                                  <  6:203160/149‑1 (MQ=255)
                                                                                aCGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTca                                                                >  8:150326/1‑148 (MQ=255)
                                                                                        tttGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCta                                                       >  6:168399/1‑149 (MQ=255)
                                                                                        tttGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCta                                                       >  1:67878/1‑149 (MQ=255)
                                                                                            ccACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGc                                                   <  7:319383/149‑1 (MQ=255)
                                                                                             cACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGc                                                   >  5:236538/1‑148 (MQ=255)
                                                                                             cACCGGCTTGCGTGATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCctt                                                                                                                      >  8:357343/1‑81 (MQ=255)
                                                                                              aCCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCa                                                 >  5:69015/1‑149 (MQ=255)
                                                                                                 ggCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCAttt                                              <  5:168399/149‑1 (MQ=255)
                                                                                                   cTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGt                                            <  6:25301/149‑1 (MQ=255)
                                                                                                       cGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGa                                                                             >  1:296852/1‑112 (MQ=255)
                                                                                                       cGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGa                                                                             <  2:296852/112‑1 (MQ=255)
                                                                                                          tattaAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCTGTTGCCCGCCCGATGCCGCGACTTccacca                                                                    <  7:71896/118‑1 (MQ=255)
                                                                                                          tattaAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTcc                                      >  7:356447/1‑148 (MQ=255)
                                                                                                          tattaAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTcc                                      >  3:119505/1‑148 (MQ=255)
                                                                                                          tattaAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGACACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCACCATGCCGCGACTGccacca                                                                    >  8:71896/1‑118 (MQ=255)
                                                                                                                gTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGaaa                                >  5:95549/1‑148 (MQ=255)
                                                                                                                     tcaCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaa                          >  4:40459/1‑149 (MQ=255)
                                                                                                                     tcaCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaa                          <  4:196132/149‑1 (MQ=255)
                                                                                                                     tcaCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaa                          >  7:59287/1‑149 (MQ=255)
                                                                                                                     tcaCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaa                          >  5:176934/1‑149 (MQ=255)
                                                                                                                      caCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaat                         >  2:388487/1‑149 (MQ=255)
                                                                                                                       aCTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAaattaat                         <  2:50602/148‑1 (MQ=255)
                                                                                                                                  cgcgTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTAt             <  6:69015/149‑1 (MQ=255)
                                                                                                                                      tGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGt                                                                                                        >  7:60241/1‑54 (MQ=255)
                                                                                                                                      tGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGt                                                                                                        <  8:60241/54‑1 (MQ=255)
                                                                                                                                      tGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCttt         <  2:67878/149‑1 (MQ=255)
                                                                                                                                      tGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCttt         <  7:150326/149‑1 (MQ=255)
                                                                                                                                        gAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGc                                                                                  >  8:141193/1‑74 (MQ=255)
                                                                                                                                        gAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGc                                                                                  <  7:141193/74‑1 (MQ=255)
                                                                                                                                        gAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTga       <  8:356447/149‑1 (MQ=255)
                                                                                                                                        gAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTg        <  1:270729/148‑1 (MQ=255)
                                                                                                                                              ttGCTGCTAATTAACCGCCACCGTATACCCTGCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCAGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTg  <  5:204398/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGCTAGCGGTTCATCGTGCTGATCGATTGCTGTAAACATCGCAACGACCTGGTTTTCGTCGCTGATATCCGCCTGGAGCACGAATGCTTTGCCACCGGCTTGCGTTATTAAGTTCATCACTTCCTGCGCCGCGTGGAGGTTTTGCTGATAATTAACCGCCACCGTATACCCTTCTTGCGCCAACAGTAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTG  >  NC_000913/3045658‑3045947

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