Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,075,435 C→T 5.7% A1150T (GCA→ACA)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,075,4350CT 5.7% 226.6 / 3.7 70A1150T (GCA→ACA) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base T (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

TCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGC  >  NC_000913/1075287‑1075581
                                                                                                                                                    |                                                                                                                                                  
tcatcGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTtgctgc                                                                                                                                                    >  3:420367/1‑149 (MQ=255)
   tcGTCAATACACAACACTCGCTCACGTGCCCTCACGGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCACCACTCGTTGTATTCCTTTCTGCGCCACCTCGCCACCTGGCGTACATGACGCGTTCCATCCTGGACCATTTGCCGCCCa                                                                                                                                                 >  4:475934/1‑149 (MQ=255)
      gCAATACACAACACGCGCTCACACGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTc                                                                                                                                              >  2:692945/1‑149 (MQ=255)
         aTACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCg                                                                                                                                           >  4:50548/1‑149 (MQ=255)
          tACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTc                                                                                                                                              >  3:580578/1‑145 (MQ=255)
          tACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTc                                                                                                                                              <  4:580578/145‑1 (MQ=255)
             acaacaCGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGt                                                                                                                                       >  4:666975/1‑149 (MQ=255)
             acaacaCGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGt                                                                                                                                       <  2:721570/149‑1 (MQ=255)
                 cacGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCAt                                                                                                                                   >  1:412205/1‑149 (MQ=255)
                        tCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCgg                                                                                                                            <  1:284587/149‑1 (MQ=255)
                        tCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCgg                                                                                                                            >  1:305684/1‑149 (MQ=255)
                          aCGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCt                                                                                                                          >  2:792551/1‑149 (MQ=255)
                            gcgGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGa                                                                                                                        >  2:256487/1‑149 (MQ=255)
                            gcgGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGa                                                                                                                        <  4:616527/149‑1 (MQ=255)
                            gcgGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGa                                                                                                                        <  1:522671/149‑1 (MQ=255)
                                    agcGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGc                                                                                                                >  4:589323/1‑149 (MQ=255)
                                    agcGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGc                                                                                                                <  1:574505/149‑1 (MQ=255)
                                    agcGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGc                                                                                                                <  1:464093/149‑1 (MQ=255)
                                    agcGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAAGGc                                                                                                                <  1:656678/149‑1 (MQ=255)
                                       gTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGc                                                                                                             <  3:475934/149‑1 (MQ=255)
                                       gTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGc                                                                                                             <  3:666975/149‑1 (MQ=255)
                                        tCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCt                                                                                                            <  1:431492/149‑1 (MQ=255)
                                           aGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTc                                                                                                         <  3:50548/149‑1 (MQ=255)
                                           aGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTc                                                                                                         >  4:198141/1‑149 (MQ=255)
                                            ggTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCa                                                                                                        >  4:599173/1‑149 (MQ=255)
                                             gtgtTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCaa                                                                                                       >  2:114111/1‑149 (MQ=255)
                                                  gCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTgcg                                                                                                  >  3:387551/1‑149 (MQ=255)
                                                  gCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTgcg                                                                                                  <  1:50647/149‑1 (MQ=255)
                                                      tCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATc                                                                                              <  4:630269/149‑1 (MQ=255)
                                                        aCCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGTTGCCCATTCCCGCAGt                                                                                                                                       <  1:713788/106‑1 (MQ=255)
                                                        aCCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGTTGCCCATTCCCGCAGt                                                                                                                                       >  2:713788/1‑106 (MQ=255)
                                                          ccGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGAcc                                                                                          >  1:640701/1‑149 (MQ=255)
                                                             tCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGa                                                                                       <  3:599173/149‑1 (MQ=255)
                                                                     ggCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACttt                                                                                 >  4:247168/1‑147 (MQ=255)
                                                                     ggCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACttt                                                                                 <  3:247168/147‑1 (MQ=255)
                                                                           aaTCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACtt                                                                                  >  3:75348/1‑140 (MQ=255)
                                                                           aaTCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACtt                                                                                  <  4:75348/140‑1 (MQ=255)
                                                                             tCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCCGAc                                                                       <  1:692945/149‑1 (MQ=255)
                                                                               gTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGAcgc                                                                     >  2:236004/1‑149 (MQ=255)
                                                                                ttGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGAcgcg                                                                    >  1:807652/1‑149 (MQ=255)
                                                                                ttGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGAcgcg                                                                    >  4:76452/1‑149 (MQ=255)
                                                                                     ttCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGc                                                               <  1:792551/149‑1 (MQ=255)
                                                                                      tCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGTCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTc                                                                                                         >  4:91902/1‑106 (MQ=255)
                                                                                      tCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGTCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTc                                                                                                         <  3:91902/106‑1 (MQ=255)
                                                                                          gcctgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCAc                                                          >  2:574376/1‑149 (MQ=255)
                                                                                             tgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGc                                                       <  1:114111/149‑1 (MQ=255)
                                                                                               cgcCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCa                                                     >  2:812636/1‑149 (MQ=255)
                                                                                                     cTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCAc                                               >  2:556185/1‑149 (MQ=255)
                                                                                                        gCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACttt                                            >  2:754653/1‑149 (MQ=255)
                                                                                                                gCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGc                                                                                                                              >  3:104702/1‑59 (MQ=255)
                                                                                                                gCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGc                                                                                                                              <  4:104702/59‑1 (MQ=255)
                                                                                                                gCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGa                                    <  2:640701/149‑1 (MQ=255)
                                                                                                                  gTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGAtt                                  <  3:589323/149‑1 (MQ=255)
                                                                                                                     cATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGc                               >  1:471096/1‑149 (MQ=255)
                                                                                                                      aTAACGCCTGCAATTCTGGACGATTTGCTGTCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAg                                                                >  1:534367/1‑115 (MQ=255)
                                                                                                                      aTAACGCCTGCAATTCTGGACGATTTGCTGTCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAg                                                                <  2:534367/115‑1 (MQ=255)
                                                                                                                      aTAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGcc                              >  4:505491/1‑149 (MQ=255)
                                                                                                                        aaCGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCcag                            >  1:415599/1‑149 (MQ=255)
                                                                                                                          cGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCcagca                          <  1:236004/149‑1 (MQ=255)
                                                                                                                          cGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCcagca                          <  4:514478/149‑1 (MQ=255)
                                                                                                                          cGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCcagca                          <  3:233393/149‑1 (MQ=255)
                                                                                                                             cTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGAc                       >  4:453693/1‑149 (MQ=255)
                                                                                                                                   ttCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTc                                 >  2:4081/1‑133 (MQ=255)
                                                                                                                                   ttCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTc                                 >  2:246928/1‑133 (MQ=255)
                                                                                                                                   ttCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTc                                 <  1:246928/133‑1 (MQ=255)
                                                                                                                                   ttCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTc                                 <  1:4081/133‑1 (MQ=255)
                                                                                                                                     cTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAgaga               <  3:76452/149‑1 (MQ=255)
                                                                                                                                      tGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTagag                >  3:157533/1‑147 (MQ=255)
                                                                                                                                      tGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTagag                <  4:157533/147‑1 (MQ=255)
                                                                                                                                       ggACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATa             <  4:336065/149‑1 (MQ=255)
                                                                                                                                           gATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGc         >  3:759888/1‑149 (MQ=255)
                                                                                                                                                  tgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGc  >  3:596851/1‑149 (MQ=255)
                                                                                                                                                  tgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGc  >  4:595695/1‑149 (MQ=255)
                                                                                                                                                   gcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCTCGCGCGTc                                                            <  1:418074/90‑1 (MQ=255)
                                                                                                                                                   gcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCTCGCGCGTc                                                            >  2:418074/1‑90 (MQ=255)
                                                                                                                                                    cccATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCa                                                >  3:778786/1‑101 (MQ=255)
                                                                                                                                                    cccATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCa                                                <  4:778786/101‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGC  >  NC_000913/1075287‑1075581

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