Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,078,235 A→G 5.9% T216T (ACT→ACC)  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,078,2350AG 5.9% 260.7 / 3.7 85T216T (ACT→ACC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base G (2/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGA  >  NC_000913/1078087‑1078382
                                                                                                                                                    |                                                                                                                                                   
cgATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCa                                                                                                                                                     <  6:11529/149‑1 (MQ=255)
cgATGGGTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCAATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCa                                                                                                                                                     <  6:95426/149‑1 (MQ=255)
cgATGGAGTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCa                                                                                                                                                     <  8:450078/149‑1 (MQ=255)
          gCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAAcag                                                                                                                                             >  1:124754/1‑147 (MQ=255)
          gCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGGGGAAACCAGTTTGCCAGTAAAcagca                                                                                                                                           <  2:124754/149‑1 (MQ=255)
          gCCAGTGACGAACTGCGCCCCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAAcagca                                                                                                                                           <  3:237051/149‑1 (MQ=255)
           ccAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAAcagcag                                                                                                                                          <  5:19389/149‑1 (MQ=255)
            cAGTGAAGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACagcagc                                                                                                                                         >  1:272590/1‑149 (MQ=255)
                gACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCa                                                                                                                                                     >  6:104121/1‑133 (MQ=255)
                gACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCa                                                                                                                                                     <  5:104121/133‑1 (MQ=255)
                  cGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAg                                                                                                                                   <  4:157565/149‑1 (MQ=255)
                       tGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTggc                                                                                                                              >  4:294385/1‑149 (MQ=255)
                         cTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTggcgg                                                                                                                            >  8:363433/1‑149 (MQ=255)
                          tCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCTCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTggcgg                                                                                                                            >  4:388258/1‑148 (MQ=255)
                          tCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGcc                                                                                                                                        <  7:178759/136‑1 (MQ=255)
                          tCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGcc                                                                                                                                        >  8:178759/1‑136 (MQ=255)
                          tCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTggcgg                                                                                                                            >  7:150037/1‑148 (MQ=255)
                           cACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCa                                                                                                                          >  2:251072/1‑149 (MQ=255)
                            aCCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCGGTAAACAGCAGCCCCCAGGTGGCGGCAt                                                                                                                         >  7:184396/1‑149 (MQ=255)
                              ccATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCAtt                                                                                                                        >  3:89068/1‑148 (MQ=255)
                              ccATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCAtt                                                                                                                        >  3:162003/1‑148 (MQ=255)
                              ccATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCAt                                                                                                                         <  2:269809/147‑1 (MQ=255)
                              ccATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCAt                                                                                                                         >  1:269809/1‑147 (MQ=255)
                                 tCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCGGTAAACAGCAGCCCCCAGGTGGCGGCATTaaca                                                                                                                    <  6:319460/149‑1 (MQ=255)
                                 tCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTaac                                                                                                                     >  4:246639/1‑148 (MQ=255)
                                  cAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAAcagca                                                                                                                                           >  1:457781/1‑125 (MQ=255)
                                  cAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAAcagca                                                                                                                                           <  2:457781/125‑1 (MQ=255)
                                   aGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGATCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCATTTTGCCAGTAAACAGCAGCCCCCAGGTGGAGGCATTaacaaa                                                                                                                  <  2:182100/149‑1 (MQ=255)
                                   aGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGGTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTaacaa                                                                                                                   <  7:363433/148‑1 (MQ=255)
                                    ggCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGTGGGGCTGGCTTCGGTATGGGGGGAAACGAGTTTGGCAGTAAACAGCAGCCCCCAGGGGGCGGCGTTagcagcc                                                                                                                 >  8:353851/1‑146 (MQ=255)
                                                 atCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAACAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTAc                                                                                                    >  4:235465/1‑149 (MQ=255)
                                                   ccACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAACCCAGTTTGCCCGTAAATAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGa                                                                                                  >  1:433132/1‑149 (MQ=255)
                                                   ccACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGa                                                                                                  >  2:273509/1‑149 (MQ=255)
                                                    cacaCCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGa                                                                                                  <  4:89068/148‑1 (MQ=255)
                                                    cacaCCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGAc                                                                                                 <  5:204799/149‑1 (MQ=255)
                                                     acacCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGAc                                                                                                 >  4:384899/1‑148 (MQ=255)
                                                     acacCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGAc                                                                                                 >  2:298884/1‑148 (MQ=255)
                                                       acCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCaa                                                                                              >  6:208423/1‑149 (MQ=255)
                                                       acCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCaa                                                                                              <  3:297413/149‑1 (MQ=255)
                                                         cTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGTAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAAt                                                                                             <  2:237900/148‑1 (MQ=255)
                                                         cTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGTAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAAt                                                                                             >  1:237900/1‑148 (MQ=255)
                                                          tttGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAAtgt                                                                                           >  1:321008/1‑149 (MQ=255)
                                                            tGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCGGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAAtgtg                                                                                          <  8:184396/148‑1 (MQ=255)
                                                            tGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAAtgt                                                                                           <  3:246639/147‑1 (MQ=255)
                                                              cGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGa                                                                                       >  8:135620/1‑149 (MQ=255)
                                                                    aGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGTGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGTTGTCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGc                                                                                   <  1:349349/147‑1 (MQ=255)
                                                                    aGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGCCGGGCTACGACCAATGTGTGACTGc                                                                                   >  6:78138/1‑147 (MQ=255)
                                                                     gCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAg                                                                                >  4:269124/1‑149 (MQ=255)
                                                                           cACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCACCAGGTGGCGGCATTaacaaaca                                                                                                                >  2:88739/1‑111 (MQ=255)
                                                                             ccGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGtt                                                                        >  1:18859/1‑149 (MQ=255)
                                                                                   ttACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGa                                                                   <  2:321008/148‑1 (MQ=255)
                                                                                     aCCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGAtt                                                                 <  4:162003/148‑1 (MQ=255)
                                                                                      ccGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGAttt                                                                <  2:6259/148‑1 (MQ=255)
                                                                                       cGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTg                                                              >  7:154580/1‑149 (MQ=255)
                                                                                          tAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTcg                                                           <  1:251072/149‑1 (MQ=255)
                                                                                            aTGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTcgc                                                          <  7:157530/148‑1 (MQ=255)
                                                                                             tGCGGTTCCGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTcgcg                                                         >  6:144726/1‑148 (MQ=255)
                                                                                               cGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAAt                                                      <  7:135620/149‑1 (MQ=255)
                                                                                                     agcgagcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAAcgc                                                <  3:294385/149‑1 (MQ=255)
                                                                                                     agcgagcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAAcg                                                 <  1:298884/148‑1 (MQ=255)
                                                                                                          gcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCGGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAAGAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTTGCGAATTAACGCGTc                                             <  6:193232/147‑1 (MQ=255)
                                                                                                          gcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCGGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTTGCGAATTAACGCGTc                                             >  5:193232/1‑147 (MQ=255)
                                                                                                            gggAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCgg                                         <  3:384899/149‑1 (MQ=255)
                                                                                                                   gctggctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTcgc                                           <  3:111337/140‑1 (MQ=255)
                                                                                                                   gctggctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTcgc                                           >  4:111337/1‑140 (MQ=255)
                                                                                                                       gctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTCGCTTTGTc                              >  4:197481/1‑149 (MQ=255)
                                                                                                                        ctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGGTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCg                             <  2:272590/149‑1 (MQ=255)
                                                                                                                         ttCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            >  3:398932/1‑149 (MQ=255)
                                                                                                                         ttCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            >  4:217764/1‑149 (MQ=255)
                                                                                                                         ttCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            >  7:353866/1‑149 (MQ=255)
                                                                                                                         ttCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCAACCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAAAGCGTCGCGGGTGGCTTTGTCgg                            >  3:240284/1‑149 (MQ=255)
                                                                                                                           cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCggg                           <  8:353866/148‑1 (MQ=255)
                                                                                                                           cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            <  5:78138/147‑1 (MQ=255)
                                                                                                                           cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            <  3:217764/147‑1 (MQ=255)
                                                                                                                           cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            <  3:388258/147‑1 (MQ=255)
                                                                                                                           cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                            <  7:353851/147‑1 (MQ=255)
                                                                                                                             ttATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAAt                        <  8:150037/149‑1 (MQ=255)
                                                                                                                                  ggTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCttt                                 <  1:63946/135‑1 (MQ=255)
                                                                                                                                  ggTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCttt                                 >  2:63946/1‑135 (MQ=255)
                                                                                                                                   gTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaa                  >  6:107061/1‑149 (MQ=255)
                                                                                                                                   gTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaa                  >  7:159001/1‑149 (MQ=255)
                                                                                                                                   gTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaa                  <  5:107061/149‑1 (MQ=255)
                                                                                                                                   gTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTCGCGGCATTAACAAACAGGGACGGGCTACGACCAATGTGTGCCTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaa                  >  6:92098/1‑149 (MQ=255)
                                                                                                                                      gAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaacaa               >  7:346900/1‑149 (MQ=255)
                                                                                                                                      gAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGcaacaa               >  5:451764/1‑149 (MQ=255)
                                                                                                                                       aaaCCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGAGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCGTTGTCGGGAATACGCaacaac              <  2:18859/149‑1 (MQ=255)
                                                                                                                                        aaCCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACg             >  7:397329/1‑149 (MQ=255)
                                                                                                                                                   cAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGca                   >  5:432666/1‑132 (MQ=255)
                                                                                                                                                   cAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGca                   <  6:432666/132‑1 (MQ=255)
                                                                                                                                                   cAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGa  <  3:269124/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGA  >  NC_000913/1078087‑1078382

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: ‑