Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,804,926 (A)6→5 coding (112/1203 nt) proV → glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,804,9210A.98.2% 262.1 / ‑1.5 58coding (107/1203 nt)proVglycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (26/30);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAA  >  NC_000913/2804771‑2805066
                                                                                                                                                      |                                                                                                                                                 
gggTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGa                                                                                                                                                   <  2:15791/151‑1 (MQ=255)
      tatatCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGaaaaact                                                                                                                                             >  1:499641/1‑149 (MQ=255)
                tAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGc                                                                                                                                  >  2:159130/1‑151 (MQ=255)
                tAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGc                                                                                                                                  >  1:59729/1‑151 (MQ=255)
                     aaaTAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGc                                                                                                                             >  2:699123/1‑151 (MQ=255)
                         aaaGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTaa                                                                                                                         <  1:97693/151‑1 (MQ=255)
                                   ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGcc                                                                                                                 >  1:319412/1‑149 (MQ=255)
                                    tCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGt                                                                                                              >  1:80023/1‑151 (MQ=255)
                                      tATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCt                                                                                                            >  1:132289/1‑151 (MQ=255)
                                      tATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCt                                                                                                            >  1:45640/1‑151 (MQ=255)
                                       aTTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCt                                                                                                            <  1:667331/150‑1 (MQ=255)
                                       aTTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTg                                                                                                           >  2:290419/1‑151 (MQ=255)
                                        ttGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTgg                                                                                                          >  1:366127/1‑151 (MQ=255)
                                           cATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGcc                                                                                                        >  1:595901/1‑150 (MQ=255)
                                              ggCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTg                                                                                                    >  2:334074/1‑151 (MQ=255)
                                               gCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTga                                                                                                   <  1:561754/151‑1 (MQ=255)
                                                 aattaaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTga                                                                                                   >  1:422828/1‑149 (MQ=255)
                                                  attaaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTga                                                                                                   >  1:625802/1‑148 (MQ=255)
                                                          aGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCgag                                                                                         <  1:685087/150‑1 (MQ=255)
                                                            aaaTTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGata                                                                                      <  1:681262/151‑1 (MQ=255)
                                                             aaTTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGatat                                                                                     >  2:621201/1‑151 (MQ=255)
                                                                tAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTg                                                                                  >  1:21245/1‑151 (MQ=255)
                                                                 aaaaaTCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTg                                                                                  <  2:675572/150‑1 (MQ=255)
                                                                          tataAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATggg                                                                        >  1:361928/1‑151 (MQ=255)
                                                                          tataAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATggg                                                                        >  2:126233/1‑151 (MQ=255)
                                                                           ataAAATATTTCGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAATAACTTTCAAAAGAACAAATTCT‑GCAAAAACTAAACTATCTCTTGTCGTAAAAGACGCCAGTCTTTCCATTCAAAAAGACAATATATTTGTcatcatttta                                                                       <  1:369247/151‑5 (MQ=255)
                                                                            taAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                      >  1:187945/1‑151 (MQ=255)
                                                                            taAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                      >  2:307437/1‑151 (MQ=255)
                                                                             aaaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                      >  1:589105/1‑151 (MQ=255)
                                                                             aaaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAtt                                                                     <  2:127220/151‑1 (MQ=255)
                                                                               aaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTAt                                                                   <  1:89735/151‑1 (MQ=255)
                                                                               aaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTAt                                                                   >  1:305721/1‑151 (MQ=255)
                                                                                  atTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATcc                                                                 <  1:22114/150‑1 (MQ=255)
                                                                                  atTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCg                                                                >  1:405645/1‑151 (MQ=255)
                                                                                      ggCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCt                                                             <  2:718885/150‑1 (MQ=255)
                                                                                          aGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGt                                                        <  1:138652/151‑1 (MQ=255)
                                                                                           gCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTa                                                       <  1:21580/151‑1 (MQ=255)
                                                                                            cATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTaa                                                      <  2:255969/151‑1 (MQ=255)
                                                                                             aTCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTaa                                                      >  1:195828/1‑150 (MQ=255)
                                                                                               ccACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATc                                                   <  2:687476/151‑1 (MQ=255)
                                                                                                 acaGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCca                                                 <  2:507684/151‑1 (MQ=255)
                                                                                                   agcgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCcaca                                               <  2:77104/151‑1 (MQ=255)
                                                                                                     cgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAAt                                             <  2:289684/151‑1 (MQ=255)
                                                                                                        gcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGACTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGt                                          <  2:588205/151‑1 (MQ=255)
                                                                                                          gTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTAc                                        <  2:106054/151‑1 (MQ=255)
                                                                                                                  atatCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTtct                                 <  1:639214/150‑1 (MQ=255)
                                                                                                                    atCGAACAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCaa                              <  1:728732/151‑1 (MQ=255)
                                                                                                                          cAAGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCt                        <  1:415374/151‑1 (MQ=255)
                                                                                                                           aaGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTg                       <  2:386222/151‑1 (MQ=255)
                                                                                                                           aaGGACTTTCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTg                       <  2:617577/151‑1 (MQ=255)
                                                                                                                                 tttCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGa                  >  1:247682/1‑150 (MQ=255)
                                                                                                                                   tCAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacc               >  2:351557/1‑151 (MQ=255)
                                                                                                                                    cAAAAGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAaccc              <  2:90817/151‑1 (MQ=255)
                                                                                                                                     aaaaGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGTGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAaccc              >  1:455950/1‑150 (MQ=255)
                                                                                                                                     aaaaGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAaccca             <  1:438346/151‑1 (MQ=255)
                                                                                                                                      aaaGAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacccac            >  1:622101/1‑151 (MQ=255)
                                                                                                                                         gAACAAATTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCcg         <  1:266627/151‑1 (MQ=255)
                                                                                                                                               aTTCTGG‑AAAAACTGTGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCAGGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCTCCCGCGGGCa   <  1:520509/151‑1 (MQ=255)
                                                                                                                                               aTTCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCa   <  1:257038/151‑1 (MQ=255)
                                                                                                                                                 tCTGG‑AAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCaa  <  1:103970/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GGGTAATATATCGACATAGACAAATAAAGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAA  >  NC_000913/2804771‑2805066

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