Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 770,824 T→G intergenic (+213/‑634) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,8240TG98.4% 226.5 / ‑3.3 62intergenic (+213/‑634)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (26/35);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGG  >  NC_000913/770676‑770968
                                                                                                                                                    |                                                                                                                                                
ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGaa                                                                                                                                                <  1:135139/151‑1 (MQ=255)
 cGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAAc                                                                                                                                               <  1:63739/151‑1 (MQ=255)
  ggATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGaa                                                                                                                                                >  2:583378/1‑149 (MQ=255)
   gATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACa                                                                                                                                              >  1:182144/1‑150 (MQ=255)
   gATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACa                                                                                                                                              >  1:541437/1‑150 (MQ=255)
    aTACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACaaa                                                                                                                                            <  1:464077/151‑1 (MQ=255)
       cGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAaata                                                                                                                                         <  2:626557/151‑1 (MQ=255)
       cGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAaata                                                                                                                                         <  2:52722/151‑1 (MQ=255)
             gAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAAc                                                                                                                                    <  1:735574/150‑1 (MQ=255)
              aaCGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAAccc                                                                                                                                  <  1:51527/151‑1 (MQ=255)
                 gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTgg                                                                                                                               >  2:314762/1‑151 (MQ=255)
                 gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTgg                                                                                                                               <  2:271869/151‑1 (MQ=255)
                    ttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGt                                                                                                                             <  1:6610/150‑1 (MQ=255)
                      aTTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAg                                                                                                                          <  2:448254/151‑1 (MQ=255)
                        tCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGt                                                                                                                         <  2:210868/150‑1 (MQ=255)
                        tCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGt                                                                                                                         <  1:742992/150‑1 (MQ=255)
                         cGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGtt                                                                                                                        <  2:864464/150‑1 (MQ=255)
                               aCAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATtata                                                                                                                 <  2:528088/151‑1 (MQ=255)
                                  aaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATtataata                                                                                                              <  1:734628/151‑1 (MQ=255)
                                   aaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATtataatat                                                                                                             >  2:192525/1‑151 (MQ=255)
                                          cAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATtata                                                                                                      <  2:828060/151‑1 (MQ=255)
                                          cAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATtata                                                                                                      <  2:529124/151‑1 (MQ=255)
                                             aTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGt                                                                                                   >  2:344961/1‑151 (MQ=255)
                                             aTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGt                                                                                                   <  2:383814/151‑1 (MQ=255)
                                                    atatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTa                                                                                             <  2:784171/150‑1 (MQ=255)
                                                    atatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTa                                                                                             <  2:632328/150‑1 (MQ=255)
                                                         gCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACCCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGt                                                                                        <  1:61682/150‑1 (MQ=255)
                                                          cAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTa                                                                                      <  1:737481/151‑1 (MQ=255)
                                                            gagaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAat                                                                                    >  2:245583/1‑151 (MQ=255)
                                                             agaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAat                                                                                    >  2:259768/1‑150 (MQ=255)
                                                              gatgatGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAatat                                                                                  <  2:473313/151‑1 (MQ=255)
                                                               atgatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATAtt                                                                                 <  2:502313/151‑1 (MQ=255)
                                                                 gatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTgg                                                                               >  1:879735/1‑151 (MQ=255)
                                                                        acacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCgggggg                                                                        >  2:778539/1‑149 (MQ=255)
                                                                        acacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTg                                                                        >  2:182896/1‑151 (MQ=255)
                                                                                        cgcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTTTGCCTttatt                                                        >  1:763861/1‑151 (MQ=255)
                                                                                        cgcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTATGCCTttatt                                                        >  2:314887/1‑151 (MQ=255)
                                                                                          cATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGGCtttttt                                                        >  1:639250/1‑149 (MQ=255)
                                                                                          cATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTttatta                                                       <  1:37325/150‑1 (MQ=255)
                                                                                          cATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAg                                                      <  1:851429/151‑1 (MQ=255)
                                                                                          cATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAg                                                      <  2:591877/151‑1 (MQ=255)
                                                                                           aTCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTATGCCTTTATTAGt                                                     >  2:82522/1‑151 (MQ=255)
                                                                                              aGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTTTTAGTaa                                                   >  1:368388/1‑150 (MQ=255)
                                                                                               ggTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTTTGCCTttattaat                                                     >  1:211557/1‑145 (MQ=255)
                                                                                               ggTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAAt                                                  >  2:283944/1‑150 (MQ=255)
                                                                                               ggTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTATGCCTTTATTAGTaa                                                   >  1:118964/1‑149 (MQ=255)
                                                                                                     tttGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGGTTTATGCCTTTATTATTGATCCTGaa                                            >  2:136181/1‑150 (MQ=255)
                                                                                                       tGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATAGTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGaaa                                           <  2:700508/149‑1 (MQ=255)
                                                                                                        gCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGGTTTATTCCTTTATTAGTAATCCCGaa                                            >  2:835523/1‑147 (MQ=255)
                                                                                                        gCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAAct                                         <  2:719915/150‑1 (MQ=255)
                                                                                                          gTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTg                                      <  2:224202/151‑1 (MQ=255)
                                                                                                                cTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTATTAATCCTGAAACTCTGc                                     >  1:50911/1‑146 (MQ=255)
                                                                                                                      tctcAAGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATTCCTTTATTAGTAATCCTGAAAATCTGCGTCGTATTTGc                          >  1:260286/1‑151 (MQ=255)
                                                                                                                        tcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTATGCCTTTATTTGTAATCCTGAAAATCTGCGGCGTATTAGcc                         >  1:896980/1‑150 (MQ=255)
                                                                                                                           aGCCACGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTg                     <  1:290869/151‑1 (MQ=255)
                                                                                                                                cGTATATGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTTTTTGTAATCCTGAAACTCTTCGTCGTATTTGCCAGTGAcc                  >  2:544757/1‑149 (MQ=255)
                                                                                                                                    tatGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTATGCCTTTTTTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCaaaaa             >  1:452994/1‑150 (MQ=255)
                                                                                                                                     atGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCaaaaaaa           <  2:325332/151‑1 (MQ=255)
                                                                                                                                     atGTGGCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCaaaaaa            <  2:355909/150‑1 (MQ=255)
                                                                                                                                         ggCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGaa        <  1:904517/150‑1 (MQ=255)
                                                                                                                                         ggCTTTATTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAAt       <  1:463657/151‑1 (MQ=255)
                                                                                                                                               aTTTTGAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGGCTTTATTAGTAATCCTGAAACTCTGCGTCGGATTAGCCAGTGACCAAAAAAAAAATTAAgg  >  1:346805/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGG  >  NC_000913/770676‑770968

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