Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,053,520 A→T 7.1% Q17L (CAG→CTG)  PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,053,5200AT 7.1% 206.7 / 5.9 70Q17L (CAG→CTG) PP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base T (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGGG  >  NC_002947/4053371‑4053668
                                                                                                                                                     |                                                                                                                                                    
aTACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAACCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCa                                                                                                                                                      <  2:71916/150‑1 (MQ=255)
aTACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCa                                                                                                                                                      >  1:2365633/1‑150 (MQ=255)
 tACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAg                                                                                                                                                     >  2:162331/1‑150 (MQ=255)
    tATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGcc                                                                                                                                                  >  1:1782790/1‑150 (MQ=255)
          cGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTg                                                                                                                                            <  2:1828491/150‑1 (MQ=255)
           gCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCgggg                                                                                                                                           >  1:942613/1‑150 (MQ=255)
            ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGc                                                                                                                                          <  2:1396262/150‑1 (MQ=255)
             cATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCgcc                                                                                                                                            <  2:1949412/147‑3 (MQ=255)
             cATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCgcc                                                                                                                                            >  1:1949412/1‑145 (MQ=255)
               tttCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCAt                                                                                                                                       <  2:1199775/150‑1 (MQ=255)
               tttCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCAt                                                                                                                                       >  2:2246472/1‑150 (MQ=255)
                 tCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGa                                                                                                                                     <  1:398737/150‑1 (MQ=255)
                   cACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGc                                                                                                                                   >  2:1442823/1‑150 (MQ=255)
                      ggCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGaa                                                                                                                                <  2:2195675/150‑1 (MQ=255)
                         aTCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGcc                                                                                                                                         >  2:308097/1‑138 (MQ=255)
                         aTCCCTAGTCTGGCTCCAGGCAACGGCGGTCCCCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGcc                                                                                                                                         <  1:308097/138‑1 (MQ=255)
                              tAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAACAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGCCCCCGGCCAGGCCGGCGTGGCCATCACCGAACACCCGAg                                                                                                                        >  2:877868/1‑150 (MQ=255)
                                        ccAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTGGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGc                                                                                                              <  2:83255/150‑1 (MQ=255)
                                           ggCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTg                                                                                                           >  2:373553/1‑150 (MQ=255)
                                            gCAAGGGCGGGCCTACAAGCGGACCGGAGAAAAACCGCCTGAAGAGGGGCGCGGGCGTGGCATGGGGGAGGTTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCACGAGCGAACACCTGAGCAAGCTTGGCGTGc                                                                                                          <  1:877868/150‑1 (MQ=255)
                                            gCAACGGCGGTCCTCCAAGCCGACCTGAAAAAAGCTGTCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGc                                                                                                          <  2:188655/150‑1 (MQ=255)
                                            gCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTGGTAGGCTCTGGCTAAGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGc                                                                                                          >  2:663620/1‑150 (MQ=255)
                                                  gcggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAcc                                                                                                    <  2:247446/150‑1 (MQ=255)
                                                       cctccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTa                                                                                               <  1:1008573/150‑1 (MQ=255)
                                                         tccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACt                                                                                             >  2:900442/1‑150 (MQ=255)
                                                            agaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGa                                                                                          <  2:599947/148‑1 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGa                                                                                          <  2:1782790/150‑1 (MQ=255)
                                                            aaaCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGa                                                                                          <  2:2106627/150‑1 (MQ=255)
                                                              aCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAgc                                                                                        <  1:2246472/150‑1 (MQ=255)
                                                                   ccTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAc                                                                                   <  1:1320301/150‑1 (MQ=255)
                                                                    cTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAcc                                                                                  >  2:1632860/1‑150 (MQ=255)
                                                                        cAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACCCCTGAGCAAGATTGGCGTGCCGTACTTGGTACTGGAGCGCAACCGCAt                                                                              <  2:2365633/150‑1 (MQ=255)
                                                                         aaaaaCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATc                                                                             >  1:1039399/1‑150 (MQ=255)
                                                                                gCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAgg                                                                      >  1:1271184/1‑150 (MQ=255)
                                                                                gCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAgg                                                                      <  1:1442823/150‑1 (MQ=255)
                                                                                  cTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAAGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGcc                                                                    >  2:1136067/1‑150 (MQ=255)
                                                                                     aaCAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTgg                                                                 >  1:2359686/1‑150 (MQ=255)
                                                                                      aCAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGc                                                                >  2:289246/1‑150 (MQ=255)
                                                                                            gCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTCCCg                                                          >  2:139147/1‑150 (MQ=255)
                                                                                             cTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCgg                                                         >  2:1085977/1‑150 (MQ=255)
                                                                                                   cATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAATACTGCCGCTCAGGCCTGGCGTACCGGCCGCTg                                                   <  1:663620/150‑1 (MQ=255)
                                                                                                   cATGACACGGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTg                                                   <  1:511832/150‑1 (MQ=255)
                                                                                                      gACACTGCACAACCTCGAAATCGTCACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGa                                                <  1:1071525/150‑1 (MQ=255)
                                                                                                         acTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCTCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCt                                                                                                                           <  2:605476/72‑1 (MQ=255)
                                                                                                         acTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCTCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCt                                                                                                                           >  1:605476/1‑72 (MQ=255)
                                                                                                         acTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                             >  2:2449673/1‑150 (MQ=255)
                                                                                                         acTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                             >  1:949325/1‑150 (MQ=255)
                                                                                                         acTGAACAACCTCGAAATCGACACCATCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTc                                             >  1:1105068/1‑150 (MQ=255)
                                                                                                             aacaacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGc                                                                     >  1:579590/1‑122 (MQ=255)
                                                                                                             aacaacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGc                                                                     <  2:579590/122‑1 (MQ=255)
                                                                                                                aacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTc                                      <  2:2359686/150‑1 (MQ=255)
                                                                                                                   cTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTa                                                             <  2:1306180/124‑1 (MQ=255)
                                                                                                                   cTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTa                                                             >  1:1306180/1‑124 (MQ=255)
                                                                                                                    tCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTg                                         >  2:873029/1‑143 (MQ=255)
                                                                                                                    tCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTg                                         <  1:873029/143‑1 (MQ=255)
                                                                                                                     cGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCaa                                 <  1:1085977/150‑1 (MQ=255)
                                                                                                                     cGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCaa                                 <  1:2449673/150‑1 (MQ=255)
                                                                                                                          tCGACACTCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCGTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcc                            >  2:1486090/1‑150 (MQ=255)
                                                                                                                          tCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcc                            >  1:1340703/1‑150 (MQ=255)
                                                                                                                           cGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggccc                           <  1:1451927/150‑1 (MQ=255)
                                                                                                                            gACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcccg                          >  2:180568/1‑150 (MQ=255)
                                                                                                                             acacCCTCGTCGTCGGCGCCGGCCAGGCCGGTGTGGCCATGAGCGAACACCTGAGCACGCTTGGCGTGCCGAACCTGGTACTGGAGCGCAACCGCATCGCGCAGTCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcccgg                         <  2:942613/150‑1 (MQ=255)
                                                                                                                                 ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGa                                                                                          >  1:875579/1‑81 (MQ=255)
                                                                                                                                 ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGa                                                                                          <  2:875579/81‑1 (MQ=255)
                                                                                                                                 ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTg                     >  2:2096291/1‑150 (MQ=255)
                                                                                                                                  cTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTc                                      >  1:2223579/1‑132 (MQ=255)
                                                                                                                                  cTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTc                                      <  2:2223579/132‑1 (MQ=255)
                                                                                                                                        gtcgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGAc              >  1:1096843/1‑150 (MQ=255)
                                                                                                                                           cgcgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGc           >  1:1063949/2‑150 (MQ=255)
                                                                                                                                           cacgcTGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcggcac                           <  2:2170813/132‑3 (MQ=255)
                                                                                                                                           cacgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGc           >  1:1833869/3‑150 (MQ=255)
                                                                                                                                                    cAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCgggg  >  1:1964927/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGGG  >  NC_002947/4053371‑4053668

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