Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,962,544 T→A 6.5% L233Q (CTG→CAG)  PP_3492 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,962,5440TA 6.5% 248.0 / 7.7 77L233Q (CTG→CAG) PP_3492acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base A (2/3);  total (38/39)
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.15e-01

GCGTTCCTGGTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATT  >  NC_002947/3962396‑3962692
                                                                                                                                                    |                                                                                                                                                    
gCGTTCCTGGTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTg                                                                                                                                                     <  1:1841532/150‑1 (MQ=255)
        ggTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGc                                                                                                                                             <  2:449538/150‑1 (MQ=255)
        ggTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCACCTGGCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGc                                                                                                                                             >  1:610772/1‑150 (MQ=255)
         gTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGcc                                                                                                                                            >  1:502708/1‑150 (MQ=255)
             ccACCGACACCCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGGGGCAATGGCCAGGCCCTTGccacgc                                                                                                                                        <  1:274268/150‑5 (MQ=255)
             ccACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGt                                                                                                                                        >  1:333187/1‑150 (MQ=255)
              cACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                       >  1:1704283/1‑150 (MQ=255)
                        ccTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGaa                                                                                                                             <  2:1704283/150‑1 (MQ=255)
                         cTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAg                                                                                                                            <  2:562147/150‑1 (MQ=255)
                             gTTCAATGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAggcgg                                                                                                                        >  1:1285356/1‑150 (MQ=255)
                               tCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGTCCATTGCCCTGTCCAACCTCGAAGGCGGcc                                                                                                                      <  2:800361/150‑1 (MQ=255)
                               tCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGcc                                                                                                                      >  1:69695/1‑150 (MQ=255)
                               tCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGcc                                                                                                                      <  2:865728/150‑1 (MQ=255)
                                         ccGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCAt                                                                                                            >  1:115423/1‑150 (MQ=255)
                                          cGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGcc                                                                                                                                            <  2:1822935/117‑1 (MQ=255)
                                          cGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGcc                                                                                                                                            >  1:1822935/1‑117 (MQ=255)
                                            cAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCCCACCTCCCCCGGGCAATGGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTgc                                                                                                         <  2:50195/150‑1 (MQ=255)
                                                   cacaAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACGTCGAAGGCGGCCGTATCGGCATTGCCGCGCAg                                                                                                  >  1:475125/1‑150 (MQ=255)
                                                   cacaAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAg                                                                                                  <  1:1796444/150‑1 (MQ=255)
                                                    acaAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAgg                                                                                                                           <  1:1703867/124‑1 (MQ=255)
                                                    acaAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAgg                                                                                                                           >  2:1703867/1‑124 (MQ=255)
                                                      aaGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGcc                                                                                               <  1:1526759/150‑1 (MQ=255)
                                                      aaGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGcc                                                                                               >  2:1940156/1‑150 (MQ=255)
                                                          ggggCATTCGCGGCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTGCAGCCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTgg                                                                                           <  2:1728989/149‑1 (MQ=255)
                                                             gCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCa                                                                                        >  2:1199305/1‑150 (MQ=255)
                                                                ttCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAATGGCCAGGccctt                                                                                                                                               >  2:1666035/1‑89 (MQ=255)
                                                                ttCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAATGGCCAGGCCCtt                                                                                                                                               <  1:1666035/92‑1 (MQ=255)
                                                                ttCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCg                                                                                     >  2:529254/1‑150 (MQ=255)
                                                                        tCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGcgcc                                                                             >  1:1418538/1‑150 (MQ=255)
                                                                          cGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAgg                                                                                                                           >  1:1531132/1‑102 (MQ=255)
                                                                          cGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAgg                                                                                                                           <  2:1531132/102‑1 (MQ=255)
                                                                            aTACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTgg                                                                                           <  2:1582881/132‑1 (MQ=255)
                                                                            aTACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTgg                                                                                           >  1:1582881/1‑132 (MQ=255)
                                                                            aTACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCAt                                                                         <  1:144462/150‑1 (MQ=255)
                                                                                 tGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGaa                                                                    >  2:1524988/1‑150 (MQ=255)
                                                                                    gcgGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAgcg                                                                 >  2:765503/1‑150 (MQ=255)
                                                                                      ggTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAATGGCCAGGCCCTTGc                                                                                                                                             <  1:871858/72‑1 (MQ=38)
                                                                                      ggTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAATGGCCAGGCCCTTGc                                                                                                                                             >  2:871858/1‑72 (MQ=38)
                                                                                        tCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGgcgc                                                             <  2:502708/150‑1 (MQ=255)
                                                                                          aCCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTg                                                           >  1:117165/1‑150 (MQ=255)
                                                                                            cTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTgg                                                          <  2:365323/149‑1 (MQ=255)
                                                                                            cTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTgg                                                          >  1:365323/1‑149 (MQ=255)
                                                                                            cTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGt                                                         >  2:1659963/1‑150 (MQ=255)
                                                                                            cTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGt                                                         <  2:1596634/150‑1 (MQ=255)
                                                                                            cTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGt                                                         <  2:1863296/150‑1 (MQ=255)
                                                                                             ttCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGCTc                                                        <  1:1940156/150‑1 (MQ=255)
                                                                                                  cAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTc                                                   <  1:854473/150‑1 (MQ=255)
                                                                                                     cTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCa                                                                                                   <  1:486290/99‑1 (MQ=255)
                                                                                                     cTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCa                                                                                                   >  2:486290/1‑99 (MQ=255)
                                                                                                     cTGCCGCATTCCTGCCGCCAATATACTCGGTGAGCGTGGCAAGGGACTGGCCATTGACCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCAGTCGAAGCGGCGCTGGTCTACTcgcg                                                >  1:1858407/1‑150 (MQ=255)
                                                                                                        ccGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGc                                                                                    >  2:1541145/1‑111 (MQ=255)
                                                                                                        ccGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGc                                                                                    <  1:1541145/111‑1 (MQ=255)
                                                                                                          gCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGAcc                                           >  1:342196/1‑150 (MQ=255)
                                                                                                          gCATTCCTGACGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAg                                                                                                  <  2:1392995/95‑1 (MQ=255)
                                                                                                          gCATTCCTGACGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAg                                                                                                  >  1:1392995/1‑95 (MQ=255)
                                                                                                                cTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTg                                                                                            >  2:1802269/1‑95 (MQ=255)
                                                                                                                cTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTg                                                                                            <  1:1802269/95‑1 (MQ=255)
                                                                                                                cTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAgc                                                                  >  2:128548/1‑121 (MQ=255)
                                                                                                                cTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAgc                                                                  <  1:128548/121‑1 (MQ=255)
                                                                                                                   ccgccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGt                                  >  1:1223027/1‑150 (MQ=255)
                                                                                                                      ccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTg                               <  2:610772/150‑1 (MQ=255)
                                                                                                                       cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAg                                   >  1:1215821/1‑145 (MQ=255)
                                                                                                                       cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAg                                   <  2:1215821/145‑1 (MQ=255)
                                                                                                                          tatCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCaa                           <  2:1230842/150‑1 (MQ=255)
                                                                                                                          tatCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCaa                           >  1:1230842/1‑150 (MQ=255)
                                                                                                                            tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGc                         <  2:1223027/150‑1 (MQ=255)
                                                                                                                             ccTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGcc                        <  2:1418538/150‑1 (MQ=255)
                                                                                                                              cTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCg                       >  1:1915833/1‑150 (MQ=255)
                                                                                                                                 ggTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATc                    >  2:1304907/1‑150 (MQ=255)
                                                                                                                                   tGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCaa                  >  1:388649/1‑150 (MQ=255)
                                                                                                                                    gAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAAc                 <  1:1659963/150‑1 (MQ=255)
                                                                                                                                          ggCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAacac           >  1:1130058/1‑150 (MQ=255)
                                                                                                                                               gggCCTGGCCATCGCCCTGTCCAACCTCCAAGGCCGCCATATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGATGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCagag      <  2:475125/150‑1 (MQ=255)
                                                                                                                                                 gCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCAtt                                                                        <  2:327998/82‑1 (MQ=255)
                                                                                                                                                 gCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCAtt                                                                        >  1:327998/1‑82 (MQ=255)
                                                                                                                                                 gCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCa    <  1:1199305/150‑1 (MQ=255)
                                                                                                                                                   cTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCAtt  <  2:388649/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGTTCCTGGTGCCCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAAGCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATT  >  NC_002947/3962396‑3962692

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