Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,500 T→C intergenic (+145/‑210) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5000TC96.7% 211.6 / NA 61intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/36);  minor base A (2/0);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCA  >  NC_000913/2767365‑2767634
                                                                                                                                          |                                                                                                                                      
ttctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑cccc                                                                                                                                     <  2:567187/139‑1 (MQ=255)
         aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTaa                                                                                                                            <  2:1670581/139‑1 (MQ=255)
                    aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCTCCCCCACAGCATTCTGCAAAATGAAAACAGTTACCCCAGTCTAACTACGATTGAGT‑T‑CACACACCTTTAAACATATCaata                                                                                                                 >  1:1224162/1‑139 (MQ=255)
                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACCACGCCCTGCATAAGGAACACGTTTTTCACTGTGCACCCACGACTAATT‑C‑TCCCCCCCTATAAACTTCTAAATAAGGt                                                                                                             >  1:667613/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCACAAAACGTCCTGCATAGGGAACACGTTTTCCCCTTTGCACCCCCGACTAAAA‑T‑CCCCCCCCTCTAAAcacaccaacaattt                                                                                                             >  1:1618308/1‑126 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATCAGGAACACGTTTGCACCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTc                                                                                                            >  1:1693609/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTATTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTc                                                                                                            >  1:939129/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACCTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCACCCCCTTTAAACTTCTCAATAATGTc                                                                                                            >  1:977375/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCACACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACAAAAT‑T‑ACCCCCCCTTCAAACTTCTCAACAATATc                                                                                                            >  2:1360917/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CACCCCCTCACAACGTCCTGCATTGGGGACACGTTTTCCCCCGTGCACCCACCACTAAATAA‑CCCCCCCCTTTTAACTTCCCAATCATGt                                                                                                             >  1:37350/1‑139 (MQ=255)
                                       atatTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                                                              <  2:1224162/139‑1 (MQ=255)
                                         atTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                                                            <  2:335298/139‑1 (MQ=255)
                                         agagacaggTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                                                            <  2:1365020/131‑1 (MQ=255)
                                           tAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                                                          <  2:868449/139‑1 (MQ=255)
                                                 gTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                                                               <  2:506889/128‑1 (MQ=255)
                                                 gTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAAGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCAACAATTTTCTgg                                                                                               >  1:506889/1‑128 (MQ=255)
                                                  tCCAATATGATCAAATAGCATTAA‑CCCCCTCTGCCCACGTCCTGCATAGGGAACACGTTTTCACCTGTGCACCCACGACAAATT‑T‑ACCCCCCCTTTAATCTTCTCAAATATGTCACCAAGTGTCTGGAAATCGTatgat                                                                                   >  1:109079/1‑139 (MQ=255)
                                                  tCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                                                                      <  2:805932/136‑1 (MQ=255)
                                                  tCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                                                                      >  1:805932/1‑136 (MQ=255)
                                                         tGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                            <  2:875106/139‑1 (MQ=255)
                                                            tCAAATAGCATTAA‑CCCGCCCAAACAACGTCGTTCAGCAGGAACAGGTCTTCCTCTTAGAACCCTCTGATAAAT‑C‑CCCCGCCCTTAGATCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTGTctct                                                                         >  1:1365020/1‑136 (MQ=255)
                                                            tCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACATCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                                                         >  2:699593/1‑139 (MQ=255)
                                                                 tAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGAGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTAGGATGCCTTTCCCGGTaaa                                                                    <  1:992768/139‑1 (MQ=255)
                                                                  aGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                                                                   <  2:939129/139‑1 (MQ=255)
                                                                               cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACGCCCGACTAAAT‑T‑CCCCCCCCCTTTAACGCTTCCATAATGTCACTAATTGTCTGGAAATCCTTTGATGCCTTTCACGATAAACGGTGCtgttcgt                                                       >  1:840380/1‑135 (MQ=255)
                                                                                ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt                                                      <  2:374972/139‑1 (MQ=255)
                                                                                  tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                        <  1:1350530/135‑1 (MQ=255)
                                                                                  tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                        >  2:1350530/1‑135 (MQ=255)
                                                                                     caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaa                                                 >  1:1556902/1‑139 (MQ=255)
                                                                                       aCGTCCTGCATTGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg                                               <  1:1360917/139‑1 (MQ=255)
                                                                                        cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              >  2:185650/1‑139 (MQ=255)
                                                                                        cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              <  2:1556902/139‑1 (MQ=255)
                                                                                        cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTGTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              <  1:311101/139‑1 (MQ=255)
                                                                                          tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                                                       >  2:1441815/1‑96 (MQ=255)
                                                                                          tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                                                       <  1:1441815/96‑1 (MQ=255)
                                                                                          tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt                                            >  2:499023/1‑139 (MQ=255)
                                                                                            cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCt                                          >  2:212804/1‑139 (MQ=255)
                                                                                             tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                         <  2:1618308/139‑1 (MQ=255)
                                                                                              gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        <  1:212804/139‑1 (MQ=255)
                                                                                               cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                            <  2:75199/102‑1 (MQ=255)
                                                                                               cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                            >  1:75199/1‑102 (MQ=255)
                                                                                               cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCa                                       <  1:499023/139‑1 (MQ=255)
                                                                                                aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCaa                                      <  1:699593/139‑1 (MQ=255)
                                                                                                     gAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTAt                                 <  2:840380/139‑1 (MQ=255)
                                                                                                         acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           >  1:1165611/1‑109 (MQ=255)
                                                                                                         acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           <  2:1165611/109‑1 (MQ=255)
                                                                                                             tttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg                                               >  2:75630/1‑117 (MQ=255)
                                                                                                             tttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg                                               <  1:75630/117‑1 (MQ=255)
                                                                                                              ttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           <  1:1310594/104‑1 (MQ=255)
                                                                                                              ttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           >  2:1310594/1‑104 (MQ=255)
                                                                                                              ttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           <  1:498674/104‑1 (MQ=255)
                                                                                                              ttCCCCTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           >  2:498674/1‑104 (MQ=255)
                                                                                                                 cccTGTGCACCCACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGaaaa                     <  2:977375/139‑1 (MQ=255)
                                                                                                                           ccACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                        >  2:1614157/1‑94 (MQ=255)
                                                                                                                           ccACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                        <  1:1614157/94‑1 (MQ=255)
                                                                                                                            ctcGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCTATAATGTTACCAATTTGCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat          <  1:167028/137‑1 (MQ=255)
                                                                                                                            cACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat          <  2:1531856/139‑1 (MQ=255)
                                                                                                                            cACGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat          <  2:1703801/139‑1 (MQ=255)
                                                                                                                             aCGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                                   >  2:396043/1‑97 (MQ=255)
                                                                                                                             aCGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                                   <  1:396043/97‑1 (MQ=255)
                                                                                                                             aCGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagata         >  2:1111809/1‑139 (MQ=255)
                                                                                                                              cGACTAAAT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagatag        >  2:1283514/1‑139 (MQ=255)
                                                                                                                                aCTAAAT‑TCCCCCCCCCTTTTAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga       <  2:37350/139‑1 (MQ=255)
                                                                                                                                    aaT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCa  <  2:109079/139‑1 (MQ=255)
                                                                                                                                    aaT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCa  <  2:667613/139‑1 (MQ=255)
                                                                                                                                    aaT‑T‑CCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCa  <  1:185650/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCA  >  NC_000913/2767365‑2767634

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