Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,624,859 C→T 5.7% I29I (ATC→ATT ydeJ → inactive PncC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,8590CT 5.7% 220.2 / 2.7 70I29I (ATC→ATTydeJinactive PncC family protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base T (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

ATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGC  >  NC_000913/1624712‑1625003
                                                                                                                                                   |                                                                                                                                                
atTGGACGAGCTTACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATcg                                                                                                                                                 >  8:254718/1‑149 (MQ=255)
           aTACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTa                                                                                                                                      <  8:328837/149‑1 (MQ=255)
            tACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTAc                                                                                                                                     >  3:95986/1‑149 (MQ=255)
            tACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTAAGTCAAAGACTTATCGCGGATCAATTAc                                                                                                                                     <  3:53610/149‑1 (MQ=255)
              cAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGc                                                                                                                                   <  7:237621/149‑1 (MQ=255)
               aGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCTGATCAATTACGCt                                                                                                                                  >  8:111709/1‑149 (MQ=255)
               aGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCt                                                                                                                                  >  5:342724/1‑149 (MQ=255)
                   aaaGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAAc                                                                                                                              <  5:373455/149‑1 (MQ=255)
                            cACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTAccgccg                                                                                                                      >  8:162814/1‑148 (MQ=255)
                             aCCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGa                                                                                                                                             <  7:280357/124‑1 (MQ=255)
                             aCCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGa                                                                                                                                             >  8:280357/1‑124 (MQ=255)
                                tttACATATAAAAGCGTCGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATc                                                                                                                  >  8:305632/1‑148 (MQ=255)
                                         ataaGCTTGGACAATTATCAGTGATCATTAATAGTGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGTCATCCGTGTTGAGTCAATGACTTATCGCGGATCATTTAGGCTTTACTACCGCCGAATCATGCACcgt                                                                                                         <  4:337961/146‑2 (MQ=255)
                                         aaaaGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcggc                                                                                                        >  8:275650/1‑149 (MQ=255)
                                          aaaGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcggcg                                                                                                       >  5:210701/1‑149 (MQ=255)
                                                    cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGc                                                                                                                                   >  4:251635/1‑111 (MQ=255)
                                                    cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGc                                                                                                                                   <  3:251635/111‑1 (MQ=255)
                                                    cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGc                                                                                             <  7:311929/149‑1 (MQ=255)
                                                      aaaaaCAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGc                                                                                             >  1:17917/1‑147 (MQ=255)
                                                          aCAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATcg                                                                                                                                                 >  1:337712/1‑91 (MQ=255)
                                                          aCAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATcg                                                                                                                                                 <  2:337712/91‑1 (MQ=255)
                                                            aGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATc                                                                                                                  <  6:345323/120‑1 (MQ=255)
                                                            aGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATc                                                                                                                  >  5:345323/1‑120 (MQ=255)
                                                             gTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAg                                                                                                  <  3:141027/135‑1 (MQ=255)
                                                             gTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAg                                                                                                  >  4:141027/1‑135 (MQ=255)
                                                              tGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtg                                                                                   >  8:269750/1‑149 (MQ=255)
                                                               gAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgt                                                                                  >  4:127014/1‑149 (MQ=255)
                                                                aaCATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTg                                                                                               <  8:225266/135‑1 (MQ=255)
                                                                aaCATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTg                                                                                               >  7:225266/1‑135 (MQ=255)
                                                                aaCATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgt                                                                                  <  6:306890/148‑1 (MQ=255)
                                                                  cATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgtgtg                                                                               <  7:275650/149‑1 (MQ=255)
                                                                  cATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgtgtg                                                                               <  7:254718/149‑1 (MQ=255)
                                                                     tAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAg                                                                            >  7:262682/1‑149 (MQ=255)
                                                                          gagaTAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGaa                                                                       <  2:17917/149‑1 (MQ=255)
                                                                          gagaTAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGa                                                                        <  2:146639/148‑1 (MQ=255)
                                                                                aaaTTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGGTGAGGCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATaca                                                                 <  3:125289/149‑1 (MQ=255)
                                                                                      tACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCaa                                                            <  4:368202/148‑1 (MQ=255)
                                                                                            ttGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCGGTGTGCAGCTGAAGATACACCCAAATTTTAc                                                     >  2:285586/1‑149 (MQ=255)
                                                                                              gTAGATACCGATACTATTGAAATCCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTAc                                                     <  6:342724/147‑1 (MQ=255)
                                                                                              gTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACgg                                                   >  4:175404/1‑149 (MQ=255)
                                                                                                 gATACCGAGACTAGTGAAAACCTGACATCCGCGGTGAGTCAAAGACTTATCGCGGATCAAGTACGCTTAACTACCGCCGCCTCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCCAATTTTACGGTg                                                 >  4:6648/1‑148 (MQ=255)
                                                                                                   tACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAg                                              <  7:68973/149‑1 (MQ=255)
                                                                                                         tACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTg                                        <  7:305632/149‑1 (MQ=255)
                                                                                                            tATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATTGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGtt                                      >  1:199630/1‑148 (MQ=255)
                                                                                                            tATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATTGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTa                                     >  1:205758/1‑149 (MQ=255)
                                                                                                            tATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTa                                     >  3:266186/1‑149 (MQ=255)
                                                                                                              ttGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACt                                   >  4:332243/1‑149 (MQ=255)
                                                                                                                gAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACtt                                  <  6:210701/148‑1 (MQ=255)
                                                                                                                 aaaaCCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACttt                                 >  2:257133/1‑148 (MQ=255)
                                                                                                                    aCCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCAcc                             <  7:269750/149‑1 (MQ=255)
                                                                                                                         aCATCCGCGTTGAGTCAAAGACTTATTGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACTGTGCAGGCTTTGTTACTTTCACCGATc                         <  2:199630/148‑1 (MQ=255)
                                                                                                                         aCATCCGCGTTGAGTCAAAGACTTATTGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATc                         <  2:205758/148‑1 (MQ=255)
                                                                                                                         aCATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATc                         <  2:428228/148‑1 (MQ=255)
                                                                                                                           aTCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCaaa                                                           >  8:250303/1‑112 (MQ=255)
                                                                                                                           aTCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCaaa                                                           <  7:250303/112‑1 (MQ=255)
                                                                                                                              cgcgTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCaa                   <  7:162814/149‑1 (MQ=255)
                                                                                                                                      gTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGaa            >  7:162191/1‑148 (MQ=255)
                                                                                                                                       tCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTGTGTTACTTTCACCGATCAGGCAAAGATGaaaa          >  1:334550/1‑149 (MQ=255)
                                                                                                                                        cAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACg                                                    <  2:120088/106‑1 (MQ=255)
                                                                                                                                        cAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACg                                                    >  1:120088/1‑106 (MQ=255)
                                                                                                                                        cAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAAt         >  1:234392/1‑149 (MQ=255)
                                                                                                                                         aaaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCt                                                                                    >  3:196040/1‑73 (MQ=255)
                                                                                                                                         aaaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCt                                                                                    <  4:196040/73‑1 (MQ=255)
                                                                                                                                         aaaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATc        >  3:342475/1‑149 (MQ=255)
                                                                                                                                         aaaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATc        <  3:127014/149‑1 (MQ=255)
                                                                                                                                          aaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGa                <  7:400462/140‑1 (MQ=255)
                                                                                                                                          aaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGa                >  8:400462/1‑140 (MQ=255)
                                                                                                                                            gACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCaaa                  <  4:203442/136‑1 (MQ=255)
                                                                                                                                            gACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCaaa                  >  3:203442/1‑136 (MQ=255)
                                                                                                                                               ttATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGc  >  3:194826/1‑149 (MQ=255)
                                                                                                                                                 aTCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGc  <  8:262682/147‑1 (MQ=255)
                                                                                                                                                  tCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                            <  1:243684/56‑1 (MQ=255)
                                                                                                                                                  tCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                            >  2:243684/1‑56 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGC  >  NC_000913/1624712‑1625003

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