Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,570,548 G→C 7.3% P106R (CCG→CGG)  sad‑II ← NAD+‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,570,5480GC 7.3% 210.1 / 5.8 68P106R (CCG→CGG) sad‑IINAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (24/39);  new base C (2/3);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGTGTTACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAGG  >  NC_002947/3570400‑3570694
                                                                                                                                                    |                                                                                                                                                  
ggTGTTACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccg                                                                                                                                                    >  8:249269/1‑149 (MQ=255)
    ttACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCCCGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccgggcc                                                                                                                                                >  6:157234/1‑149 (MQ=255)
        ccGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGc                                                                                                                                            >  2:19893/1‑149 (MQ=255)
         cGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCt                                                                                                                                           >  2:13788/1‑149 (MQ=255)
          gCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGCTACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTc                                                                                                                                          >  1:196399/1‑149 (MQ=255)
          gCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCACGGCCCGGCCATGCTc                                                                                                                                          <  8:300365/149‑1 (MQ=255)
          gCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCACGGCCCGGCCATGCTc                                                                                                                                          <  5:84809/149‑1 (MQ=255)
          gCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCACGGCCCGGCCATGCTc                                                                                                                                          >  6:84809/1‑149 (MQ=255)
           ccAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCg                                                                                                                                         <  3:62812/149‑1 (MQ=255)
           ccAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTCTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCg                                                                                                                                         <  1:245637/149‑1 (MQ=255)
            cAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCgg                                                                                                                                        <  2:262784/149‑1 (MQ=255)
              gAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCCGGATCGGGCCCCGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCg                                                                                                                                      >  3:89028/1‑149 (MQ=255)
                aTCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGt                                                                                                                                     <  4:128192/148‑1 (MQ=255)
                 tCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTAc                                                                                                                                   <  3:242142/149‑1 (MQ=255)
                        cACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccggg                                                                                                                                                  <  3:62403/127‑1 (MQ=255)
                        cACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATggccggg                                                                                                                                                  >  4:62403/1‑127 (MQ=255)
                         aCGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGc                                                                                                                            >  1:52566/1‑148 (MQ=255)
                            gcgcCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAg                                                                                                                        <  7:249269/149‑1 (MQ=255)
                                  cgcAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTGGGCGTACCAGCGGCAAAGGCCTac                                                                                                                  <  7:48727/149‑1 (MQ=255)
                                      gcaCCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTacacac                                                                                                              <  2:142675/149‑1 (MQ=255)
                                          cTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGc                                                                                                          <  2:52566/149‑1 (MQ=255)
                                             ccAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                       <  6:49926/149‑1 (MQ=255)
                                             ccAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGa                                                                                                       <  3:311550/149‑1 (MQ=255)
                                             ccAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGGCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCt                                                                                                         <  1:266503/147‑1 (MQ=255)
                                                  ccGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCa                                                                                                  <  1:13788/149‑1 (MQ=255)
                                                   cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTctc                                                                                                 >  1:321787/1‑147 (MQ=255)
                                                   cGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCAc                                                                                                 >  8:250461/1‑149 (MQ=255)
                                                      gAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACcgc                                                                                              <  1:19893/149‑1 (MQ=255)
                                                        aGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGTCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGa                                                                                            <  2:215836/149‑1 (MQ=255)
                                                          ttCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGc                                                                                          <  8:172535/149‑1 (MQ=255)
                                                            ccACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTacacac                                                                                                              <  7:69327/127‑1 (MQ=255)
                                                            ccACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTacacac                                                                                                              >  8:69327/1‑127 (MQ=255)
                                                             cACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTg                                                                                       >  6:166105/1‑149 (MQ=255)
                                                             cACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTg                                                                                       >  4:302443/1‑149 (MQ=255)
                                                              aCGGCATCACGGCTAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTgg                                                                                      >  5:69788/1‑149 (MQ=255)
                                                                ggCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGc                                                                                    <  1:135597/149‑1 (MQ=255)
                                                                         gCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGTCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCtt                                                                           <  3:222623/149‑1 (MQ=255)
                                                                          ccAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCTCCGCGAGCCTGGGCGATCGGCTTg                                                                          <  2:321787/149‑1 (MQ=255)
                                                                           cAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGc                                                                         >  3:152232/1‑149 (MQ=255)
                                                                           cAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGc                                                                         <  5:84321/149‑1 (MQ=255)
                                                                                     cccAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTcgcg                                                               <  3:302443/149‑1 (MQ=255)
                                                                                     cccAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTcgcg                                                               <  7:250461/149‑1 (MQ=255)
                                                                                        aGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGACAATTTCGCGGCt                                                            >  3:172634/1‑149 (MQ=255)
                                                                                            ggcggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGAt                                                         >  4:54210/1‑148 (MQ=255)
                                                                                              cctgtggcGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCAAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCAt                                                      <  8:29971/142‑1 (MQ=255)
                                                                                               ggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCtt                                                                           <  1:320040/127‑1 (MQ=255)
                                                                                               ggTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCtt                                                                           >  2:320040/1‑127 (MQ=255)
                                                                                                   ttCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTgc                                                  <  4:89028/148‑1 (MQ=255)
                                                                                                      gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGcgcg                                               >  7:179748/1‑148 (MQ=255)
                                                                                                      gATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGCCTGATCATGTGcgcg                                               >  7:247380/1‑148 (MQ=255)
                                                                                                        tACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTgcgc                                                <  8:227495/145‑1 (MQ=255)
                                                                                                        tACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTgcgc                                                >  7:227495/1‑145 (MQ=255)
                                                                                                        tACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGaaa                                            <  4:16815/149‑1 (MQ=255)
                                                                                                           gggCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCt                                         >  5:267992/1‑149 (MQ=255)
                                                                                                             gCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTc                                       <  4:172634/149‑1 (MQ=255)
                                                                                                              cTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCa                                      <  7:210057/149‑1 (MQ=255)
                                                                                                                     gACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTgg                               >  3:217761/1‑149 (MQ=255)
                                                                                                                     gACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTg                                <  1:52774/148‑1 (MQ=255)
                                                                                                                     gACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGCCAAAGGCCTACACACTTGCTGACTTCACCGCCAGCCTGGGCGTTCGGCTTGCCAATTTCGCGGCCGATCATGTGCGCGAAAGCTTCAGCCTTg                                <  1:20263/148‑1 (MQ=255)
                                                                                                                           ggTCGGCTGGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAATGGCCTACACACTTGCTGAGTTCACCGCGAGCCTTGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGa                         <  5:157234/149‑1 (MQ=255)
                                                                                                                           ggTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGa                         <  4:152232/149‑1 (MQ=255)
                                                                                                                            gTCGGCTCGGCGGCGAGGATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTGAGCCTTGGCTTCGAg                        <  3:54210/149‑1 (MQ=255)
                                                                                                                               ggctcggcGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGgt                     <  6:267992/149‑1 (MQ=255)
                                                                                                                                 ctcggcGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCtt                                 <  2:195874/135‑1 (MQ=255)
                                                                                                                                 ctcggcGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCtt                                 >  1:195874/1‑135 (MQ=255)
                                                                                                                                    ggcggcGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCg                 <  4:119003/148‑1 (MQ=255)
                                                                                                                                         cgAGCACGGCCCGGCCATGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagg           >  6:266561/1‑149 (MQ=255)
                                                                                                                                           aGCACGGCCCGGCCATGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggc         <  5:266561/149‑1 (MQ=255)
                                                                                                                                            gCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggcc        >  8:223550/1‑149 (MQ=255)
                                                                                                                                            gCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACAGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggc         >  4:247198/1‑148 (MQ=255)
                                                                                                                                              aTGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggccag      <  6:69788/149‑1 (MQ=255)
                                                                                                                                              aTGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCggccagggccag      >  3:29477/1‑149 (MQ=255)
                                                                                                                                                  ccgggccgTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAgg  <  2:196399/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGTTACCCGCCAGAATCGCCGGCACGGCGCCGCGCAGCACCTGCCAGACCGGGAAGTTCCACGGCATCACGGCCAGGATCGGGCCCAGCGGGCGGTATTCGATACGGGCTTTTTCGACCTGGGTCGGCTCGGCGGCGAGCATGGCCGGGCCGTGCTCGGCGTACCAGCGGCAAAGGCCTACACACTTGCTGACTTCACCGCGAGCCTGGGCGATCGGCTTGCCAATTTCGCGGCTGATCATGTGCGCGAAAGCTTCAGCCTTGGCTTCGAGGGTGTCGGCCAGGGCCAGCAGG  >  NC_002947/3570400‑3570694

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