Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,004,613 G→A 11.3% intergenic (+60/+90) PP_1787 → / ← PP_1788 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,004,6130GA11.3% 189.0 / 12.6 71intergenic (+60/+90)PP_1787/PP_1788hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (4/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TTTGCGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCTTT  >  NC_002947/2004464‑2004759
                                                                                                                                                     |                                                                                                                                                  
tttGCGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAg                                                                                                                                                    >  2:1919578/1‑150 (MQ=255)
tttGCGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAg                                                                                                                                                    <  1:848904/150‑1 (MQ=255)
       ccGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCa                                                                                                                                             >  2:853610/1‑150 (MQ=255)
       ccGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATACCCCAAAAATTAAATAGCTTATTGCAGATTTGCa                                                                                                                                             >  2:803495/1‑150 (MQ=255)
         gTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAAt                                                                                                                                           >  1:1368546/1‑150 (MQ=255)
            ttttCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAg                                                                                                                                        <  1:50044/150‑1 (MQ=255)
              ttCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAGCt                                                                                                                                      <  2:242402/150‑1 (MQ=255)
                      ttCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAGCTATTTcatt                                                                                                                              <  1:853610/150‑4 (MQ=255)
                            gTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGt                                                                                                                        <  1:182334/150‑1 (MQ=255)
                            gTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGAAATACGCCCTTTCTGAATTTGt                                                                                                                        >  2:136695/1‑150 (MQ=255)
                             ttCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCaa                                                                                                                                            >  2:1061418/1‑129 (MQ=255)
                             ttCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCaa                                                                                                                                            <  1:1061418/129‑1 (MQ=255)
                                   gTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGa                                                                                                                 >  2:1605974/1‑150 (MQ=255)
                                     cGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAgg                                                                                                               <  2:652782/150‑1 (MQ=255)
                                       aaaTTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCg                                                                                                             >  2:1485052/1‑150 (MQ=255)
                                                   gAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTg                                                                                                 <  2:382609/150‑1 (MQ=255)
                                                     gTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGa                                                                                               <  1:849228/150‑1 (MQ=255)
                                                     gTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGa                                                                                               <  1:88841/150‑1 (MQ=255)
                                                          ttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGa                                                                                          <  2:1368546/150‑1 (MQ=255)
                                                             aCGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAAc                                                                                       <  1:2382098/150‑1 (MQ=255)
                                                                        aTGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATTTGAAACCGCTCGCCTCa                                                                            >  2:990658/1‑150 (MQ=255)
                                                                           aCGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGc                                                                         >  1:1831223/1‑150 (MQ=255)
                                                                             gcttaatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCtt                                                                       >  2:1776964/1‑150 (MQ=255)
                                                                             gcttaatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGATATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCtt                                                                       >  1:418799/1‑150 (MQ=255)
                                                                              cttaatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCttt                                                                      <  2:1831223/150‑1 (MQ=255)
                                                                              cttaatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCttt                                                                      <  1:803495/150‑1 (MQ=255)
                                                                                  atgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGa                                                                  <  2:1460615/150‑1 (MQ=255)
                                                                                   tgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                 <  2:418799/150‑1 (MQ=255)
                                                                                   tgattaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTCCAGATTTGCAATAAGCTATTTCTCAATTTGTTAGCTGAGGCGCTCGGTTTTTTGGATGTGAAACACCTCGCCTCAGGCCTTGTGAt                                                                 >  1:2401847/4‑150 (MQ=255)
                                                                                     cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTg                                                                                                                                               >  2:2111480/1‑70 (MQ=255)
                                                                                     cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTg                                                                                                                                               <  1:2111480/70‑1 (MQ=255)
                                                                                     cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgt                                                               <  1:1672781/150‑1 (MQ=255)
                                                                                      ttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgtg                                                              >  2:1214861/1‑150 (MQ=255)
                                                                                          taaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTAtt                                                          >  1:1316395/1‑150 (MQ=255)
                                                                                            aGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTa                                                        >  1:882011/1‑150 (MQ=255)
                                                                                             gATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTAt                                                       <  1:1214861/150‑1 (MQ=255)
                                                                                              atatTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTAATTATg                                                      >  2:42567/1‑150 (MQ=255)
                                                                                                          tACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgt                                                 <  2:1073837/143‑1 (MQ=255)
                                                                                                          tACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgt                                                 >  1:1073837/1‑143 (MQ=255)
                                                                                                             cTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTaa                                       >  2:185708/1‑150 (MQ=255)
                                                                                                              tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgt                                                               <  1:1303055/125‑1 (MQ=255)
                                                                                                              tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgt                                                               >  2:1303055/1‑125 (MQ=255)
                                                                                                              tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                 >  2:1521478/1‑123 (MQ=255)
                                                                                                              tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                 <  1:1521478/123‑1 (MQ=255)
                                                                                                              tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAg                                      <  1:1485052/150‑1 (MQ=255)
                                                                                                               ccccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGc                                     <  1:136695/150‑1 (MQ=255)
                                                                                                                cccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCt                                    >  2:1179813/1‑150 (MQ=255)
                                                                                                                 ccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCtt                                   <  1:1179813/150‑1 (MQ=255)
                                                                                                                  cTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgtgt                                               <  2:50069/137‑1 (MQ=255)
                                                                                                                  cTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgtgt                                               >  1:50069/1‑137 (MQ=255)
                                                                                                                     gACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCt                                                                                                                                <  2:1696593/53‑1 (MQ=255)
                                                                                                                     gACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCt                                                                                                                                >  1:1696593/1‑53 (MQ=255)
                                                                                                                         tGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGNCTTTGTNATGTGTATTTATNATAGTGTTGATCTAAGCTTGCAATACg                           <  1:1605974/150‑1 (MQ=255)
                                                                                                                              cGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTg                      <  1:1180779/150‑1 (MQ=255)
                                                                                                                                      aaaTAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                 <  2:463623/99‑1 (MQ=255)
                                                                                                                                      aaaTAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                 >  1:463623/1‑99 (MQ=255)
                                                                                                                                      aaaTAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAAt              >  1:2238719/1‑150 (MQ=255)
                                                                                                                                       aaTAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAg                                      <  1:2179247/125‑1 (MQ=255)
                                                                                                                                       aaTAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAg                                      >  2:2179247/1‑125 (MQ=255)
                                                                                                                                         tAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGt           <  1:856676/150‑1 (MQ=255)
                                                                                                                                           gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgt                                                 >  1:1307050/1‑110 (MQ=255)
                                                                                                                                           gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgt                                                 <  2:1307050/110‑1 (MQ=255)
                                                                                                                                            cTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTTAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGcc        >  2:1017766/1‑150 (MQ=255)
                                                                                                                                            cTTATTGCAAATCTGCAATAAGCTATTTCTAAATTTGTTAGCTGAgg                                                                                                               >  1:170183/1‑47 (MQ=18)
                                                                                                                                            cTTATTGCAAATCTGCAATAAGCTATTTCTAAATTTGTTAGCTGAgg                                                                                                               <  2:170183/47‑1 (MQ=18)
                                                                                                                                            cTTATTGCAAATCTGCAATAAGCTATTTCTAAATGTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACTGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATAGGGGTTGGCTTTAgtggtgtc        >  1:316922/1‑148 (MQ=255)
                                                                                                                                             ttATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATATGGATTGGTTTTAATGGTGCCa       >  1:865474/1‑150 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGt           <  1:1314906/141‑1 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGt           >  2:1314906/1‑141 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCttt  >  2:1970085/1‑150 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCttt  <  1:1919578/150‑1 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGAAATACGGATTGGTTTTAATGGTg          <  2:2085812/142‑1 (MQ=255)
                                                                                                                                                  gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGAAATACGGATTGGTTTTAATGGTg          >  1:2085812/1‑142 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTGCGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCTTT  >  NC_002947/2004464‑2004759

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