Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,295 T→A 9.4% intergenic (+29/‑80) accC → / → yhdT acetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,407,2950TA 9.4% 219.0 / 11.7 64intergenic (+29/‑80)accC/yhdTacetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/4);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTG  >  NC_000913/3407150‑3407437
                                                                                                                                                 |                                                                                                                                              
gACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTcc                                                                                                                                             <  2:359292/149‑1 (MQ=255)
 aCGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCg                                                                                                                                            <  2:22855/149‑1 (MQ=255)
  cGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCgg                                                                                                                                           >  4:148749/1‑149 (MQ=255)
  cGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCgg                                                                                                                                           >  2:1136092/1‑149 (MQ=255)
     tATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCt                                                                                                                                        <  1:519019/149‑1 (MQ=255)
       tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttt                                                                                                                                      <  3:276810/149‑1 (MQ=255)
       tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttt                                                                                                                                      >  1:354876/1‑149 (MQ=255)
       tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                                                                      <  4:848828/149‑1 (MQ=255)
       tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                                                                      <  4:449213/149‑1 (MQ=255)
          aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtttttt                                                                                                                                   >  4:416304/1‑149 (MQ=255)
                  gTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                                           <  4:81234/149‑1 (MQ=255)
                  gTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                                           >  1:579896/1‑149 (MQ=255)
                          gCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCt                                                                                                                                        >  2:115454/1‑128 (MQ=255)
                                ccGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                                                             <  2:354876/149‑1 (MQ=255)
                                         gaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTa                                                                                                    <  1:1136092/149‑1 (MQ=255)
                                              aCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                                           >  1:910795/1‑121 (MQ=255)
                                              aCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                                           <  2:910795/121‑1 (MQ=255)
                                                gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATc                                                                                                                      >  1:1015401/1‑124 (MQ=255)
                                                gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATc                                                                                                                      <  2:1015401/124‑1 (MQ=255)
                                                gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATccc                                                                                             >  4:258823/1‑149 (MQ=255)
                                                  gaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                                                                        <  3:891816/120‑1 (MQ=255)
                                                  gaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                                                                        >  4:891816/1‑120 (MQ=255)
                                                      ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCtt                                                                                        >  1:416430/1‑148 (MQ=255)
                                                      ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCtt                                                                                        <  2:416430/148‑1 (MQ=255)
                                                       tCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttc                                                                                      >  2:728340/1‑149 (MQ=255)
                                                            cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAc                                                                                 >  3:599934/1‑149 (MQ=255)
                                                                            tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaa                                                                 >  3:880686/1‑149 (MQ=255)
                                                                                       gagaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAg                                                                                                       <  2:102872/100‑1 (MQ=255)
                                                                                       gagaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAg                                                                                                       >  1:102872/1‑100 (MQ=255)
                                                                                       gagaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCg                                                      >  3:496154/1‑149 (MQ=255)
                                                                                         gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCa                                                                           <  4:482712/126‑1 (MQ=255)
                                                                                         gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCa                                                                           >  3:482712/1‑126 (MQ=255)
                                                                                          aaaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGttt                                                   <  2:579896/149‑1 (MQ=255)
                                                                                           aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtttt                                                  <  3:148749/149‑1 (MQ=255)
                                                                                              aaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGtt                                               <  3:258823/149‑1 (MQ=255)
                                                                                              aaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGtt                                               >  2:432120/1‑149 (MQ=255)
                                                                                               aCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTc                                              <  3:416304/149‑1 (MQ=255)
                                                                                                           ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAag                                  >  3:846253/1‑149 (MQ=255)
                                                                                                            gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAaga                                 <  1:728340/149‑1 (MQ=255)
                                                                                                                aaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTc                                                       <  3:11619/123‑1 (MQ=255)
                                                                                                                aaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTc                                                       >  4:11619/1‑123 (MQ=255)
                                                                                                                 aTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgcgc                            >  1:94878/1‑149 (MQ=255)
                                                                                                                   aaGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                                                              >  4:914514/1‑113 (MQ=255)
                                                                                                                   aaGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                                                              <  3:914514/113‑1 (MQ=255)
                                                                                                                    aGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    >  2:26752/1‑122 (MQ=255)
                                                                                                                    aGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    <  1:26752/122‑1 (MQ=255)
                                                                                                                      aCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgg                       <  1:432120/149‑1 (MQ=255)
                                                                                                                        tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    >  3:507576/1‑118 (MQ=255)
                                                                                                                        tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    >  1:1091435/1‑118 (MQ=255)
                                                                                                                        tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    <  4:507576/118‑1 (MQ=255)
                                                                                                                        tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                    <  2:1091435/118‑1 (MQ=255)
                                                                                                                                 gTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCt            >  3:1057463/1‑149 (MQ=255)
                                                                                                                                  tCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttct                                                                                     >  3:148875/1‑75 (MQ=255)
                                                                                                                                  tCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttct                                                                                     <  4:148875/75‑1 (MQ=255)
                                                                                                                                   cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGa          <  4:846253/149‑1 (MQ=255)
                                                                                                                                     aaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGa                                                                       <  4:342258/86‑1 (MQ=39)
                                                                                                                                     aaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGa                                                                       >  3:342258/1‑86 (MQ=38)
                                                                                                                                       aGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCt      >  2:139526/1‑149 (MQ=255)
                                                                                                                                          ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt   <  4:880686/149‑1 (MQ=255)
                                                                                                                                          ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt   >  3:737349/1‑149 (MQ=255)
                                                                                                                                           cGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgct                  <  4:264699/133‑1 (MQ=255)
                                                                                                                                           cGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgct                  >  3:264699/1‑133 (MQ=255)
                                                                                                                                           cGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTg  <  4:737349/149‑1 (MQ=255)
                                                                                                                                           cGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTg  >  2:1117553/1‑149 (MQ=255)
                                                                                                                                               ttttCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGAc         >  1:970713/1‑138 (MQ=255)
                                                                                                                                               ttttCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGAc         <  2:970713/138‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTG  >  NC_000913/3407150‑3407437

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