Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 949,117 (T)5→4 coding (452/909 nt) ycaN ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913949,1130T.91.9% 344.0 / 16.7 87coding (456/909 nt)ycaNLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base . (39/40);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTG  >  NC_000913/948964‑949258
                                                                                                                                                     |                                                                                                                                                  
aaaGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATTGCCACGCCTATCCTAtttt                                                                                                                                                   >  1:806984/1‑151 (MQ=255)
       tGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCttttcgac                                                                                                                                             <  2:76341/150‑5 (MQ=255)
       tGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCttttcgac                                                                                                                                             <  2:922128/150‑5 (MQ=255)
              gggATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCaa                                                                                                                                     <  2:962012/150‑1 (MQ=255)
              gggATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACCCTGATCATATCTTTTTGACGATGCaa                                                                                                                                     >  1:631425/1‑150 (MQ=255)
                gATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTa                                                                                                                                  <  1:350000/151‑1 (MQ=255)
                 aTAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTTGACGATGCAACTa                                                                                                                                  >  1:616066/1‑150 (MQ=255)
                    aCGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGTCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTc                                                                                                                              >  2:629027/1‑151 (MQ=255)
                    aCGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTc                                                                                                                              >  2:951535/1‑151 (MQ=255)
                     cGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCg                                                                                                                             <  1:453987/151‑1 (MQ=255)
                      gAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTTGACGATGCAACTAAGTCg                                                                                                                             >  1:749012/1‑150 (MQ=255)
                      gAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCg                                                                                                                             >  2:659001/1‑150 (MQ=255)
                      gAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGa                                                                                                                            <  2:1116328/151‑1 (MQ=255)
                        aaaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAAc                                                                                                                          <  2:963455/151‑1 (MQ=255)
                        aaaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAAc                                                                                                                          <  1:1014456/151‑1 (MQ=255)
                        aaaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAAc                                                                                                                          <  2:695993/151‑1 (MQ=255)
                        aaaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCCCCACGCTGATCATATCTTTTCCACGATGCAACCAAGTCGaa                                                                                                                           >  1:1017496/1‑150 (MQ=255)
                         aaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGaa                                                                                                                           >  2:45460/1‑149 (MQ=255)
                         aaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAAc                                                                                                                          >  2:1181125/1‑150 (MQ=255)
                         aaaCAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACg                                                                                                                         <  1:737434/151‑1 (MQ=255)
                          aacaacACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGaa                                                                                                                           >  2:545223/1‑148 (MQ=255)
                             aacaCATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAg                                                                                                                     <  2:1007804/151‑1 (MQ=255)
                                acaTTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCg                                                                                                                   >  2:85792/1‑150 (MQ=255)
                                acaTTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGt                                                                                                                 <  2:336047/151‑1 (MQ=255)
                                  aTTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCCCGCTGATCATATCTTTTTGACGATGCAACTAAGTCGAACGCCAACg                                                                                                                   >  1:952659/1‑148 (MQ=255)
                                      atgatTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTTGACGATGCCACTAAGTCGAACGCCAGCGTCaaa                                                                                                             >  2:722246/1‑149 (MQ=255)
                                           tAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCCCGCTGGTCATATCTTTTTGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGccc                                                                                                         >  1:13641/1‑148 (MQ=255)
                                            aaGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGtt                                                                                                     >  2:661447/1‑151 (MQ=255)
                                                 gATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACg                                                                                                >  1:451174/1‑151 (MQ=255)
                                                  aTCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACg                                                                                                >  2:788000/1‑150 (MQ=255)
                                                  aTCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGa                                                                                               <  2:257833/151‑1 (MQ=255)
                                                   tCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGGACGCCACGCTTTACCTATCTTTTCGACCATGGCACTAAGTCTGACGCCCACGTCAAAGACCCGTTGGAc                                                                                                 >  1:1152041/1‑148 (MQ=255)
                                                    cATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGat                                                                                              >  1:808593/1‑150 (MQ=255)
                                                       gAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCg                                                                                          <  1:441260/151‑1 (MQ=255)
                                                          gATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACa                                                                                       <  1:959925/151‑1 (MQ=255)
                                                             ggCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAag                                                                                     <  1:906445/150‑1 (MQ=255)
                                                              gCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAagag                                                                                   <  1:709753/151‑1 (MQ=255)
                                                                  ggTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAgtc                                                                               <  1:1155727/151‑1 (MQ=255)
                                                                    tttGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAgtcg                                                                              >  2:144787/1‑150 (MQ=255)
                                                                        ccATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATATTTTCGACGATGCAACTAAGGCGAACGCCATCGTCAAAGCCCTGTTGGACGATATCGGCAAGAGAGGCGTCAg                                                                          >  2:1201340/1‑150 (MQ=255)
                                                                         cATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGtt                                                                        <  1:977451/151‑1 (MQ=255)
                                                                         cATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGtt                                                                        <  1:98394/151‑1 (MQ=255)
                                                                          atatCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGGCGATCGCCACCCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCCGTTTGACGATATCGACAAGAGAGTCGTCCGtt                                                                        >  1:809261/1‑150 (MQ=255)
                                                                             tCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGACCGCCACGCCTATCATATCTTTTTGGAGATTCAACTTCATCGGACGCCAGCGCTCAATCCCCTTTTTGCGGTATACACAAAAGAGTCCTCCCTtgtg                                                                     >  1:577920/1‑150 (MQ=255)
                                                                              cTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTtgtg                                                                     >  2:178373/1‑149 (MQ=255)
                                                                                 gCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCt                                                                 <  2:661739/150‑1 (MQ=255)
                                                                                    aaGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGAc                                                             <  1:179874/151‑1 (MQ=255)
                                                                                    aaGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGACGAACGCCACGCTGGTCATATCCTTTCCAAGAAGCAACTAAATCGAACGCCAGCGTCAAAGCCCTGTTGGACGAAATCGACAAGAAAGTCGTCAGGTGTGAGCCCGa                                                              >  1:50988/1‑150 (MQ=255)
                                                                                     aGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCCATCGCCACGCTGCGCCAATCTTTTCGACGAATCTAAAAAATCGAACGCCAGCGCCAAAGCCCTCTTTGTGGATATCGACAATAGAGTCGTCAGTTGTGGGGTCgcgt                                                            >  1:618488/1‑148 (MQ=255)
                                                                                      gTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACtt                                                           >  2:741334/1‑151 (MQ=255)
                                                                                      gTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACtt                                                           <  2:675952/151‑1 (MQ=255)
                                                                                             ggCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACCCTGATCATATCTTTTCGACGATGCAACAAAGTCGAAAGCCAGCGTCAAAACCCTCTTTGACGATTTCGGAAAGAGAGTCGTCTGTTGTGAGGTCGACTTTTatat                                                     >  2:569729/1‑150 (MQ=255)
                                                                                              gCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCa                                                   <  1:341110/151‑1 (MQ=255)
                                                                                                 ttGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAATCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAAAGAGTCGTCAGTTGTGAGCTCGACTTTTATATCAgg                                                 >  2:143657/1‑150 (MQ=255)
                                                                                                   gctgctACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGat                                               >  2:209823/1‑150 (MQ=255)
                                                                                                   gctgctACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGat                                               <  1:464514/150‑1 (MQ=255)
                                                                                                   gctgctACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGat                                               <  2:911620/150‑1 (MQ=255)
                                                                                                      gctACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGatat                                             <  1:376688/149‑1 (MQ=255)
                                                                                                      gctACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATAtt                                            <  1:579726/150‑1 (MQ=255)
                                                                                                           acaTAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACg                                        >  2:1187834/1‑149 (MQ=255)
                                                                                                               aaTTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGaa                                  <  2:35199/151‑1 (MQ=255)
                                                                                                                aTTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGaaa                                 <  1:226403/151‑1 (MQ=255)
                                                                                                                        ggCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAaaccaacc                           <  2:800156/149‑1 (MQ=255)
                                                                                                                        ggCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAaaccaacc                           >  2:1003840/1‑149 (MQ=255)
                                                                                                                        ggCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAg                         <  2:736530/151‑1 (MQ=255)
                                                                                                                        ggCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAg                         <  1:370709/151‑1 (MQ=255)
                                                                                                                         gCGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCGACCAGt                        <  2:917107/151‑1 (MQ=255)
                                                                                                                          cGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGt                        >  2:972682/1‑150 (MQ=255)
                                                                                                                          cGGGCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTa                       <  1:609751/151‑1 (MQ=255)
                                                                                                                            ggCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTaaa                     <  2:635551/151‑1 (MQ=255)
                                                                                                                             gCCGATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAg                    <  1:478874/151‑1 (MQ=255)
                                                                                                                                gATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGAc                  >  1:1007123/1‑150 (MQ=255)
                                                                                                                                gATCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACa                 <  2:1183607/151‑1 (MQ=255)
                                                                                                                                 aTCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACAt                >  2:234696/1‑151 (MQ=255)
                                                                                                                                 aTCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACAt                <  2:862370/151‑1 (MQ=255)
                                                                                                                                  tCGCCACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACAtt               >  1:156434/1‑151 (MQ=255)
                                                                                                                                     ccACGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTaa             >  2:176964/1‑150 (MQ=255)
                                                                                                                                        cGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAat         >  2:503667/1‑151 (MQ=255)
                                                                                                                                        cGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAat         >  1:1144394/1‑151 (MQ=255)
                                                                                                                                        cGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAat         >  1:552215/1‑151 (MQ=255)
                                                                                                                                        cGCTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAat         >  2:1062334/1‑151 (MQ=255)
                                                                                                                                          cTGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAatat       >  2:719958/1‑151 (MQ=255)
                                                                                                                                           tGATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAatat       >  1:35715/1‑150 (MQ=255)
                                                                                                                                            gATCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATAc     >  1:299327/1‑151 (MQ=255)
                                                                                                                                             aTCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACg    >  2:853919/1‑151 (MQ=255)
                                                                                                                                             aTCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACg    >  2:1116121/1‑151 (MQ=255)
                                                                                                                                              tCATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGt   >  2:317236/1‑151 (MQ=255)
                                                                                                                                               cATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGta  <  1:195243/151‑2 (MQ=255)
                                                                                                                                               cATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGt   <  1:348178/150‑1 (MQ=255)
                                                                                                                                               cATATCTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTg  <  2:83921/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGCCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTG  >  NC_000913/948964‑949258

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