Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,800 C→A 5.7% A580S (GCG→TCG)  yicI ← alpha‑D‑xyloside xylohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,8000CA 5.7% 205.6 / 3.0 80A580S (GCG→TCG) yicIalpha‑D‑xyloside xylohydrolase
Reads supporting (aligned to +/- strand):  ref base C (27/48);  new base A (3/2);  total (30/50)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

TGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGA  >  NC_000913/3832658‑3832940
                                                                                                                                               |                                                                                                                                            
tGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCCTCGGCGTACCCCGCGCGTTCGc                                                                                                                                              >  1:549302/1‑143 (MQ=255)
   acaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGTAACTCCATGATCATGGCACGCATCATCGGCATACCCAGCGCATTCGCAcg                                                                                                                                           >  1:211441/1‑143 (MQ=255)
   acaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCCCGCGTTCGCAcg                                                                                                                                           >  2:254922/1‑143 (MQ=255)
     aTCGCGCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACGAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgcg                                                                                                                                         <  2:383833/143‑1 (MQ=255)
        gcgcGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAg                                                                                                                                      <  2:382845/140‑1 (MQ=255)
        gCCGGCTGCAGGGAACACCGGCGCAGTCATCACGTTGTGGGCTGACATGTATGGACGGCGAAGGTAATCACAAGCCGGATCGTCCGGGAACGCCATCATCATGGCGCGCATCATCGGGGTACCGCGCGCGTGCGAACgcgcgg                                                                                                                                      <  2:347774/143‑3 (MQ=255)
        gCCCGCTTCGGTGAACACCGGCGCCGAGTATGCGTTGTCGCCTAAGATGTCTGGACGGTCAAGGTAATCACAAGCCGGATCGTCCGCGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAg                                                                                                                                      <  2:202333/143‑1 (MQ=255)
         cccGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCCGGAACTCCATCATCATGGCCCGCATCATCGGCCTACCCCGCGCGTTCGCACGCGCAGc                                                                                                                                     >  2:153556/1‑143 (MQ=255)
         cccGCTCAAGTTAACGCCTGCGCAACCATCACGTTGCCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                                                                     <  1:300751/143‑4 (MQ=255)
           cGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCtt                                                                                                                                   >  1:506540/1‑143 (MQ=255)
             cTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                 <  2:67898/143‑1 (MQ=255)
              ttCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGCACTCCATCATCATGGCCCGCATCCTCGGCGTACCCCGCGCGTTCGCCCGCGCAGCTTCAc                                                                                                                                >  1:420658/1‑143 (MQ=255)
                  gTGAACGCCGGCGCACCCATCGCGTTGTCGCCTAACAGCTGTGGACGGTCAAGGTAATCACAAGCCGGATCGCCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                                                            <  2:810195/143‑1 (MQ=255)
                    gAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACa                                                                                                                          <  2:724836/143‑1 (MQ=255)
                    gAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTACGAGAGCGCGGCTCCGCGGTACa                                                                                                                          >  1:301163/1‑143 (MQ=255)
                      acacCGGCGCAACCCTCACGTTGTCGCCTCACGTGTAGTGGCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGa                                                                                                                        <  2:242176/143‑1 (MQ=255)
                          cgccgcAACCAGCACGTTGTCGCCCAACATGTAGTGACGGCCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCAGCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                                                                    <  1:105534/140‑1 (MQ=255)
                          cGGCGCAACCGTCACGCTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACGCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                                                                    <  2:372989/143‑1 (MQ=255)
                           ggCGCAACCATCACTTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACgg                                                                                                                   <  2:442316/143‑1 (MQ=255)
                                 aaCCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatca                                                                                                              <  1:813214/143‑1 (MQ=255)
                                     aTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCCGCTTCCCGATGCCGCTACGGCACCATGCg                                                                                                          >  1:34242/1‑143 (MQ=255)
                                     aTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCg                                                                                                          <  1:797654/143‑1 (MQ=255)
                                               cggctcaCATGTATTGACGAATATGGTTATCACGAGCCGGGTCGTCCGGGACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                                                                <  2:753569/137‑1 (MQ=255)
                                               cGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                                                                <  1:731434/143‑1 (MQ=255)
                                               cGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                                                                <  1:175701/143‑1 (MQ=255)
                                                gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGCACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGCTACAGATACGGCATCATGCGGCATTTCAGtt                                                                                               >  1:200085/1‑143 (MQ=255)
                                                  cTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                             <  2:153002/143‑1 (MQ=255)
                                                    accATGTCTGGACTGTCAAGGTAATCTCACGCCGGATCGTCCGGGACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                           <  2:641941/141‑1 (MQ=255)
                                                    aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                           <  1:265493/143‑1 (MQ=255)
                                                    aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                           <  1:141124/143‑1 (MQ=255)
                                                      cATGTATGGACGGTCAAGGGAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTga                                                                                         <  1:210120/143‑1 (MQ=255)
                                                       aTGTACTCACGCTGAATGTAATCACAGTCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCATCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaa                                                                                        <  2:659658/143‑1 (MQ=255)
                                                                  gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCCTCCTCATGGCCCGCATCATCGGCGTACCCCACGCCTTCGCACGCGCAGCTACACGAAACAGAAACGTCATCATGCGGCATTTCAGTTGCGTGCAGAAGCGcacca                                                                             >  1:652914/1‑143 (MQ=255)
                                                                     aaGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACTCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACAt                                                                          <  1:413993/143‑1 (MQ=255)
                                                                      aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATc                                                                         >  2:78984/1‑143 (MQ=255)
                                                                         tagtCACACGCCGGATCGTCCGGGAGCTCGAGCATCGTGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                      <  2:226726/140‑1 (MQ=255)
                                                                         tAATCACACGCCGGACCGTCCGGGATCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                      <  2:825741/143‑1 (MQ=255)
                                                                         tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                      <  1:392837/143‑1 (MQ=255)
                                                                         caaTCACAAGCCGGATCGTCCGGGGACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGCACGCGCAGCTGCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                      <  1:555081/142‑1 (MQ=255)
                                                                          aaTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAg                                                                     <  2:476937/143‑1 (MQ=255)
                                                                           aTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTCACCCCCTCGTTCGCACCCGCAGCCTCACGATCCCCATACGGCATCATACGGCACCTCAGTACCGCGCAGCAGCGCACCAAATCAGAgg                                                                    >  2:738812/1‑143 (MQ=255)
                                                                            tCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                   <  1:717900/143‑1 (MQ=255)
                                                                            tCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTAGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGACAAAGCGCAACACATAACACGa                                                                   >  2:641587/1‑143 (MQ=255)
                                                                              acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                                                                 >  2:389020/1‑143 (MQ=255)
                                                                               caAGCCGGATCGTCCGGGAACTCCATCATCATGGTCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATCTCAGTTGCGAGAAGACCAGCACCACATCACAGGACtc                                                                >  2:662045/1‑143 (MQ=255)
                                                                                aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAAtca                                                               >  2:449313/1‑143 (MQ=255)
                                                                                     ggATCGTCCGGGACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCatcatc                                                          <  2:313753/143‑1 (MQ=255)
                                                                                         cctCCGGCACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                      <  2:374322/141‑1 (MQ=255)
                                                                                         cGTCCGGGAGCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                      <  2:519092/143‑1 (MQ=255)
                                                                                         cGTCCGGGACGTCCATCTTCATGGCCCGCATCATCGGCCCCCCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                      <  2:581354/143‑1 (MQ=255)
                                                                                         cGTCCGGGAACTTCATTATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTGACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                      <  2:382265/143‑1 (MQ=255)
                                                                                          gTCCGGGAACGCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg                                                     <  1:696610/143‑1 (MQ=255)
                                                                                          cgccGGGATCTCCATCATCATGGCCCGCATCACCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTGCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg                                                     <  1:121925/141‑1 (MQ=255)
                                                                                           tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCCTGCGGCATTTCAGCTGCGCGAAAACGCGCCCCCCATCACAGGAATAATGATCCTAGGc                                                    >  2:791654/1‑143 (MQ=255)
                                                                                           tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGGAGGc                                                    >  1:823355/1‑143 (MQ=255)
                                                                                            ccggtcaCTCCAGCACCGTGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGcc                                                   <  1:329716/137‑1 (MQ=255)
                                                                                            ccGGGAACTCCATCATCATGGCCCGCATCAGCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATCCAGATACGGCATCATGCGGCAGTTCAGTGGCGTGAGGGCGCGCACCCCATCCCAGGACTCATCATCGCATGcc                                                   >  1:21968/1‑143 (MQ=255)
                                                                                                gagcTCCATCATCATGGACCGCATCATCGGCGTACCCCGCGCGGTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAGTTCAGTGGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACg                                               <  1:90953/140‑1 (MQ=255)
                                                                                                gaccTCCATCATCATGGCCCGCATCACCGGCGTGCCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCCAATACTTGGACTCATTATCGTAGGCCCACg                                               <  2:722636/140‑1 (MQ=255)
                                                                                                  ccTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATCCGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGc                                             <  1:65876/142‑1 (MQ=255)
                                                                                                      catcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGCTACGGCATCATGCGGCATTTCAGTCGCGTGAAGCAGCGCACCACATCACAGGACTCCCCAGCGTCGGCCCACGGcagcc                                         >  1:315362/1‑140 (MQ=255)
                                                                                                        tcatcatGGCCCGCATCATCAGCGTACCCCGAGCGTGCGCACGCGCAGCTTCACGATACAGATACGGCATCATGAGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCACACGGCACCCGa                                       >  2:627107/1‑143 (MQ=255)
                                                                                                                gCCCGCCTCACCAGCGACCCCCGCGCGTTGGCACGCGCAGCTTCCCGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAttt                               <  2:752778/143‑1 (MQ=255)
                                                                                                                gCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCATCTTCACGATACAGATACGGCATCATCCGGCATTTAAGTTGCGAGAAGAAGCGCACCAAATCACAGGACTCATCACCATACGACCCCGGCACAAGATAGGAtat                               >  2:193343/1‑141 (MQ=255)
                                                                                                                        catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgt                       <  1:757624/143‑1 (MQ=255)
                                                                                                                        catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACACAACATCGTAGGCCCACGACACACGATACGACTCGCTACCgt                       >  2:400521/1‑143 (MQ=255)
                                                                                                                          tcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtgt                     <  2:475767/143‑1 (MQ=255)
                                                                                                                          tcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtgt                     <  1:449313/143‑1 (MQ=255)
                                                                                                                             gCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGACGCGCACCACATCACAGGACTCATCATCGTCGGCCCACGGCACACGATCAGACTTGCTGCCGTGTaca                  >  1:631374/1‑141 (MQ=255)
                                                                                                                             gCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACACCACCGGCCTCATCATCGTAGGCCCACGGCACACGATGAGAGTTGCTACCGTGTaaa                  >  2:624111/1‑143 (MQ=255)
                                                                                                                              cGTACCCCGCGCGTTCGCACGCGCAGCTGCACGATACAGATACGGCGTCATGCGGCATTTCAGTTGCGTGAAGAAGGGCACCACATCTCAGGACTCAGCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAAc                 <  2:23553/143‑1 (MQ=255)
                                                                                                                                tcccccGCGTGTTCGCACGCGCAGCCTCACGACACAGCTACGGCCTCATGCGGCATTTCAGTCGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg               <  1:552101/141‑1 (MQ=255)
                                                                                                                                tACCCCGCGCGTTCGCACTCGCAGCTCCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg               <  2:506540/143‑1 (MQ=255)
                                                                                                                                   gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCGCACGATAAGACTTGCTCCCGTGTACACGGCtc            >  1:470199/4‑142 (MQ=255)
                                                                                                                                    ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATCAGATTTGCTACCGTATAAACGGCTAt           >  2:317099/1‑143 (MQ=255)
                                                                                                                                    ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTGGCGACCGTGGAAACGGCTAt           >  2:5206/1‑143 (MQ=255)
                                                                                                                                     tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCAACAACGAAGGCACAGGGCACACGATAAGATTTGCTACCGTGTAAACAGCTaag          >  2:242983/2‑141 (MQ=255)
                                                                                                                                     cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACCCATCATCGTAGGCCCACGGCACACGATAAGATCTGCTACCGTGTAAACGGCGAtc          >  1:169602/1‑142 (MQ=255)
                                                                                                                                      gcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAATTGCTACCGTGTAAACGGCTATga         >  2:278282/1‑142 (MQ=255)
                                                                                                                                      gcccgTTCGCACGCTCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATgg         <  2:592642/140‑1 (MQ=255)
                                                                                                                                             cGCACGCGCAGCTTCTCTATACAGATACGGCATCATGCGGCATTTCAGTTGCTTGAAGAAGCGCACCACATCACATTATTTATCATCGTAGGCACACGGCAAATGATAAGATTTGCGACCGTGTAAACGGCTATGGCTGgaga  <  1:626742/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGA  >  NC_000913/3832658‑3832940

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