Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,512 A→C 7.1% intergenic (+23/‑28) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5120AC 7.1% 754.9 / 39.5 224intergenic (+23/‑28)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (111/95);  new base C (10/6);  total (122/102)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 226 total aligned reads displayed.

AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTA  >  NC_000913/3020364‑3020660
                                                                                                                                                     |                                                                                                                                                    
aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCa                                                                                                                                                      <  2:957460/149‑1 (MQ=255)
 gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAg                                                                                                                                                     >  2:12636/1‑149 (MQ=255)
 gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAg                                                                                                                                                     >  1:1300578/1‑149 (MQ=255)
  cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAgg                                                                                                                                                    <  3:1546618/149‑1 (MQ=255)
    tGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGa                                                                                                                                                  >  1:443764/1‑149 (MQ=255)
       tGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAaga                                                                                                                                               >  1:1257981/1‑149 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAGGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCGC‑GCAGGGAagag                                                                                                                                              <  4:1458025/149‑1 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                              <  3:1341675/149‑1 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                              <  2:111034/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                             <  3:1446639/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                             <  1:80185/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAgg                                                                                                                                             <  3:703109/149‑1 (MQ=255)
          tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAggg                                                                                                                                            <  3:1720300/149‑1 (MQ=255)
           cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGt                                                                                                                                           <  3:1412620/149‑1 (MQ=255)
            aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                          >  1:861015/1‑149 (MQ=255)
                     aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                                 >  4:1721827/1‑149 (MQ=255)
                     aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGCCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                                 >  2:1013043/1‑149 (MQ=255)
                          cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                            <  4:499427/149‑1 (MQ=255)
                             tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                         <  4:1126804/149‑1 (MQ=255)
                              gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAggg                                                                                                                                            >  4:198862/1‑129 (MQ=255)
                              gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatg                                                                                                                        >  4:1418157/1‑149 (MQ=255)
                               aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                       >  3:1334207/1‑149 (MQ=255)
                               aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                       >  1:533636/1‑149 (MQ=255)
                                aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                      >  1:1038928/1‑149 (MQ=255)
                                aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                      <  2:1580671/149‑1 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     >  2:515459/1‑149 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  3:1710048/149‑1 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  4:938999/149‑1 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  2:1522941/149‑1 (MQ=255)
                                  aCCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                    >  3:1104871/1‑149 (MQ=255)
                                    cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                  <  3:1677849/149‑1 (MQ=255)
                                    cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                  <  1:1624457/149‑1 (MQ=255)
                                          tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                            >  1:455935/1‑149 (MQ=255)
                                            tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                          <  1:676072/149‑1 (MQ=255)
                                              aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAAt                                                                                                                                      >  1:707399/1‑119 (MQ=255)
                                              aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAAt                                                                                                                                      <  2:707399/119‑1 (MQ=255)
                                              aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                        <  4:370092/149‑1 (MQ=255)
                                               cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                                       >  4:1228743/1‑149 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                         >  1:871455/1‑130 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                         <  2:871455/130‑1 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      >  1:195208/1‑149 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      >  3:1593820/1‑149 (MQ=255)
                                                 tCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAg                                                                                                                                   >  2:436608/1‑119 (MQ=255)
                                                    aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCg                                                                                                  >  2:161478/1‑149 (MQ=255)
                                                        cGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAAt                                                                                                                                      >  3:829492/1‑109 (MQ=255)
                                                           ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                           <  4:279008/117‑1 (MQ=255)
                                                           ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                           >  3:279008/1‑117 (MQ=255)
                                                               aTTTATGGACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGCTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGt                                                                                           >  4:1085833/1‑145 (MQ=255)
                                                               aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                       >  1:282683/1‑149 (MQ=255)
                                                                tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATTGATGCGCTGGCATAAATGACAGATGCACTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                      >  3:206633/1‑149 (MQ=255)
                                                                tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                      <  3:193054/149‑1 (MQ=255)
                                                                 ttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                                     <  2:307401/149‑1 (MQ=255)
                                                                       acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                        >  2:663721/1‑124 (MQ=255)
                                                                       acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                        <  1:663721/124‑1 (MQ=255)
                                                                         aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                             >  4:443441/1‑149 (MQ=255)
                                                                         aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGAAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                             <  2:383897/149‑1 (MQ=255)
                                                                             cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                         >  4:957343/1‑149 (MQ=255)
                                                                               gAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     >  3:807562/1‑103 (MQ=255)
                                                                               gAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  4:807562/103‑1 (MQ=255)
                                                                                aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                      <  3:76250/149‑1 (MQ=255)
                                                                                aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                      <  4:1674598/149‑1 (MQ=255)
                                                                                 aaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAg                                                                     <  3:443441/149‑1 (MQ=255)
                                                                                     ccGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                                 >  4:660021/1‑149 (MQ=255)
                                                                                              gTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACaa                                                                                                               <  4:581095/94‑1 (MQ=255)
                                                                                              gTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACaa                                                                                                               >  3:581095/1‑94 (MQ=255)
                                                                                               tATGTGGCGTCGGATGGATGCGCTAGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGCAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGa                                                       >  3:1082449/1‑149 (MQ=255)
                                                                                                  gtgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGa                                                    >  3:1133252/1‑149 (MQ=255)
                                                                                                   tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                                   <  1:1086450/149‑1 (MQ=255)
                                                                                                       gTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                         <  1:100340/91‑1 (MQ=255)
                                                                                                        tCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGt                                              <  3:1228743/149‑1 (MQ=255)
                                                                                                           gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTg                                           <  1:575856/149‑1 (MQ=255)
                                                                                                                atgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                                      >  3:894340/1‑149 (MQ=255)
                                                                                                                    gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGAct                                                                                               >  1:1386690/1‑88 (MQ=255)
                                                                                                                    gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGAct                                                                                               <  2:1386690/88‑1 (MQ=255)
                                                                                                                      tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                                <  2:1201597/149‑1 (MQ=255)
                                                                                                                      tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                                >  3:851882/1‑149 (MQ=255)
                                                                                                                        gCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgc                              >  3:491254/1‑149 (MQ=255)
                                                                                                                          aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                            >  3:402328/1‑149 (MQ=255)
                                                                                                                              aTGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGATGCCGTCTGGTTCGCCGGGGGCAGCaaa                        >  3:1188912/1‑149 (MQ=255)
                                                                                                                                   aGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                   >  4:82676/1‑149 (MQ=255)
                                                                                                                                    ttaacccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt                  >  3:1278585/5‑149 (MQ=255)
                                                                                                                                    gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt                  <  4:697412/149‑1 (MQ=255)
                                                                                                                                    gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt                  <  1:165465/149‑1 (MQ=255)
                                                                                                                                     aTGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATg                 >  1:358686/1‑149 (MQ=255)
                                                                                                                                        cccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAAc                                                                          >  3:910549/1‑89 (MQ=39)
                                                                                                                                        cccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                   <  1:1449045/144‑1 (MQ=255)
                                                                                                                                        cccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                   >  2:1449045/1‑144 (MQ=255)
                                                                                                                                         ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  1:988619/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  2:234781/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  2:1671720/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  3:339627/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             <  4:878251/149‑1 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  4:103569/1‑149 (MQ=255)
                                                                                                                                          ctctTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca            <  4:139189/149‑1 (MQ=255)
                                                                                                                                          ctctTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca            >  3:139189/1‑149 (MQ=255)
                                                                                                                                          ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACACTTTTTCAGTCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCATCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca            <  2:766158/149‑1 (MQ=255)
                                                                                                                                           tctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCAGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac           <  4:1188912/149‑1 (MQ=255)
                                                                                                                                                cc‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc      >  4:1555605/1‑149 (MQ=255)
                                                                                                                                                cc‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc      >  1:41910/1‑149 (MQ=255)
                                                                                                                                                     aGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTa  <  1:1177198/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTA  >  NC_000913/3020364‑3020660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑