Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,745,036 C→T 14.9% Q108* (CAA→TAA) ‡ paoC → aromatic aldehyde dehydrogenase molybdopterin‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,745,0360CT14.9% 163.3 / 16.0 67Q108* (CAA→TAA) ‡paoCaromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (4/6);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CGTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACACCGCCGC  >  NC_002947/3744897‑3745183
                                                                                                                                           |                                                                                                                                                    
cgTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTgg                                                                                                                                             <  5:288201/149‑1 (MQ=255)
 gTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgta                                                                                                                                            <  3:360665/149‑3 (MQ=255)
    cccGGGGTGATCGCGGTGCTTGATCTTGGCCATCGTCCAGATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGGCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtaga                                                                                                                                          <  4:51020/148‑5 (MQ=255)
    cccGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtagag                                                                                                                                         <  1:84475/149‑6 (MQ=255)
    cccGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtagag                                                                                                                                         <  2:179605/149‑6 (MQ=255)
         ggTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtggtgg                                                                                                                                    <  1:92679/149‑1 (MQ=255)
           tGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGtt                                                                                                                                  <  4:114778/149‑1 (MQ=255)
                   gTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACTCAGACGCCTCGCAGGGTGCGCCGTTTAGGCCCTTGTTCGATTACCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                         <  1:88177/149‑1 (MQ=255)
                       ttGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtg                                                                                                                                        >  8:238010/1‑131 (MQ=255)
                       ttGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtg                                                                                                                                        <  7:238010/131‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTgg                                                                                                              <  1:61493/149‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTgg                                                                                                              <  3:22031/149‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTCGCCCAAACCCTCGAACTgg                                                                                                              <  1:261487/149‑1 (MQ=255)
                               ggCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGc                                                                                                             >  7:87707/1‑149 (MQ=255)
                                gCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGcc                                                                                                            >  1:16003/1‑149 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAccc                                                                                                                       >  4:57875/1‑136 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAccc                                                                                                                       <  3:57875/136‑1 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCg                                                                                                          >  2:241243/1‑149 (MQ=255)
                                         cacaTATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGc                                                                                                    <  6:181397/148‑1 (MQ=255)
                                          acaTATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAg                                                                                                  >  1:340370/1‑149 (MQ=255)
                                            atatTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                         >  6:263996/1‑124 (MQ=255)
                                            atatTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCGTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                         <  5:263996/124‑1 (MQ=255)
                                                ttCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGc                                                                                            >  2:62484/1‑149 (MQ=255)
                                                      ttATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCAAAAAACCCTCGAAATGGCCCGGTATGCAGGTTCGCTGAt                                                                                        >  1:124937/1‑148 (MQ=255)
                                                      ttATGACGACGACTAAAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCACTTGTTCGATGATCGAGTGCTCTACAACGGTCAACAGTTGGCTCTGGTGGTTGCCCAAACACTCGCACTGGCCCGGTATGCAGGTTCGCTGAt                                                                                        >  8:197924/1‑147 (MQ=255)
                                                         tGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGa                                                                                   <  3:339615/149‑1 (MQ=255)
                                                         tGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGa                                                                                   <  6:130119/149‑1 (MQ=255)
                                                           acgacgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATc                                                                                 <  6:86505/149‑1 (MQ=255)
                                                           acgacgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGCATGCAGGTTCGCTGATCAGGATc                                                                                 >  2:223979/1‑149 (MQ=255)
                                                               cgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAg                                                                              >  3:209499/1‑148 (MQ=255)
                                                                 acTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTAt                                                                           <  2:128082/149‑1 (MQ=255)
                                                                    aCAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                        <  8:158314/149‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGCTCGAGCGCTCTACAACGGTCAACCGTTGGCGCTGGTGGGTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                        <  6:169711/148‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                        <  2:214514/148‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGa                                                                               <  3:319235/141‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGa                                                                               >  4:319235/1‑141 (MQ=255)
                                                                        cgacgccgacgcCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCgag                                                                     >  2:279736/1‑148 (MQ=255)
                                                                        cgacgccgacgcCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGAGCGAGTATGCCgag                                                                     >  6:365193/1‑148 (MQ=255)
                                                                           cgccgacgccgCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGaggagg                                                                 <  8:87707/149‑1 (MQ=255)
                                                                                  gccgcGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTCGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCccac                                                            >  5:32491/1‑147 (MQ=255)
                                                                                     gcgcAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGGGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAgac                                                       <  2:274473/149‑1 (MQ=255)
                                                                                              gcgcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCAGAGGAGGCCCACCAGACCGACCTTgg                                              >  6:244749/1‑149 (MQ=255)
                                                                                               cgcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGATGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGCCCTTGgc                                             >  7:245976/1‑149 (MQ=255)
                                                                                                gcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGgcg                                            <  4:209499/149‑1 (MQ=255)
                                                                                                       aGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                      >  6:175726/1‑148 (MQ=255)
                                                                                                        ggCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCaa                                    <  1:62484/149‑1 (MQ=255)
                                                                                                          cccTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCaa                                    >  3:29065/1‑147 (MQ=255)
                                                                                                                ttCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCa                            >  4:208523/1‑149 (MQ=255)
                                                                                                                  cGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCATACCCTCGAACTGTCCCGGTATGCAGTTTCGCTGATCAGGTTCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGGGCATGCCAACAGTCGCAt                           <  7:197924/148‑1 (MQ=255)
                                                                                                                   gatgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATg                          >  4:300377/1‑148 (MQ=255)
                                                                                                                    atgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGCGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgca                        <  5:365193/149‑1 (MQ=255)
                                                                                                                    atgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgca                        <  3:208523/149‑1 (MQ=255)
                                                                                                                      gatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgcag                       >  3:263530/1‑148 (MQ=255)
                                                                                                                        tCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCAc                    <  1:241243/149‑1 (MQ=255)
                                                                                                                             tGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                      >  7:203586/1‑126 (MQ=255)
                                                                                                                             tGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                      <  8:203586/126‑1 (MQ=255)
                                                                                                                               agcgccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCAt                                         >  7:40288/6‑121 (MQ=255)
                                                                                                                               agcgccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCAt                                         <  8:40288/116‑1 (MQ=255)
                                                                                                                                tctACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCa               >  6:287472/1‑146 (MQ=255)
                                                                                                                                  tACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAgaga           >  8:219230/1‑148 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGc                                       >  3:310698/3‑120 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGc                                       <  4:310698/118‑1 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAgaga           >  3:132514/3‑148 (MQ=255)
                                                                                                                                   ccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGCGGAGGCCCACCAGACCGACCTTGGCGAATGCCAACAGTCGCATGCAGCACCCGCAGAGaca         >  7:72284/2‑149 (MQ=255)
                                                                                                                                    cAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGacac        <  4:263530/149‑1 (MQ=255)
                                                                                                                                        ggTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACAccgcc    <  2:124937/149‑1 (MQ=255)
                                                                                                                                          tCAACCGTTGGCGATGGTGGTTCCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGACGACGAGGCACACCAGACCCACCTTGGCGCATGCCCACAGTCGCATTCAGCACCCGCAGAGACAccgccgc  >  8:301884/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACACCGCCGC  >  NC_002947/3744897‑3745183

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