Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,028 G→A 9.4% intergenic (+56/‑464) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,0280GA 9.4% 192.2 / 7.3 72intergenic (+56/‑464)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base A (4/3);  total (36/36)
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.97e-01

ACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCAT  >  NC_002947/1898879‑1899174
                                                                                                                                                       |                                                                                                                                                  
accacTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAg                                                                                                                                                    >  1:1673371/1‑150 (MQ=32)
       gATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACCCCCGGAAGTAAAAAAGAAAGCGACGTTAGGGTCGCt                                                                                                                                             >  2:664942/1‑150 (MQ=25)
       gATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTCGGGTCGCt                                                                                                                                             >  2:664793/1‑150 (MQ=21)
            ataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCtttttt                                                                                                                                        >  1:1853690/1‑150 (MQ=34)
               gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAt                                                                                                                                     <  1:24022/150‑1 (MQ=37)
                cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAtt                                                                                                                                    >  1:1822252/1‑150 (MQ=38)
                 gcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                   >  1:1711455/1‑150 (MQ=39)
                   aGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                                                                                 <  1:1597949/150‑1 (MQ=40)
                      ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGca                                                                                                                              <  2:15272/150‑1 (MQ=255)
                       gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                   >  2:165507/1‑144 (MQ=39)
                       gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                   <  1:165507/144‑1 (MQ=39)
                          tttcAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCTTTCTTTtttacttccgggt                                                                                                                          <  2:309510‑M1/149‑14 (MQ=255)
                          gtacAGTCGGTTAGAAGACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGCACGCCGGCCAAAAGCTCACTAAACGAGCCCACTGGACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatcat                                                                                                                          <  1:1572108/147‑1 (MQ=14)
                          gTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCCCTTCCAGGGCCCTCTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCTTTCTTTt                                                                                                                                      <  2:1054381‑M1/138‑2 (MQ=255)
                          gTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCTTTCTTTt                                                                                                                                      >  1:1054381‑M1/1‑137 (MQ=255)
                           ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATg                                                                                                                         >  2:88823/1‑150 (MQ=255)
                                tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATc                                                                                                                    <  2:1853690/150‑1 (MQ=255)
                                   gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGtt                                                                                                                 >  1:684279/1‑150 (MQ=255)
                                          tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGccc                                                                                                          >  1:1776647/1‑150 (MQ=255)
                                          tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCTTTCTTTtttacttccgg                                                                                                                            >  1:1547116‑M1/1‑121 (MQ=255)
                                          tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCTTTCTTTtttacttccgg                                                                                                                            <  2:1547116‑M1/132‑12 (MQ=255)
                                           aCCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                         <  2:1883664/150‑1 (MQ=255)
                                             cGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGc                                                                                                       >  1:1608171/1‑150 (MQ=255)
                                              ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                                <  2:756707/108‑1 (MQ=38)
                                              ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                                >  1:756707/1‑108 (MQ=38)
                                               gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTa                                                                                                                                      >  1:1150967/1‑117 (MQ=39)
                                               gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTa                                                                                                                                      <  2:1150967/117‑1 (MQ=39)
                                                     cACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                               <  2:1776647/150‑1 (MQ=255)
                                                      aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                              >  1:1558865/1‑150 (MQ=255)
                                                          cGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAg                                                                                          >  1:265775/1‑150 (MQ=255)
                                                                tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgcgc                                                                                    <  1:1776783/150‑1 (MQ=255)
                                                                  gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCa                                                                                  <  2:1673371/150‑1 (MQ=255)
                                                                     ggTTCGAGTCCCGTCCGCTGCGCCATAATCGCAACCTGCGTGGCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCActc                                                                               <  2:151884/150‑1 (MQ=25)
                                                                       ttCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtctt                                                                             <  2:732746/150‑1 (MQ=255)
                                                                       ttCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtctt                                                                             >  1:1163866/1‑150 (MQ=255)
                                                                                cGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCa                                                                    >  1:1019635/1‑150 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCAc                                                                                 <  1:1862083/132‑1 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCAc                                                                                 >  2:1862083/1‑132 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCAc                                                                                 >  2:1862228/1‑132 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCAc                                                                                 <  1:1862228/132‑1 (MQ=255)
                                                                                      cTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCa                                                              <  2:1822252/150‑1 (MQ=255)
                                                                                         cgcCATATTCCACTTCCAGGGCCCACTGAACATGCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCa                                                           <  2:109899/150‑1 (MQ=255)
                                                                                         cgcCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCa                                                           <  1:88823/150‑1 (MQ=255)
                                                                                             atatTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCaa                                                                                                                      >  2:1552407/1‑87 (MQ=255)
                                                                                             atatTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCaa                                                                                                                      <  1:1552407/87‑1 (MQ=255)
                                                                                                ttCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                                    >  1:1437520/1‑150 (MQ=255)
                                                                                                 tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                   <  2:1019635/150‑1 (MQ=255)
                                                                                                  ccACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAcc                                                  >  1:1474144/1‑150 (MQ=255)
                                                                                                        tCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTcc                                             <  2:1474144/150‑1 (MQ=255)
                                                                                                                  cACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCAcc                                    >  2:1862645/1‑150 (MQ=255)
                                                                                                                   aCTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCg                                   <  2:684279/150‑1 (MQ=255)
                                                                                                                     tGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAGACCCGTCCTCGATCACCGGc                                 >  1:1737071/1‑150 (MQ=255)
                                                                                                                       aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCGACTGTTCTTCCCGCAAGATCAGCGCCTCAGACCCGTCCTCGATGACCGGCAt                               >  2:1039779/1‑150 (MQ=255)
                                                                                                                            ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                          >  2:94688/1‑150 (MQ=255)
                                                                                                                            ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                          >  2:545101/1‑150 (MQ=255)
                                                                                                                            ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                          >  2:1219657/1‑150 (MQ=255)
                                                                                                                              ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTgg                        <  1:1398738/150‑1 (MQ=255)
                                                                                                                                   tAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTgttacttccgggtgttcagtgggccctggaagtggaatatggcgcagcggacgggactcgaacccgcgacccccggcgtgacaggccggtattctaaccgactgaactaccgctgc                   >  1:458669‑M1/1‑34 (MQ=255)
                                                                                                                                     aaaaaGAAAGCGACCTGAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATgg                 <  2:173093/150‑1 (MQ=255)
                                                                                                                                       aaagaaagCGACCTTAGGGTCGCTTTTTTTtttacttccgggtgttcagtgggccctggaagtggaatatggcgcagcggacgggactcgaacccgcgacccccggcgtgacaggccggtattctaaccgactgaactacc                        <  2:1305507‑M1/141‑112 (MQ=255)
                                                                                                                                       aaagaaagCGACCTTAGGGTCGCTTTTTTTtttacttccgggtgttcagtgggccctggaagtggaatatggcgcagcggacgggactcgaacccgcgacccccggcgtgacaggccggtattctaaccgactgaactacc                        >  1:1305507‑M1/1‑30 (MQ=255)
                                                                                                                                       aaagaaagCGACCTTAGGGTCGCTTTTTTTATTTTTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGa               <  2:600533/150‑1 (MQ=255)
                                                                                                                                       aaagaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGa               <  2:1558865/150‑1 (MQ=255)
                                                                                                                                        aagaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAGACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAg              <  2:1737071/150‑1 (MQ=255)
                                                                                                                                        aagaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAg              >  2:1684230/1‑150 (MQ=255)
                                                                                                                                         agaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAgg             <  2:1569343/150‑1 (MQ=255)
                                                                                                                                         agaaagCGACCTTAGAGTCGCTTTTTTTgttacttccgggtgttcagtgggccctggaagtggaatatggcgcagc                                                                                       <  1:900291‑M1/76‑49 (MQ=255)
                                                                                                                                         agaaagCGACCTTAGAGTCGCTTTTTTTgttacttccgggtgttcagtgggccctggaagtggaatatggcgcagc                                                                                       >  2:900291‑M1/1‑28 (MQ=255)
                                                                                                                                          aaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg            >  2:1842437/2‑150 (MQ=255)
                                                                                                                                           aaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGTCACTCTGCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGc           >  1:231971/1‑150 (MQ=255)
                                                                                                                                            aagCGACCCTAAGGTCGCTTTCTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCa          >  1:1243129/1‑150 (MQ=255)
                                                                                                                                               cGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTg       <  1:1684230/150‑1 (MQ=255)
                                                                                                                                                    ttAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAATCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCAt  <  2:265775/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
ACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCAT  >  NC_002947/1898879‑1899174

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: