Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,769 A→G 5.8% L351L (TTA→TTG insQ → putative insertion element transposase InsQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,7690AG 5.8% 212.8 / 2.8 69L351L (TTA→TTGinsQputative insertion element transposase InsQ
Reads supporting (aligned to +/- strand):  ref base A (41/24);  new base G (2/2);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTC  >  NC_000913/1504622‑1504918
                                                                                                                                                   |                                                                                                                                                     
gcttgctGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAgc                                                                                                                                                     >  1:816728/1‑150 (MQ=255)
gcttgctGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAgc                                                                                                                                                     <  1:1052926/150‑1 (MQ=255)
gcttgctCTTACGCCAGAGTATACAATCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAACTAAATACAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAgc                                                                                                                                                     >  2:753356/1‑150 (MQ=255)
   tgctGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAgcggc                                                                                                                                                  >  2:1083036/1‑150 (MQ=255)
          ccGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACAGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCAc                                                                                                                                           >  2:1388233/1‑150 (MQ=255)
                 cGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTc                                                                                                                                    <  2:654232/150‑1 (MQ=255)
                    ataCAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTGGCGGCGGGGCACGCCGTTCTTg                                                                                                                                 >  1:75472/1‑150 (MQ=255)
                      aCAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGTCGGGGCACGCCGTTCTTGcc                                                                                                                               <  2:1355636/150‑1 (MQ=255)
                        aaGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCtg                                                                                                                             >  2:1392416/1‑150 (MQ=255)
                         aGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTg                                                                                                                                 <  1:257777/145‑1 (MQ=255)
                         aGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTg                                                                                                                                 >  2:257777/1‑145 (MQ=255)
                          gCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTGGCGGCGGGGCACGCCGTTCTTGCCtgtg                                                                                                                           >  1:1099360/1‑150 (MQ=255)
                               cGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGgagag                                                                                                                      <  1:1083036/150‑1 (MQ=255)
                               cGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGATCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGACTGTGgagag                                                                                                                      >  1:609348/1‑150 (MQ=255)
                                 ttGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCATTCTTGCCTGTGGAGAGAt                                                                                                                    <  2:1300220/150‑1 (MQ=255)
                                 ttGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCATTCTTGCCTGTGGAGAGAt                                                                                                                    >  1:831576/1‑150 (MQ=255)
                                   gcgcGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCATTCTTGCCTGTGGAGAGATgg                                                                                                                  <  2:831576/150‑1 (MQ=255)
                                      cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGc                                                                                                               <  1:1144141/150‑1 (MQ=255)
                                           tgtgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCa                                                                                                          <  1:720908/150‑1 (MQ=255)
                                            gtgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAg                                                                                                         >  2:247911/1‑150 (MQ=255)
                                             tgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAgg                                                                                                        <  2:1135494/150‑1 (MQ=255)
                                              ggTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTGGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGgc                                                                                                       <  2:75472/150‑1 (MQ=255)
                                               gTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGgcc                                                                                                      >  1:818844/1‑150 (MQ=255)
                                                  aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGgccgcc                                                                                                   >  1:1389645/1‑150 (MQ=255)
                                                       gCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTg                                                                                              <  1:676029/150‑1 (MQ=255)
                                                         gAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGaa                                                                                            >  1:758526/1‑150 (MQ=255)
                                                                 aCCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAAc                                                                                    >  2:924110/1‑150 (MQ=255)
                                                                    gCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCa                                                                                 >  2:542569/1‑150 (MQ=255)
                                                                      cTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCAcc                                                                               <  1:1388233/150‑1 (MQ=255)
                                                                        gTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGa                                                                             >  1:295346/1‑150 (MQ=255)
                                                                        gTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGa                                                                             >  2:866624/1‑150 (MQ=255)
                                                                             aaaGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGa                                                                        >  1:711656/1‑150 (MQ=255)
                                                                              aaGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGAt                                                                       <  1:924110/150‑1 (MQ=255)
                                                                                gTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTc                                                                     <  2:758526/150‑1 (MQ=255)
                                                                                  aaaTTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAg                                                                   <  2:895198/150‑1 (MQ=255)
                                                                                    aTTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGc                                                                 >  1:1137778/1‑150 (MQ=255)
                                                                                     ttCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCg                                                                >  1:1399872/1‑150 (MQ=255)
                                                                                       cAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGAc                                                              >  2:174000/1‑150 (MQ=255)
                                                                                       cAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGAc                                                              >  2:1020604/1‑150 (MQ=255)
                                                                                            gCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCt                                                         >  2:547358/1‑150 (MQ=255)
                                                                                             ccAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTg                                                        <  2:818844/150‑1 (MQ=255)
                                                                                             ccAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTg                                                        >  2:82619/1‑150 (MQ=255)
                                                                                              cAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGa                                                       >  1:34699/1‑150 (MQ=255)
                                                                                                    tgtgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAg                                                 >  1:603571/1‑150 (MQ=255)
                                                                                                    tgtgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAg                                                 <  1:1020604/150‑1 (MQ=255)
                                                                                                     gtgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGc                                                >  1:1012962/1‑150 (MQ=255)
                                                                                                       ggATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAg                                              >  1:214646/1‑150 (MQ=255)
                                                                                                        gATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAgg                                             <  2:816728/150‑1 (MQ=255)
                                                                                                          tataCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGa                                           >  1:938025/1‑150 (MQ=255)
                                                                                                          tataCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGa                                           >  2:736983/1‑150 (MQ=255)
                                                                                                            taCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATc                                         >  2:1388495/1‑150 (MQ=255)
                                                                                                              cAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCa                                       >  1:944567/1‑150 (MQ=255)
                                                                                                              cAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCa                                       >  1:789386/1‑150 (MQ=255)
                                                                                                              cAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCa                                       >  2:565092/1‑150 (MQ=255)
                                                                                                               aGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCAc                                      >  2:1012730/1‑150 (MQ=255)
                                                                                                                gCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACg                                     <  1:452426/150‑1 (MQ=255)
                                                                                                                              aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCAcc                                                                               <  2:474365/94‑1 (MQ=255)
                                                                                                                              aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCAcc                                                                               >  1:474365/1‑94 (MQ=255)
                                                                                                                              aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                       >  1:662976/1‑150 (MQ=255)
                                                                                                                                aTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgga                     <  1:247911/150‑1 (MQ=255)
                                                                                                                                   gcgcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggagga                  >  1:1147329/1‑150 (MQ=255)
                                                                                                                                   gcgcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggatcc                  <  2:857455/150‑4 (MQ=255)
                                                                                                                                           aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAata          >  2:417986/1‑150 (MQ=255)
                                                                                                                                           aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAata          >  2:103064/1‑150 (MQ=255)
                                                                                                                                           aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAata          >  1:1034197/1‑150 (MQ=255)
                                                                                                                                             cATTTTGGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATc        <  2:1099360/150‑1 (MQ=255)
                                                                                                                                              aTTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCg       <  2:295346/150‑1 (MQ=255)
                                                                                                                                               ttttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCgg      <  1:547358/150‑1 (MQ=255)
                                                                                                                                                tttAGCGGCGGGGCACCCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGCACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGt     >  2:497500/1‑150 (MQ=255)
                                                                                                                                                   aGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTc  <  2:945746/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTC  >  NC_000913/1504622‑1504918

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