Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,321 G→C 11.5% L458V (CTG→GTG)  PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,3210GC11.5% 206.3 / 17.2 70L458V (CTG→GTG) PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base C (4/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

TTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTG  >  NC_002947/3189176‑3189466
                                                                                                                                                 |                                                                                                                                                 
ttGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGc                                                                                                                                               >  1:838568/1‑150 (MQ=255)
   ggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                            <  2:3200/150‑1 (MQ=255)
     ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCTGTGa                                                                                                                                          >  1:14084/1‑150 (MQ=255)
       atgaggGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCCCGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgtgggg                                                                                                                                        >  2:191688/1‑147 (MQ=255)
         gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                      >  2:470673/1‑150 (MQ=255)
             gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAgcaaagg                                                                                                                                   <  2:230955/149‑6 (MQ=255)
             gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAgcaaagg                                                                                                                                   >  1:230955/1‑144 (MQ=255)
             gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCTGTGAgcaaaggc                                                                                                                                  <  2:14084/150‑7 (MQ=255)
              ttCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTTATTTCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGGGCCCa                                                                                                                                 >  1:78221/1‑150 (MQ=255)
                  cATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACTTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                          >  2:188303/1‑137 (MQ=255)
                  cATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACTTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                          <  1:188303/137‑1 (MQ=255)
                        aCCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCATTTGCTGTGAGCGCGCCCAATCCGTcgcg                                                                                                                       >  2:1082236/1‑150 (MQ=255)
                                   tCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGt                                                                                                            <  1:969110/150‑1 (MQ=255)
                                             gTCGGACTCGAAGCTTGCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAg                                                                                                  <  1:96600/150‑1 (MQ=255)
                                                gAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcg                                                                                               >  1:813194/1‑150 (MQ=255)
                                                    tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGa                                                                                           >  2:140563/1‑150 (MQ=255)
                                                             tcatcaCTGTGTTCACGCCCGGTCAAACCATCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGTTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCTGCCGCTCGAAGGCGCTGAAGATCGCCt                                                                                  >  1:1081348/1‑150 (MQ=255)
                                                              catcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCta                                                                                 <  2:856675/150‑2 (MQ=255)
                                                                  aCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGata                                                                             <  2:773208/150‑1 (MQ=255)
                                                                   cTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGatat                                                                            >  1:538917/1‑150 (MQ=255)
                                                                    tgtgTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATc                                                                           <  1:547922/150‑1 (MQ=255)
                                                                    tgtgTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATc                                                                           >  2:936355/1‑150 (MQ=255)
                                                                      tgtTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGt                                                                         >  2:132422/1‑150 (MQ=255)
                                                                       gtTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTc                                                                        <  2:185187/150‑1 (MQ=255)
                                                                           acgcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                    >  2:535145/1‑150 (MQ=255)
                                                                               acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtgg                                                                <  1:132422/150‑1 (MQ=255)
                                                                               acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggt                                                                  >  2:428054/1‑148 (MQ=255)
                                                                               acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGGCCTGggt                                                                  <  1:428054/148‑1 (MQ=255)
                                                                                  gTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTaa                                                             <  1:936355/150‑1 (MQ=255)
                                                                                   tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAAt                                                            >  1:110653/1‑150 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                      >  2:575964/1‑70 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                      <  1:575964/70‑1 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtg                                                                 <  2:116384/139‑1 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtg                                                                 >  1:116384/1‑139 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgg                                                      <  2:813194/150‑1 (MQ=255)
                                                                                          cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCCGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                     >  2:329949/1‑150 (MQ=255)
                                                                                          cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                     >  1:659752/1‑150 (MQ=255)
                                                                                          cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                     >  2:853376/1‑150 (MQ=255)
                                                                                           tCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCCGCCAGCGTCCGCCCGAAGGCGCTGAAGCTCGCCTCGATATCGTCCTGGGTGGGAATCCACggag                                                    >  1:1020359/1‑150 (MQ=255)
                                                                                             aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGTAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                  >  2:1130188/1‑150 (MQ=255)
                                                                                                   gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCGGGGGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAg                                            <  1:396986/150‑1 (MQ=255)
                                                                                                    ttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagc                                                                                                                  <  1:427188/79‑1 (MQ=255)
                                                                                                    ttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagc                                                                                                                  >  2:427188/1‑79 (MQ=255)
                                                                                                    ttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAgg                                           <  1:1011037/150‑1 (MQ=255)
                                                                                                         ttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAATGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                      <  2:538917/150‑1 (MQ=255)
                                                                                                         ttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                      <  2:838568/150‑1 (MQ=255)
                                                                                                            cGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCaa                                                         >  1:41629/1‑128 (MQ=255)
                                                                                                            cGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCaa                                                         <  2:41629/128‑1 (MQ=255)
                                                                                                             gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                  >  2:97677/1‑150 (MQ=255)
                                                                                                             gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                  <  1:140563/150‑1 (MQ=255)
                                                                                                             gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                  >  1:775349/1‑150 (MQ=255)
                                                                                                              tGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCCACGGAGGACAGAAGGCGATCGTGTc                                 <  1:906860/150‑1 (MQ=255)
                                                                                                                 gTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTa                              <  1:678442/150‑1 (MQ=255)
                                                                                                                     ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGACGAGGCTGAAGCGAGCCTCGAAATCTTCGTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATg                          >  1:570091/1‑150 (MQ=255)
                                                                                                                     ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                    >  1:898786/1‑108 (MQ=255)
                                                                                                                     ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                    <  2:898786/108‑1 (MQ=255)
                                                                                                                       aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCaa                                                         >  2:14031/1‑117 (MQ=255)
                                                                                                                       aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCaa                                                         <  1:14031/117‑1 (MQ=255)
                                                                                                                       aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGgc                        >  1:750661/1‑150 (MQ=255)
                                                                                                                            aTGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGaa                   <  1:470673/150‑1 (MQ=255)
                                                                                                                                 gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc              <  1:529577/150‑1 (MQ=255)
                                                                                                                                  ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGcc                                                                                   >  1:797919/1‑80 (MQ=255)
                                                                                                                                  ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGcc                                                                                   <  2:797919/80‑1 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                           <  1:737760/129‑1 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg      >  2:143217/1‑150 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGGCGTAGAt                           >  2:737760/1‑129 (MQ=255)
                                                                                                                                         tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATc                                                                                      >  2:536199/1‑70 (MQ=39)
                                                                                                                                         tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATc                                                                                      <  1:536199/70‑1 (MQ=39)
                                                                                                                                           acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                    >  1:51243/1‑86 (MQ=255)
                                                                                                                                           acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                    <  2:51243/86‑1 (MQ=255)
                                                                                                                                            cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCtt   <  1:853376/150‑1 (MQ=255)
                                                                                                                                            cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCtt   <  2:750661/150‑1 (MQ=255)
                                                                                                                                             agcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtg  >  1:1085881/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTG  >  NC_002947/3189176‑3189466

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