Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,786 T→A 8.8% E166D (GAA→GAT) ‡ feaR ← transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,7860TA 8.8% 190.4 / 5.7 72E166D (GAA→GAT) ‡feaRtranscriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base T (22/43);  new base A (5/2);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 9.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGA  >  NC_000913/1446638‑1446932
                                                                                                                                                    |                                                                                                                                                  
ggCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGttt                                                                                                                                                   >  2:1317802/1‑150 (MQ=255)
ggCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGACTGTAGCGCAGCTTCACTTTCTGttt                                                                                                                                                   >  1:4342/1‑150 (MQ=255)
 gCGTAATATCTATTCGCGAATATTATCGGCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACTCCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTc                                                                                                                                                  <  2:253671/150‑1 (MQ=255)
  cGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCa                                                                                                                                                 <  2:466330/150‑1 (MQ=255)
     aaTATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGaa                                                                                                                                              <  1:618943/150‑1 (MQ=255)
      atatCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAAATGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTACACTTTCTGTTTCAGaaa                                                                                                                                             >  1:407735/1‑150 (MQ=255)
       tatCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAg                                                                                                                                            >  1:429177/1‑150 (MQ=255)
         tctcTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTg                                                                                                                                          <  2:382801/150‑1 (MQ=255)
         tctcTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTg                                                                                                                                          <  1:1221782/150‑1 (MQ=255)
            ctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCg                                                                                                                                       >  1:1081722/1‑150 (MQ=255)
             ttCGCGAATATTATCGTCTATCAACTTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaagc                                                                                                                                      <  2:112089/150‑5 (MQ=255)
                 cgAATATTATCCTCTATCCACGTAACCACTTTTTGCAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCTGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGatta                                                                                                                                  <  2:79416/150‑1 (MQ=255)
                   aaTATTATGCTATATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGattatt                                                                                                                                <  2:123328/150‑1 (MQ=255)
                        tatCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGc                                                                                                                           <  1:1155361/150‑1 (MQ=255)
                               tATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTg                                                                                                                    >  1:1202116/1‑150 (MQ=255)
                                aTCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGt                                                                                                                   <  1:1084574/150‑1 (MQ=255)
                                  cAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTCCCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTaa                                                                                                                 <  2:355860/150‑1 (MQ=255)
                                    aCGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACa                                                                                                               >  2:899997/1‑150 (MQ=255)
                                      gTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCTCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAgg                                                                                                             <  1:344358/150‑1 (MQ=255)
                                       tAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGTCTCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTACAGAAAGTGCCGGATTATTCATGCTCTCCAGTAACAGGc                                                                                                            >  1:462613/1‑150 (MQ=255)
                                        aaCCACTTTTTGAAACTGACGTTCACGACGAGGGTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCg                                                                                                           <  2:649236/150‑1 (MQ=255)
                                           cACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATg                                                                                                        >  2:583277/1‑150 (MQ=255)
                                             cTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTATGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGAc                                                                                                      >  1:881816/1‑150 (MQ=255)
                                             atttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGAc                                                                                                      <  1:613345/149‑1 (MQ=255)
                                                 ttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAg                                                                                                  <  2:462992/150‑1 (MQ=255)
                                                    aaaCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTg                                                                                               <  2:779549/150‑1 (MQ=255)
                                                     aaCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGATAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGc                                                                                              >  1:607467/1‑150 (MQ=255)
                                                      aCTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCa                                                                                             <  1:1316017/150‑1 (MQ=255)
                                                       cTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCAc                                                                                            <  2:498326/150‑1 (MQ=255)
                                                           cGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                        <  2:873853/150‑1 (MQ=255)
                                                           cGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                        <  2:1017738/150‑1 (MQ=255)
                                                             ttCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATggg                                                                                      <  2:1258615/150‑1 (MQ=255)
                                                              tCACGACGAGGTTGAACAGATTCCCGCGGATGAAGTTCCGGGCTCAGCAGACACACCATCGCCGGTAGCGCAGCGTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGt                                                                                     <  1:176127/150‑1 (MQ=255)
                                                                  gacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGCTTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGt                                                                                 >  1:627782/1‑150 (MQ=255)
                                                                   acgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTc                                                                                <  1:1041840/150‑1 (MQ=255)
                                                                   acgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTc                                                                                <  1:1047243/150‑1 (MQ=255)
                                                                   acgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTc                                                                                <  1:97335/150‑1 (MQ=255)
                                                                     gaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAg                                                                              <  1:957972/150‑1 (MQ=255)
                                                                     gaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAg                                                                              >  1:617096/1‑150 (MQ=255)
                                                                      aGGTTGAACATATTCCCGCTGATGAAGTACCGGGCGCAGCAGACTCAGCATAGCCCGTAGCGCAGCTTCACTTTCTGATTCAGAAACTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGc                                                                             <  1:100864/150‑1 (MQ=255)
                                                                      aGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCTGCAGACAAACCATCGCATGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGATTTTCACAATGGGTAATTCAGc                                                                             >  1:354410/1‑150 (MQ=255)
                                                                      aGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGc                                                                             <  1:76811/150‑1 (MQ=255)
                                                                       ggTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCg                                                                            >  1:89663/1‑150 (MQ=255)
                                                                          tGAACAGATTCCAGCTGATCAAGTACCGGGCGCAGCAGACACACCAACGCCTCTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTCCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGAAGCACCATGGGTAAGTAAGCGTcc                                                                         >  2:600120/1‑150 (MQ=255)
                                                                            aaCAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCGCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAg                                                                       >  1:961334/1‑150 (MQ=255)
                                                                              cAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTc                                                                     <  2:369308/150‑1 (MQ=255)
                                                                                  ttCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTc                                                                 >  2:414887/1‑150 (MQ=255)
                                                                                        cTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCa                                                           <  2:1131605/150‑1 (MQ=255)
                                                                                         tgatgaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGATTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAg                                                          <  1:1135088/150‑1 (MQ=255)
                                                                                               aGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGTTTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGt                                                    >  2:542342/1‑150 (MQ=255)
                                                                                                       gcgcCGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATg                                            <  2:173764/150‑1 (MQ=255)
                                                                                                         gcagcaGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATggg                                          <  1:1111609/150‑1 (MQ=255)
                                                                                                           agcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGa                                        <  2:131163/150‑1 (MQ=255)
                                                                                                           agcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGa                                        <  2:493081/150‑1 (MQ=255)
                                                                                                           agcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTATCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACAATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGa                                        >  1:269312/1‑150 (MQ=255)
                                                                                                                   cacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTg                                >  1:580050/1‑150 (MQ=255)
                                                                                                                       aTCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTcc                            <  2:917661/150‑1 (MQ=255)
                                                                                                                        tCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCca                           >  2:416262/1‑150 (MQ=255)
                                                                                                                           ccTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagca                        <  2:827247/150‑1 (MQ=255)
                                                                                                                            cTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagcag                       >  1:258046/1‑150 (MQ=255)
                                                                                                                             tGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCAAGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTACTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCaga                      >  1:981416/1‑150 (MQ=255)
                                                                                                                              gTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGTGTTGCACCATGGGTAAGTCAGCGTACAGTATTTCTTCGCAGATAGGTTTTTGTTGGGGAAAATATTGTTCCCGCagag                     >  2:153635/1‑150 (MQ=255)
                                                                                                                              gTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCagag                     <  2:1011503/150‑1 (MQ=255)
                                                                                                                              gTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCagag                     <  1:466728/150‑1 (MQ=255)
                                                                                                                                    atccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGGGTg               >  1:715491/7‑150 (MQ=255)
                                                                                                                                     aGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTgg              <  2:761426/150‑1 (MQ=255)
                                                                                                                                      gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc             <  2:1316044/150‑1 (MQ=255)
                                                                                                                                      gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTAATCATGCTCTCCTGTAACAGGCGATGACTGATTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc             <  2:558987/150‑1 (MQ=255)
                                                                                                                                       cggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa            >  2:443635/4‑150 (MQ=255)
                                                                                                                                       cggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa            >  1:496839/4‑150 (MQ=255)
                                                                                                                                       cTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa            <  1:1125932/150‑1 (MQ=255)
                                                                                                                                        ggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaa           >  1:547782/3‑150 (MQ=255)
                                                                                                                                         tCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCGTGTAACTGGCGATGACTGAGTTGCACCATGGGTAAGGCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaaa          <  2:607467/150‑1 (MQ=255)
                                                                                                                                                 ggTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGAAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGGCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGa  <  2:1193084/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGA  >  NC_000913/1446638‑1446932

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