Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,208,205 G→A 9.5% G61D (GGT→GAT)  ttdT → L‑tartrate/succinate antiporte

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,208,2050GA 9.5% 210.3 / 10.1 72G61D (GGT→GAT) ttdTL‑tartrate/succinate antiporte
Reads supporting (aligned to +/- strand):  ref base G (28/37);  new base A (2/5);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

CGCCGCTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTT  >  NC_000913/3208058‑3208353
                                                                                                                                                   |                                                                                                                                                    
cgccgcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTggg                                                                                                                                                      <  1:208202/148‑1 (MQ=255)
cgccgcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGt                                                                                                                                                     >  3:357751/1‑149 (MQ=255)
cgccgcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGt                                                                                                                                                     <  3:627896/149‑1 (MQ=255)
    gcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTatat                                                                                                                                                 <  4:168517/149‑1 (MQ=255)
        gCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCca                                                                                                                                              <  1:111520/148‑1 (MQ=255)
            tCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                          <  1:782108/148‑1 (MQ=255)
              aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGc                                                                                                                                       <  4:124625/149‑1 (MQ=255)
              aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCg                                                                                                                                        <  2:31840/148‑1 (MQ=255)
              aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCg                                                                                                                                        <  2:390913/148‑1 (MQ=255)
              aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCg                                                                                                                                        <  2:800797/148‑1 (MQ=255)
                 gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGCTGGGTATATCCATCATCGCCa                                                                                                                                     <  1:552980/148‑1 (MQ=255)
                  ccATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCAt                                                                                                                                    >  1:106560/1‑148 (MQ=255)
                   cATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                                   >  2:18756/1‑148 (MQ=255)
                   cATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTCATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                                   >  1:653399/1‑148 (MQ=255)
                     ttattGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATAct                                                                                                                                 >  1:188637/1‑148 (MQ=255)
                     ttattGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGCCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCGGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATAct                                                                                                                                 <  1:19636/148‑1 (MQ=255)
                          gCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTcc                                                                                                                           <  4:357751/149‑1 (MQ=255)
                          gCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTcc                                                                                                                           <  4:543987/149‑1 (MQ=255)
                           ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTcc                                                                                                                           <  1:302587/148‑1 (MQ=255)
                           ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTcc                                                                                                                           <  1:18756/148‑1 (MQ=255)
                            tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                         <  3:529252/149‑1 (MQ=255)
                            tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                         <  4:449752/133‑1 (MQ=255)
                            tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                         >  3:449752/1‑133 (MQ=255)
                            tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGATGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                         >  4:50047/1‑149 (MQ=255)
                                     ggTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGctgctg                                                                                                                 <  1:357630/148‑1 (MQ=255)
                                       ttCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtt                                                                                                               >  1:634720/1‑148 (MQ=255)
                                        tCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTc                                                                                                              <  1:696241/148‑1 (MQ=255)
                                        tCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCa                                                                                                             <  4:92960/149‑1 (MQ=255)
                                           cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                          >  3:38355/1‑149 (MQ=255)
                                           cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                          <  3:331337/149‑1 (MQ=255)
                                                    ggAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcagc                                                                                                 <  4:130210/149‑1 (MQ=255)
                                                         aCCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcgct                                                                                            >  4:89525/1‑149 (MQ=255)
                                                             tACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAg                                                                                         >  1:58226/1‑148 (MQ=255)
                                                              aCCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                        >  2:719971/1‑148 (MQ=255)
                                                                 tGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAg                                                                                     <  1:401657/148‑1 (MQ=255)
                                                                  ggCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCGGCCAgg                                                                                    <  2:106560/148‑1 (MQ=255)
                                                                  ggCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAgg                                                                                    <  2:749365/148‑1 (MQ=255)
                                                                     tctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCttt                                                                                >  3:329963/1‑149 (MQ=255)
                                                                      ctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTc                                                                                                              <  2:108305/118‑1 (MQ=255)
                                                                      ctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTc                                                                                                              >  1:108305/1‑118 (MQ=255)
                                                                         cTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                             >  1:102257/1‑148 (MQ=255)
                                                                          tttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGc                                                                                   <  2:276698/141‑1 (MQ=255)
                                                                          tttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGc                                                                                   >  1:276698/1‑141 (MQ=255)
                                                                            tGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                         <  4:38355/149‑1 (MQ=255)
                                                                             gCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCTGGCTTTAAATTCAc                                                                        <  4:777887/149‑1 (MQ=255)
                                                                             gCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                         <  2:634720/148‑1 (MQ=255)
                                                                             gCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                        <  4:430189/149‑1 (MQ=255)
                                                                                 tttttACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCa                                                                    >  4:492255/1‑149 (MQ=255)
                                                                                       cTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATcc                                                              <  3:821970/149‑1 (MQ=255)
                                                                                         ggCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATccc                                                             >  1:893752/1‑148 (MQ=255)
                                                                                           cGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGCGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCctc                                                           >  2:572706/1‑148 (MQ=255)
                                                                                                  gTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGc                                                   >  4:457100/1‑149 (MQ=255)
                                                                                                    tGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGt                                                       >  1:69757/1‑143 (MQ=255)
                                                                                                    tGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGt                                                       <  2:69757/143‑1 (MQ=255)
                                                                                                         tGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTc                                             >  1:766942/1‑148 (MQ=255)
                                                                                                          gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                            <  1:719971/148‑1 (MQ=255)
                                                                                                            tCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCgg                                          >  1:569752/1‑148 (MQ=255)
                                                                                                              cTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGttt                                       >  4:850791/1‑149 (MQ=255)
                                                                                                                cGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttt                                      <  1:53264/148‑1 (MQ=255)
                                                                                                                cGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttt                                      >  2:302918/1‑148 (MQ=255)
                                                                                                                 gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTc                                    <  4:674666/149‑1 (MQ=255)
                                                                                                                  aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTc                                                                          >  1:541026/1‑110 (MQ=255)
                                                                                                                  aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTc                                                                          <  2:541026/110‑1 (MQ=255)
                                                                                                                      cGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGTTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAAtt                               <  3:50047/149‑1 (MQ=255)
                                                                                                                      cGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAAtt                               <  4:778420/149‑1 (MQ=255)
                                                                                                                       gTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTc                              <  3:366862/149‑1 (MQ=255)
                                                                                                                        tGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCg                             <  4:281628/149‑1 (MQ=255)
                                                                                                                         gCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCgg                            <  3:457100/149‑1 (MQ=255)
                                                                                                                            gggTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTa                         >  4:65428/1‑149 (MQ=255)
                                                                                                                              gTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATc                       <  3:89525/149‑1 (MQ=255)
                                                                                                                                 ccGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTgg                    >  3:784539/1‑149 (MQ=255)
                                                                                                                                     tatggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGa                >  3:221473/3‑149 (MQ=255)
                                                                                                                                      gtggCGATGGTGGGTATATCCAGCATCGCCATGCTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAt               >  4:775945/1‑149 (MQ=255)
                                                                                                                                                  ggTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCt    >  2:188159/1‑148 (MQ=255)
                                                                                                                                                   gTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCttt  >  3:897623/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCGCTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTT  >  NC_000913/3208058‑3208353

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