Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,673,113 G→A 5.1% A234T (GCC→ACC)  yhjC → LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,673,1130GA 5.1% 262.8 / 2.8 78A234T (GCC→ACC) yhjCLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

TTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGG  >  NC_000913/3672968‑3673256
                                                                                                                                                 |                                                                                                                                               
ttCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGcc                                                                                                                                               >  6:64640/1‑148 (MQ=255)
 tCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCg                                                                                                                                              >  3:30220/1‑148 (MQ=255)
      cGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTg                                                                                                                                         <  4:99065/148‑1 (MQ=255)
      cGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGGTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTg                                                                                                                                         <  3:12186/148‑1 (MQ=255)
        cAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTggg                                                                                                                                       >  4:10659/1‑148 (MQ=255)
        cAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTggg                                                                                                                                       <  4:30220/148‑1 (MQ=255)
        cAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTggg                                                                                                                                       <  2:37149/148‑1 (MQ=255)
          aTAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGa                                                                                                                                     >  3:124115/1‑148 (MQ=255)
             gTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGatt                                                                                                                                   <  1:115238/147‑1 (MQ=255)
               gCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGattatt                                                                                                                                >  1:8558/1‑148 (MQ=255)
                cGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGattatt                                                                                                                                <  6:127683/147‑1 (MQ=255)
                cGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGattatt                                                                                                                                <  4:124115/147‑1 (MQ=255)
                cGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGGGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTc                                                                                                                               <  3:92338/148‑1 (MQ=255)
                 gTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCa                                                                                                                              <  6:179875/148‑1 (MQ=255)
                    ataCACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGattatt                                                                                                                                >  8:205364/1‑143 (MQ=255)
                     taCACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGCCGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATc                                                                                                                          >  3:73926/1‑148 (MQ=255)
                       cacCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGGTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAACCTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATcc                                                                                                                         <  2:272365/147‑1 (MQ=255)
                          cGCCCCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCcgcg                                                                                                                     >  6:34490/1‑148 (MQ=255)
                           gccccTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCgcgc                                                                                                                    <  4:42710/145‑1 (MQ=255)
                           gCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCcgcg                                                                                                                     >  2:167266/1‑147 (MQ=255)
                           gCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCgcgc                                                                                                                    <  7:205364/148‑1 (MQ=255)
                            cACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAACCTATCTCACCGCCGGTCTTACCGGTCTGGGGATTATTCAGATCCCGCGCa                                                                                                                   >  1:76164/1‑148 (MQ=255)
                            cACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCa                                                                                                                   <  1:226628/148‑1 (MQ=255)
                                tGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGc                                                                                                               >  1:163377/1‑148 (MQ=255)
                                      tACATCCGTAAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTgcg                                                                                                         <  1:272291/148‑1 (MQ=255)
                                              ttAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTACCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGccc                                                                                                 >  3:260398/1‑148 (MQ=255)
                                                  gTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCg                                                                                             <  1:167266/148‑1 (MQ=255)
                                                    tttGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTg                                                                                           >  7:160720/1‑148 (MQ=255)
                                                          gTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGggc                                                                                     <  6:144710/148‑1 (MQ=255)
                                                                aGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTa                                                                               >  3:85288/1‑148 (MQ=255)
                                                                 gCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTAt                                                                              <  8:160720/148‑1 (MQ=255)
                                                                   gTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTACCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTg                                                                            <  2:76164/148‑1 (MQ=255)
                                                                    ttAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTATCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGa                                                                           <  2:57032/148‑1 (MQ=255)
                                                                    ttAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTTAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGa                                                                           <  5:64640/148‑1 (MQ=255)
                                                                    ttAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGa                                                                           <  2:8558/148‑1 (MQ=255)
                                                                       aTGGCGTCCAGTGGTTTAAGGCTGGCGGCATGTTGACGGGAAACAGTAGCGAAAACTATCTCACCGCCGGGCTTGCCGGTCGGGGGATTATTCAGATCCCGCTCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGt                                                                        <  5:34490/148‑1 (MQ=255)
                                                                        tGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTa                                                                       >  8:3149/1‑148 (MQ=255)
                                                                           cGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTAtt                                                                     >  4:196419/1‑147 (MQ=255)
                                                                           cGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTACCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTa                                                                    <  4:260398/148‑1 (MQ=255)
                                                                           cGTCCAGTGGGTTAAGTCTGGCGGCATGTTGACGGAAAACCGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCAATGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTa                                                                    <  5:34372/148‑1 (MQ=255)
                                                                             tCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTAcc                                                                  >  1:132635/1‑148 (MQ=255)
                                                                              ccAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTa                                                                    >  4:278079/1‑145 (MQ=255)
                                                                              ccAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTa                                                                    <  3:278079/145‑1 (MQ=255)
                                                                                 gTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTAcc                                                                  <  1:27522/144‑1 (MQ=255)
                                                                                 gTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTAcc                                                                  >  2:27522/1‑144 (MQ=255)
                                                                                  tGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTTTCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTa                                                                       >  5:152170/1‑138 (MQ=255)
                                                                                  tGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCggg                                                                                      <  8:98255/123‑1 (MQ=255)
                                                                                  tGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCggg                                                                                      >  7:98255/1‑123 (MQ=255)
                                                                                  tGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTa                                                                       <  6:152170/138‑1 (MQ=255)
                                                                                          gTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGcc                                                     >  6:15340/1‑148 (MQ=255)
                                                                                                     tGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCtctc                                          >  3:71810/1‑148 (MQ=255)
                                                                                                      gTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTcc                                         >  2:275578/1‑148 (MQ=255)
                                                                                                      gTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTcc                                         >  6:83667/1‑148 (MQ=255)
                                                                                                          aCGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCttt                                     >  1:220386/1‑148 (MQ=255)
                                                                                                          aCGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCttt                                     <  4:85288/148‑1 (MQ=255)
                                                                                                           cGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTc                                    >  6:195182/1‑148 (MQ=255)
                                                                                                              tAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCt                                 <  8:257376/148‑1 (MQ=255)
                                                                                                              tAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCt                                 >  6:195682/1‑148 (MQ=255)
                                                                                                                aaCAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTgg                               <  3:10659/148‑1 (MQ=255)
                                                                                                                  cAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGtt                             >  6:109683/1‑148 (MQ=255)
                                                                                                                        cGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCtt                                      >  7:147997/1‑133 (MQ=255)
                                                                                                                        cGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCtt                                      <  8:147997/133‑1 (MQ=255)
                                                                                                                        cGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCg                       >  8:158961/1‑148 (MQ=255)
                                                                                                                         gAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGc                      <  2:132635/148‑1 (MQ=255)
                                                                                                                           aaaCTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAg                    >  3:227677/1‑148 (MQ=255)
                                                                                                                           aaaCTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCACTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAg                    >  7:177012/1‑148 (MQ=255)
                                                                                                                            aaCTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGc                   >  6:304216/1‑148 (MQ=255)
                                                                                                                                aTCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCg               >  8:335381/1‑148 (MQ=255)
                                                                                                                                 tctcACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCg               <  5:83667/147‑1 (MQ=255)
                                                                                                                                 gctcACCGCCGGTCTTGCCGGGCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCGACCGTGCCGAGCCGCTCTCCCTTTCGCTGGGTTATCCGCAGCGTCg               <  1:235908/146‑1 (MQ=255)
                                                                                                                                  ctcACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCGACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCggg             >  3:245333/1‑148 (MQ=255)
                                                                                                                                     aCCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGc          >  8:167889/1‑148 (MQ=255)
                                                                                                                                        gccgGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCttt       >  1:109673/1‑148 (MQ=255)
                                                                                                                                        gccgGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCAGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATTCGCAGCGTCGGGAGCttt       <  2:163377/148‑1 (MQ=255)
                                                                                                                                          cgGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTcc     <  1:85664/148‑1 (MQ=255)
                                                                                                                                            gTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCg   >  7:54956/1‑148 (MQ=255)
                                                                                                                                             tCTTGCCGGTCGGGGGATTCTTCAGATCCCGCGCATGGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCgg  <  4:73926/148‑1 (MQ=255)
                                                                                                                                             tCTTGCCGGTATGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCgg  >  2:88965/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAAGTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGG  >  NC_000913/3672968‑3673256

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: