Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,696,909 G→A G151E (GGA→GAA)  bcsE → c‑di‑GMP‑binding protein BcsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,696,9090GA95.2% 138.5 / ‑1.1 42G151E (GGA→GAA) bcsEc‑di‑GMP‑binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (22/18);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

AATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGT  >  NC_000913/3696666‑3697147
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
aaTTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAaataat                                                                                                                                                                                                                                         >  1:158751/1‑251 (MQ=255)
             tcAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAAt                                                                                                                                                                                                                            >  2:519944/1‑251 (MQ=255)
                         aaaaaGAAAATTATTTTCAATGCAAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGAGTTCTAGTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATGAATCCCGAAAATAAGAACGATAAACAATTTTCATTGTTGc                                                                                                                                                                                                                <  2:76118/251‑1 (MQ=255)
                          aaaataaaatTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTTTTGCt                                                                                                                                                                                                               <  1:417832/251‑1 (MQ=255)
                            aataaaatTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGAGCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTg                                                                                                                                                                                                             <  2:147938/251‑1 (MQ=255)
                               aaaatTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAgg                                                                                                                                                                                                          >  1:315070/1‑251 (MQ=255)
                               aaaatTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAgg                                                                                                                                                                                                          <  1:358293/251‑1 (MQ=255)
                                                 aaaTCATGAAAAAGGTCTATACTATTTGACCCGTGAGTTGCAGTGTTCTATTGATCCCCATAATTTCCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCGGCCGAGCGGATTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGGTGCTTGAGGAATACCGTTCACttttt                                                                                                                                                                                         <  1:206076/250‑1 (MQ=255)
                                                   aTCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGt                                                                                                                                                                                      <  1:76690/251‑1 (MQ=255)
                                                               gTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGTAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCATTTTg                                                                                                                                                                          >  2:361859/1‑251 (MQ=255)
                                                                         tttGACCCGTGATATGCAGCGTTCTATTGTTCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATACCGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGCGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACACTTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAgg                                                                                                                                                                >  2:299265/1‑251 (MQ=255)
                                                                             aCCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTGCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGa                                                                                                                                                            <  1:571478/251‑1 (MQ=255)
                                                                                cGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGAcc                                                                                                                                                          >  1:113536/1‑250 (MQ=255)
                                                                                                   ttGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGcc                                                                                                                                      >  1:452310/1‑251 (MQ=255)
                                                                                                       tCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCt                                                                                                                                  >  1:577/1‑251 (MQ=255)
                                                                                                           cATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACAGTTGCTGGATATTGCCTTCTGGTg                                                                                                                              <  2:572520/251‑1 (MQ=255)
                                                                                                                    cTTTTTATTCTTTTTTGCGCAAATAACGCATGGCTCAACATTCCTGCCGAGCTGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGACGGTTAAACCATTGTTCGCTTTTGGAAATTAATCCCGAAAATAATAACGATAAACAATTTTCCTTGTTGCTTGAGGAATCCCGTTCTCTTTTTGGTCTTGCCAGTTTGCGTATTCAGGGTGACCAACTTTTGCTGGATATTGCCTTCTGGTGCAACGaaa                                                                                                                      >  2:364372/1‑250 (MQ=255)
                                                                                                                               tGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAgcg                                                                                                          >  1:70927/1‑251 (MQ=255)
                                                                                                                                ttttGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAgcg                                                                                                          <  2:452310/249‑1 (MQ=255)
                                                                                                                                       cAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTc                                                                                                    >  1:315696/1‑249 (MQ=255)
                                                                                                                                        aaaTAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTcagc                                                                                                 >  1:525813/1‑251 (MQ=255)
                                                                                                                                        aaaTAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGATTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTcagc                                                                                                 >  1:307495/1‑251 (MQ=255)
                                                                                                                                         aaTAACGCAGGGAAAAACATTCCTGCCGGGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGGGCAACGAAAAAGGGGTCAGCGCCCGTcagca                                                                                                <  2:241743/251‑1 (MQ=255)
                                                                                                                                                     cAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCTCTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCGTCTGGGGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCa                                                                                    <  1:326806/251‑1 (MQ=255)
                                                                                                                                                     cAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCa                                                                                    <  1:396370/251‑1 (MQ=255)
                                                                                                                                                       aaaCATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGc                                                                                  <  1:135837/251‑1 (MQ=255)
                                                                                                                                                                 gCCGAGCGGCTTCGCTCATGTTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATgg                                                                        <  1:527044/251‑1 (MQ=255)
                                                                                                                                                                 gCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAAt                                                                          <  1:161209/249‑1 (MQ=255)
                                                                                                                                                                  ccGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGCTAAACCTTTGTTCGCTTTTGGTAATAAATCCCGACAATCATAACGAAAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGTAACAGAGGTCAGCGCCCGTCAGCAGCTTAGCGTTCCGCAACAAAATga                                                                        >  1:530721/1‑249 (MQ=255)
                                                                                                                                                                    gAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTAt                                                                     <  2:158751/251‑1 (MQ=255)
                                                                                                                                                                               cTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACAtt                                                            >  2:294654/1‑249 (MQ=255)
                                                                                                                                                                                tCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGt                                                         >  2:66546/1‑251 (MQ=255)
                                                                                                                                                                                    ggTTGGATAAAATGAATCAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCCTTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGACAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAACATGGTATCTGGACATTAGTTCa                                                      >  2:35909/1‑250 (MQ=255)
                                                                                                                                                                                     attGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCAGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGttaaa                                                     >  2:301690/2‑247 (MQ=255)
                                                                                                                                                                                                       aTGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTATCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATcc                                  <  1:343511/251‑1 (MQ=255)
                                                                                                                                                                                                        tGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGCTAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAACTGGTATCTGGACATTCGTTCAAAGCGAAGAGGCGGAGATcc                                  >  1:34759/1‑250 (MQ=255)
                                                                                                                                                                                                        tGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGAGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATcc                                  >  1:214294/1‑250 (MQ=255)
                                                                                                                                                                                                           aGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAAcc                              >  2:219912/1‑251 (MQ=255)
                                                                                                                                                                                                           aGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAAcc                              <  1:415434/251‑1 (MQ=255)
                                                                                                                                                                                                                                cTTTTGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCAt         >  1:261131/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tGGTAATTAATCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTg     >  1:500777/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aattaatCCCGAAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGt  >  1:563635/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
AATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGATCCCCATAATTACCTTTTTATTCTTGTTTGCGCAAATAACGCATGGCAAAACATTCCTGCCGAGCGGCTTCGCTCATGGTTGGATAAAATGAATAAATGGAGCAGGTTAAACCATTGTTCGCTTTTGGTAATTAATCCCGGAAATAATAACGATAAACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGT  >  NC_000913/3696666‑3697147

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