Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,333,941 A→G 5.7% S29S (TCA→TCG cysB → N‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,333,9410AG 5.7% 225.2 / 3.6 70S29S (TCA→TCGcysBN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (39/27);  new base G (2/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

ACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGAT  >  NC_000913/1333797‑1334088
                                                                                                                                                |                                                                                                                                                   
aCCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAAc                                                                                                                                                 >  1:428791/1‑149 (MQ=255)
aCCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAAc                                                                                                                                                 <  3:497137/149‑1 (MQ=255)
      cAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGa                                                                                                                                           <  3:270032/149‑1 (MQ=255)
       aGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGAt                                                                                                                                          <  3:155342/149‑1 (MQ=255)
         gtgtGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCa                                                                                                                                        <  1:352321/149‑1 (MQ=255)
            tgACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCGCAACCCGGGATCAGTa                                                                                                                                     >  3:565285/1‑149 (MQ=255)
            tgACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCGCAACCCGGGATCAGTa                                                                                                                                     >  3:439049/1‑149 (MQ=255)
             gACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTaa                                                                                                                                    >  4:353273/1‑149 (MQ=255)
              aCGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTaaa                                                                                                                                   <  2:337458/149‑1 (MQ=255)
                gAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACa                                                                                                                                 >  1:392396/1‑149 (MQ=255)
                           gTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGc                                                                                                                      >  2:268610/1‑149 (MQ=255)
                            ttCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCt                                                                                                                     <  2:352239/149‑1 (MQ=255)
                            ttCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCt                                                                                                                     <  3:275006/149‑1 (MQ=255)
                            ttCTGACGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCt                                                                                                                     >  2:413536/1‑149 (MQ=255)
                               tGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGa                                                                                                                  >  3:187045/1‑149 (MQ=255)
                                 aTGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAg                                                                                                                >  4:129116/1‑149 (MQ=255)
                                   ggCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGAc                                                                                                              <  3:263050/149‑1 (MQ=255)
                                         tAAGTGGATGGTTTAACATGAAATTACATCAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTa                                                                                                        >  2:351477/1‑149 (MQ=255)
                                          aaGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       <  1:52289/149‑1 (MQ=255)
                                          aaGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       <  3:172148/149‑1 (MQ=255)
                                                aTGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTc                                                                                                 >  1:146392/1‑149 (MQ=255)
                                                aTGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTc                                                                                                 >  1:562976/1‑149 (MQ=255)
                                                     ttAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAAtt                                                                                            <  4:501450/149‑1 (MQ=255)
                                                          aTGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTc                                                                                       >  3:3652/1‑149 (MQ=255)
                                                          aTGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTc                                                                                       >  1:403800/1‑149 (MQ=255)
                                                           tGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTcc                                                                                      >  1:46656/1‑149 (MQ=255)
                                                               aTTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCGCAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                                  <  4:439049/149‑1 (MQ=255)
                                                               aTTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCGCAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                                  <  4:565285/149‑1 (MQ=255)
                                                               aTTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCGGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                                  >  2:132765/1‑149 (MQ=255)
                                                               aTTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                                  <  1:548692/149‑1 (MQ=255)
                                                                   caacaaCTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                              >  3:538775/1‑149 (MQ=255)
                                                                          ttCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAc                                                                       <  1:126208/149‑1 (MQ=255)
                                                                            cGCTATATTGTTGAGGTGGTCAATCGTAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCt                                                                     >  3:512798/1‑149 (MQ=255)
                                                                              cTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGa                                                                   <  3:316988/149‑1 (MQ=255)
                                                                               tataTTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       >  4:509345/1‑112 (MQ=255)
                                                                               tataTTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       <  2:395968/112‑1 (MQ=255)
                                                                               tataTTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       >  1:395968/1‑112 (MQ=255)
                                                                               tataTTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAg                                                                                                       <  3:509345/112‑1 (MQ=255)
                                                                                 tatTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGc                                                                >  3:52386/1‑149 (MQ=255)
                                                                                        gAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAAc                                                         >  3:389461/1‑149 (MQ=255)
                                                                                        gAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAAc                                                         <  4:435856/149‑1 (MQ=255)
                                                                                          ggtggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGc                                                       >  2:111564/1‑149 (MQ=255)
                                                                                             ggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                                    >  1:54440/1‑149 (MQ=255)
                                                                                             ggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                                    >  3:62619/1‑149 (MQ=255)
                                                                                             ggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                                    >  3:119194/1‑149 (MQ=255)
                                                                                                 aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAgg                                                >  2:226620/1‑149 (MQ=255)
                                                                                                 aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAgg                                                <  3:380790/149‑1 (MQ=255)
                                                                                                  aTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAggg                                               >  2:424794/1‑149 (MQ=255)
                                                                                                       aaCCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGCATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATACAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGa                                          <  4:389461/149‑1 (MQ=255)
                                                                                                         ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGaaa                                        >  4:96476/1‑149 (MQ=255)
                                                                                                         ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGaaa                                        <  2:428791/149‑1 (MQ=255)
                                                                                                            gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                                     <  3:115320/149‑1 (MQ=255)
                                                                                                            gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                                     <  3:96476/149‑1 (MQ=255)
                                                                                                               tGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAACTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                                  <  3:305972/149‑1 (MQ=255)
                                                                                                                   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGTTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAt                              >  3:327120/1‑149 (MQ=255)
                                                                                                                   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAt                              >  1:137489/1‑149 (MQ=255)
                                                                                                                   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAt                              >  1:187598/1‑149 (MQ=255)
                                                                                                                   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAt                              >  2:567103/1‑149 (MQ=255)
                                                                                                                                   gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTg              >  2:514214/1‑149 (MQ=255)
                                                                                                                                   gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTg              <  1:351477/149‑1 (MQ=255)
                                                                                                                                     cTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTc            <  2:46656/149‑1 (MQ=255)
                                                                                                                                         acacATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGaaa        <  4:538441/149‑1 (MQ=255)
                                                                                                                                         acacATCACAACCCAGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGaca        <  4:358032/149‑3 (MQ=255)
                                                                                                                                           acaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGt      >  1:40917/1‑149 (MQ=255)
                                                                                                                                           aaaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGt      >  1:463739/3‑149 (MQ=255)
                                                                                                                                             aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTGACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCg    >  3:20841/1‑149 (MQ=255)
                                                                                                                                             aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCg    >  4:288009/1‑149 (MQ=255)
                                                                                                                                             aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCg    >  1:509573/1‑149 (MQ=255)
                                                                                                                                             aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCg    >  3:109093/1‑149 (MQ=255)
                                                                                                                                              tCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGa   >  4:254096/1‑149 (MQ=255)
                                                                                                                                               cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAGGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGAt  <  2:166154/149‑1 (MQ=255)
                                                                                                                                               cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGAt  >  2:333849/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGTGGATGGTTTAACATGAAATTACAACAACTTCGCTATATTGTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGAT  >  NC_000913/1333797‑1334088

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